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AVANT-PROPOS 


Ce document s'adresse aux étudiants de la formation d'ingénieur, de Licence et de Master dans le 
cadre des programmes officiels. Mais bien entendu il peut être étudié par tous ceux en 1° cycle, en 2ème cycle, 
ou même en post-graduation, qui désirent approfondir leurs connaissances ou avoir un document de base en 
matière d’asservissement. 


La commande et l'interprétation du comportement de procédés industriels ou de phénomènes 
physiques naturels font partie des tâches qui incombent à l'ingénieur. Ce dernier est confronté à une réalité 
qu'il lui faudra domestiquer et/ou comprendre pour en tirer le meilleur parti. Au centre de cette connaissance se 
trouve le concept de système, concept que l’on retrouve dans un grand nombre de disciplines et techniques : 
contrôle de procédé, techniques d'optimisation, traitement du signal, filtrage, mathématique des équations 
différentielles, etc. 


Dans le cadre de ce cours, nous nous intéressons principalement à l'étude des "systèmes" à la fois 
continus et linéaires, qui sont représentés sous forme de fonction de transfert (représentation externe, dite 
encore de la “boîte noire") ; ces trois conditions volontairement limitatives permettent d'introduire de façon 
simple les principaux concepts de l’automatique. 


Ce texte constitue la transcription fidèle d'un cours oral magistral traité au département 
d'Electrotechnique, à l’université de Sidi Bel-Abbès pour les étudiants en formation de Licence et Master 
(différents parcours). Le but recherché n’est donc pas d'épuiser le sujet, mais d'essayer d'en dégager les idées 
essentielles, simplifiées, quand cela est nécessaire, dans un but didactique. En conséquence, l'accent est mis 
sur les explications physiques et les exemples, plutôt que sur les démonstrations, mais celles-ci sont 
également traitées en détail surtout lorsqu'elles sont indispensables à la bonne compréhension du résultat. 


Pr. FELLAH Mohammed-Karim 
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Chapitre 1 : INTRODUCTION AUX ASSERVISSEMENTS 


1- 1 - Introduction à l’automatique 


L’automatique est généralement définie comme la science qui traite des ensembles qui se suffisent à 
eux-mêmes et où l'intervention humaine est limitée à l'alimentation en énergie et en matière première. 


L'objectif de l’automatique est de remplacer l’homme dans la plupart des tâches (tâches répétitives, 
pénibles, dangereuses, trop précises, trop rapides) qu'il réalise dans tous les domaines sans intervention 
humaine. 


Les systèmes automatiques permettent donc : 


de réaliser des opérations trop complexes ou délicates ne pouvant être confiés à l'homme, 
de se substituer à l'opérateur pour des tâches répétitives, 
d'accroître la précision, 


d'améliorer la stabilité d'un système et sa rapidité. 


De tels dispositifs se rencontrent fréquemment dans la vie courante, depuis les mécanismes 
biologiques du corps humain jusqu'aux usines entièrement automatisées. 


Une telle science englobe un grand nombre de disciplines et, par conséquent, un automaticien devrait 
être à la fois : 


Mathématicien 
Electricien 
Mécanicien 


Economiste 


1- 1.1 - Exemple 


Nous sommes entourés d'un grand nombre de systèmes automatiques, machine à laver, ascenseur, 
distributeur de boisson, robot, suivi de trajectoire d’un missile. 


1- 1.2 - Classification 


Le domaine des applications de l’automatique est très vaste et varié, mais l'observation de l’industrie 
contemporaine conduit à une certaine classification qui se résume en deux grandes familles selon les données 
que traitent ces systèmes : 


* Les automatismes séquentiels 


* Les asservissements 


Ces deux parties de l’automatique sont nettement différentes, elles s'appuient sur des notions 
théoriques qui n’ont que de lointains rapports entre elles et les techniques qui permettent de les réaliser sont, 
aussi, très différentes. 
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1- 1.2.a - Les automatismes séquentiels 


C'est la branche de l’automatique qui organise le déroulement des différentes opérations relatives au 
fonctionnement d'un ensemble complexe. 


Un automatisme à séquence impose l’ordre dans lequel les opérations se déroulent, s'assure que 
chaque opération est bien terminée avant d'aborder la suivante, décide de la marche à suivre en cas 
d'incidents. 


` 


Bien entendu, un automatisme séquentiel peut avoir à contrôler des asservissements et des 
régulateurs (voir § 1- 1.2.b) parmi les ensembles qu'il gère. 


Ce type ďd’'automatisme est utilisé par exemple dans la mise en route et l'arrêt d'installations complexes 
(centrales automatiques), sur les machines outils et, en général, dans presque toutes unités de production 
automatisées. 


Il faut noter également que toutes les séquences d'alarme et de sécurité industrielle font partie des 
applications de ce type d'automatisme. 


Les automatismes sont des systèmes logiques qui ne traitent que des données logiques (0/1, vrai/faux, 
marche/arrêt,...). Ils utilisent les moyens de commutation offerts par l'électronique (circuit logique) et la 
mécanique (logique pneumatique). Le calcul de ces automatismes impose de connaître l'algèbre de Boole et la 
théorie des circuits séquentiels. 


Ils sont classés en 2 branches : 


* 


Systèmes combinatoires : les sorties du système ne dépendent que des variables d'entrées. 


* 


Systèmes séquentiels : les sorties dépendent bien sûr de l’évolution des entrées mais aussi de 
l’état précédent des sorties. 


Exemple : Machine à laver, manipulateur pneumatique, ascenseur, distributeur de boissons. 


1- 1.2.b - Les asservissements 


Un système asservi est un système qui prend en compte, durant son fonctionnement, l’évolution de ses 
sorties pour les modifier et les maintenir conforme à une consigne. 


Cette branche de l’automatique se décompose en deux autres sous branches (séparées artificiellement 
par l'usage) : 


* Régulation : maintenir une variable déterminée, constante et égale à une valeur, dite de 
consigne, sans intervention humaine. Exemple : Régulation de température d'une pièce. 


Systèmes asservis : faire varier une grandeur déterminée suivant une loi imposée par un 
élément de comparaison. Exemple : Régulation de la vitesse d'un moteur, Suivi de trajectoire d’un 
missile. 


L'asservissement est essentiellement analogique et utilise la partie analogique des trois moyens de 
base dont on dispose : mécanique, électrotechnique et électronique. La théorie des asservissements nécessite 
une bonne base mathématique classique. 
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1- 1.3 - Systèmes continus et invariants 


* _ Système continu : un système est dit continu lorsque les variations des grandeurs physiques le 
caractérisant sont des fonctions du type f(t), avec t une variable continue, le temps en général. On 
oppose les systèmes continus aux systèmes discrets (ou échantillonnés), par exemple les 
systèmes informatiques. 


* _ Système invariant : On dit qu'un système est invariant lorsque les caractéristiques de 
comportement ne se modifient pas avec le temps. 


1- 1.4 - Evolution de l'automatique 


Ces dernières années, l'’automatique s’est considérablement modernisée, surtout depuis l'avènement 
des calculateurs numériques. Les systèmes automatiques conduits par calculateurs assurent la quasi-totalité 
des tâches : 


* ils collectent et traitent les informations issues des capteurs qui fournissent l'ensemble des 
variables d'entrée. 


* ces variables d'entrée constituent les données sur lesquelles des calculs numériques seront 
effectués. Ils correspondent à la résolution numérique de systèmes d'équations qui constituent le 
"modèle mathématique". 


* le résultat de ce traitement fourni en binaire est converti en variables continues et est injecté dans 
le processus, afin de modifier son évolution dans un sens désiré. 


En plus de ces tâches qui sont classiques en automatique, le calculateur joue un rôle optimalisateur. 
C'est-à-dire qu'il exécute le travail à faire aux meilleures conditions économiques en minimisant les déchets, en 
tenant compte du carnet de commande, etc. Cet aspect, lui, est nouveau. Ce genre de problème était traité 
séparément. Ce procédé permet de tenir compte d’un nombre considérable de variables, donc de traiter des 
problèmes jusqu'alors impossibles. En plus, il fait intervenir directement les variables économiques au niveau 
de chaque organe (moteur, pompe, etc ...). Or, jusqu'à présent, les variables économiques n'intervenaient que 
globalement. Il permet donc de traiter ce problème de façon beaucoup plus rationnelle. 


Les systèmes automatiques conduits par calculateurs nécessitent une bonne connaissance de la 
programmation en langage machine, de fortes connaissances mathématiques (pour élaborer le modèle) et 
surtout une connaissance parfaite du processus à réguler, ce qui est le plus délicat. Ceci nécessite encore de 
bonnes connaissances en théorie de l'information, en statistique et en recherche opérationnelle. 


1- 2 - Boucle de régulation 
1-2.1- Notion d'asservissement 


L'objectif d'un système automatisé est de remplacer l'homme dans une tâche donnée. Nous allons, 
pour établir la structure d’un système automatisé, commencer par étudier le fonctionnement d’un système dans 
lequel l'homme est la " partie commande ". 


Exemple : conducteur au volant d’un véhicule 


Le conducteur doit suivre la route. Pour cela, Il observe la route et son 
environnement et évalue la distance qui sépare son véhicule du bord de la route. Il 
détermine, en fonction du contexte, langle qu'il doit donner au volant pour suivre la route. Il 
agit sur le volant (donc sur le système) ; puis de nouveau, il recommence son observation 
pendant toute la durée du déplacement. Si un coup de vent dévie le véhicule, après avoir 
observé et mesuré l'écart, il agit pour s'opposer à cette perturbation. 


Si l'on veut qu'un asservissement remplace l'homme dans diverses tâches, il devra avoir un 
comportement et des organes analogues à ceux d'un être humain. C'est-à-dire qu'il devra être capable 
d'apprécier, de comparer et d'agir. 
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Exemple : ouverture de porte pour accès à une maison. 


Un autre exemple d’asservissement très simple est celui d'un homme qui veut entrer 
| dans une maison : à chaque instant, ses yeux "mesurent" l'écart qui existe entre sa position 
| et la porte. Son cerveau commande alors aux jambes d'agir, en sorte que cet écart diminue, 
| puis s’annule. 


Les yeux jouent alors le rôle d'organes de mesure (ou de capteurs), le cerveau celui de comparateur 
et les jambes celui d'organe de puissance. 


Tout asservissement comportera ces trois catégories d'éléments qui remplissent les 3 grandes 
fonctions nécessaires à sa bonne marche (fig. 1—1) : 


* Mesure (ou observation) 


* Comparaison entre le but à atteindre et la position actuelle (Réflexion) 


Observation 


* 


Action de puissance 


Tâche à réaliser Réflexion Effet de l’action 


Fig. 1—1 : Concept général d’un asservissement 


1- 2.2 - Systèmes bouclés et non bouclés 
1- 2.2.a - Exemple 1 : Tir au canon 


Pour mieux saisir la notion de système bouclé, prenons un exemple avec 2 cas. Dans le premier, nous 
considérons un système non bouclé et nous mettrons en évidence ses faiblesses. Dans le second, nous 
montrerons les avantages qu'apporte le bouclage. 


Premier cas : tir au canon sur une cible. 


On considère une cible à détruire et un canon. Pour atteindre le but que l’on s’est 
proposé, on règle langle de tir du canon et la charge de poudre de lobus en fonction des 
coordonnées de la cible et d'autres paramètres connus à l'instant du tir. Une fois lobus parti, 

| si ces paramètres extérieurs viennent à changer, par exemple si la cible se déplace, on ne 
| peut plus agir sur sa direction : l'obus est abandonné à lui-même. 


Deuxième cas : tir au canon sur une cible avec une fusée téléguidée et un radar. 


Considérons la même cible et une fusée téléguidée. Dans ce cas, même si la cible se 
déplace ou un vent latéral fait dévier la fusée de sa trajectoire initiale, elle atteindra quand 
même son but. En effet, à chaque instant, un radar donnera les positions respectives de la 
fusée et de la cible. Il suffira de les comparer pour en déduire l'erreur de trajectoire et agir sur 
les gouvernes de la fusée pour rectifier cette erreur. Dans ce cas, le système n'est plus 
abandonné à lui-même car il comporte une boucle de retour qui est constituée par le radar, 
qui "mesure" la position de la fusée et qui en informe l'opérateur, et par une télétransmission 
qui permet de modifier la trajectoire par action sur les gouvernes. 


| La boucle de retour apporte donc, au prix d'une complication certaine, un gain de précision 
| énorme. 
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1- 2.2.b - Exemple 2 : Asservissement de vitesse d’une voiture 


Supposons que l’on veuille maintenir constante la vitesse (V) d'une voiture. A la valeur (V) de la vitesse 
correspond une valeur (e) de la course de l'accélérateur. Il suffirait donc, en principe, de maintenir (e) constant 
pour que (V) le soit. Chacun sait que la réalité est différente. 


En effet, le vent, les variations de pente et le mauvais état de la route modifient (V). Ces paramètres 
extérieurs qui influent sur la vitesse sont appelés grandeurs perturbatrices ou perturbations. Si elles 
n'existaient pas, la boucle de régulation serait inutile. 


Pour que la vitesse reste constante, il faut utiliser un tachymètre qui mesure la vitesse réelle. Le 
chauffeur compare à tout instant cette vitesse réelle et la vitesse prescrite; Il en déduit un écart plus ou moins 
grand et enfonce plus ou moins l'accélérateur en fonction de cet écart. 


Si on appelle grandeur de sortie (ou sortie) la vitesse réelle et grandeur d'entrée (ou entrée) la vitesse 
imposée, le chauffeur et le tachymètre assurent une liaison entre l'entrée et la sortie, ils constituent donc une 
chaîne de retour. 


On peut donner un schéma très simple pour illustrer cet exemple (fig. 1—2) : 


aaea De HO- 


Entrée Sortie 
vitesse , vitesse (V) 
imposée Perturbations réelle 


Tachymètre 


Fig. 1-2 : Exemple dď’'asservissement de vitesse d’un véhicule 
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1- 2.3 - Définitions - Constitutions élémentaires 


On peut donc définir un asservissement comme un système bouclé ou à boucle fermée comportant 
une amplification de puissance, une mesure et une comparaison. 


A partir de ces 3 notions, on peut définir un schéma fonctionnel valable pour tous les systèmes 
présentant ces caractéristiques (fig. 1—3) : 


* 


Le triangle : représente la fonction amplification de puissance. 


* 


Le cercle : représente la fonction comparaison (qui s'effectue en faisant une différence). 


* 


Le rectangle : représente la fonction mesure et transformation. 


Fonction comparaison 


Fonction amplification de puissance 


Fonction mesure et 
transformation 


Fig. 1-3 : Schéma fonctionnel 
d'un asservissement 


- La sortie régulée représente le phénomène physique que doit régler le système, 
Grandeur c'est la raison d’être du système. Il peut s’agir d’une tension, d'un déplacement, d’un angle 
de rotation, d’un niveau, d'une vitesse, etc... 


de sortie 

Ẹ 
Grandeur La consigne, est l'entrée d'action, c'est la grandeur réglante du système. Sa 
d'entrée nature peut être différente de celle de (S). Seule importe sa valeur numérique. Si (E) et (S) 


REA sont de natures différentes, il suffit de définir une correspondance numérique entre ces 
ou référence | deux grandeurs. Par exemple, on dira qu'un volt à l'entrée représente 100 tours/mn. 


ou consigne 


On appelle écart ou erreur, la différence entre la consigne et la sortie. Cette 
€ mesure ne peut être réalisée que sur des grandeurs comparables, on la réalisera donc en 
général entre la consigne et la mesure de la sortie. Elle est fournie par le comparateur et 
est proportionnelle à la différence ( E-S’ ). Elle peut être de nature différente. Par 


exemple, E et S’ étant des tensions, on pourra avoir € sous forme de courant tel que £ = ( 
E-S')/R_ (R est une résistance). 


erreur 


ou écart 
entrée - sortie 


S’ Elle est fournie par la chaîne de retour, généralement après transformation. S’ doit 
obligatoirement avoir même nature physique que E. Ce qui est évident si on veut donner 
un sens à la différence ( E - Ş ). Un des rôles de la chaîne de retour est donc d'assurer la 
de la sortie conversion de la mesure de S dans la grandeur physique de E. 


Mesure 
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D'une manière générale, le système comprend (fig. 1—4) : 


Perturbations 
éventuelles 


Erreur ou 
Ecart 


\ Chaîne directe (ou d'action) 


l 
l 
Si 
l 
1 
1 
1 


Régulateur \ 
j LE m m E L | = SR LL | an —_ Li EER 


; l 
Entrée de E AA Actionneur , Sortie 
référence Correcteur + asservie 
; processus 
(consigne) — = æ e e e e e 
1 
! Mesure 
UE Capteur 


Comparateur : 


Chaîne de retour (ou d'observation) 


Fig. 1—4 : Organisation fonctionnelle d'un système asservi (schéma fonctionnel) 


Englobe tous les organes de puissance (nécessitant un apport extérieur 
. | d'énergie) et qui exécute le travail. 

Chaîne directe ne iii 7 
Comporte généralement nombreux éléments, notamment des amplificateurs. 


ou d'action E , a | 
* La nature de ces éléments n'est pas spécifiée sur le schéma, il peut s'agir 
aussi bien d'engins électriques, mécaniques, pneumatiques, etc... 
* Analyse et mesure le travail effectué et transmet au comparateur une 
Chaîne de retour grandeur physique proportionnelle à ce travail. 


ou de réaction Elle comprend généralement un capteur qui donne une mesure de la 


grandeur S, qui est ensuite amplifiée et transformée avant d'être utilisée. 


Compare le travail effectué à celui qui était à faire et délivre un signal d'erreur 


Comparateur proportionnel à la différence entre une grandeur de référence (E) et la 


P randeur physique issue de la chaîne de retour. 
ou détecteur 9 physiq 

d'écart Ce signal d'erreur, après amplification, agira sur les organes de puissance 
dans un sens tel que l'erreur tendra à s’annuler. 


Le régulateur se compose d'un comparateur qui détermine l'écart entre la 
Régulateur consigne et la mesure et d'un correcteur qui élabore à partir du signal d'erreur l’ordre de 
commande. 


C'est l'organe d'action qui apporte l'énergie au système pour produire l'effet 


Actionneur ie 
souhaité. 


Le capteur prélève sur le système la grandeur réglée (information physique) et la 
Capteur transforme en un signal compréhensible par le régulateur. La précision et la rapidité sont 
deux caractéristiques importantes du capteur. 


On appelle perturbation tout phénomène physique intervenant sur le système qui 
Perturbation | modifie l’état de la sortie. Un système asservi doit pouvoir maintenir la sortie à son niveau 
indépendamment des perturbations 
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1- 2.4 - Régulation et systèmes asservis 


Nous avons fait la distinction dans l'introduction entre régulation et asservissement. Nous pouvons 
maintenant préciser de façon nette cette différence : 


i Un régulateur : maintient l'erreur € entre l'entrée E et la sortie S nulle, quelles que soient les 
perturbations, la grandeur d'entrée E restant constante ou variant par palier. E est alors appelée 
consigne ou référence. 


i Un système asservi : maintient l'erreur € nulle ou minimale quelles que soient les variations de 
E. Généralement, E est une fonction du temps qui peut être périodique, mais qui doit toujours 
rester continue et finie. 


Il faut remarquer que les contraintes sont plus grandes pour un système asservi que pour un 
régulateur, puisque aucune contrainte de vitesse de variation n’est imposée pour E. 


1- 2.5 - Propriétés des systèmes linéaires 


Quand un système est linéaire, il jouit de propriétés importantes qui permettent une étude plus 
commode, en particulier le « principe de superposition linéaire » qui se traduit par les relations : 


Entrée Sortie 

e(t) => s(t) 

Additivité : eb(t) 5 S(t) 
ei(t)+e,(t) => s(t) +S2(t) 


où e(t) et s(t) sont les 
grandeurs d'entrée et de sortie 


s(t) 


Homogénéité : 
a.s(t) 


i 
7 o 
D | 
ax 
Uy 


| Ce principe traduit le fait que les effets sont proportionnels aux causes et que les causes 
| ajoutent leurs effets. 


1- 3 - Régimes transitoires des asservissements 
1- 3.1 - Définitions 


Entrée Entrée dun système dont l'expression, en fonction du temps, est du type 
Permanente | constante, linéaire, parabolique ou périodique 


CAN Il est atteint par un système quand, soumis à une entrée permanente, sa sortie est du 
Régime même type que l'entrée c'est-à-dire constante, linéaire, parabolique ou périodique. 

Permanent . | OON i 

Ce régime est aussi appelé régime forcé. 


Régime Il correspond au fonctionnement du système quand il passe d'un type de régime 
Transitoire || permanent à un autre. 


Pratiquement, un asservissement travaille toujours en régime transitoire ; en effet, même un régulateur 
dont l'entrée est constante doit constamment revenir au régime permanent, car des perturbations qui 
constituent des entrées secondaires l’en écartent. Il en est de même pour les asservissements. 
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L’aptitude du servomécanisme à revenir au régime permanent sera caractérisée par ses performances 
dynamiques. 


1- 3.2 - Performances d'un système asservi 


* En régime permanent : la grandeur de sortie doit être aussi voisine que possible de la valeur 
désirée. En réalité, il subsiste toujours une légère erreur. Cette erreur est appelée : 


- erreur statique ou écart permanent quand la grandeur d'entrée est une constante ; pour 
un système idéal, elle doit être nulle. 


- erreur de trafnage quand la grandeur d'entrée est une fonction linéaire du temps. 


* En régime transitoire : le système évoluant entre deux régimes permanents, le temps mis par le 
système pour aller de l’un à l’autre et la façon dont il parvient à l’état final, sont très importants. 


- Le temps de réponse est le temps au bout duquel la sortie du système a atteint, à 
+ 5% (ou + 2 % selon la précision voulue), sa valeur de régime permanent et y reste 
(Fig. 1-5). 

-  L'amortissement : la sortie du système dépasse généralement la valeur qu'elle doit avoir 
dans le régime permanent final et elle oscille quelques instants autour de cette valeur. 
Les oscillations doivent être amorties, le plus rapidement possible. L’amortissement est 
mesuré par le coefficient À de l’exponentielle enveloppe (Fig. 1—6). 


sortie sortie 
+5% 
-5% 
t > 
temps de 
reponse Fig. 1-5 : Temps de réponse Fig. 1-6 : Amortissement 


1- 4 - Mise en équation d'un système - Résolution 
1- 4.1 - Mise en équation 


Nous avons dit précédemment que nous nous bornions à l'étude des systèmes linéaires. Donc, les 
équations rencontrées seront des équations différentielles linéaires à coefficients constants. 


Considérons un système quelconque A, le plus général possible, possédant une entrée e(t) et une 


sortie s(t) (fig. 1—7). 


Fig. 1-7 : Représentation d’un système quelconque à 1 entrée — 1 sortie 


Si on applique un signal à l'entrée, on recueillera, à la sortie, un signal qui sera liée au signal d'entrée 
par une équation différentielle de type : 
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* Les coefficients ai et bj sont les paramètres du système et ils sont sensés être connus, ce qui est 
le cas dans la pratique pour la plupart des systèmes courants. lls représentent diverses 
constantes de temps et divers coefficients de proportionnalité accessibles à la mesure. 


* La difficulté de la mise en équation réside surtout au niveau de la connaissance du processus lui- 
même. En réalité, l'équation différentielle à laquelle on arrive n'est souvent qu'une approximation 
qui consiste à négliger des termes d'ordre plus élevé. Cette précision suffit dans la plupart des 
cas, bien qu'une étude plus poussée soit quelque fois nécessaire. 


* Une fois l'équation du système établie, il faut exprimer la valeur de la sortie en fonction du temps 
pour connaître les régimes permanents et transitoires. Pour cela, il existe 2 méthodes 
(voir fig. 1—8) : 


Consiste à résoudre l'équation différentielle décrivant ce système, c’est- 
Méthode à-dire trouver une réponse forcée et une réponse libre pour le système. Mais 
Classique cette méthode ne permet pas toujours de trouver une solution et peut amener à 
une difficulté de résolution dès que l’ordre de l'équation différentielle dépasse 2. 


Basée sur le calcul opérationnel ou, essentiellement, sur la transformée 
de Laplace qui mettra en relation, une fonction de la variable du temps f(t) avec 


Méthode une fonction de la variable complexe F(p) dépendant de la pulsation. 
Opérationnelle 
Notation : £i f(t) ] = F(p) avec p=a + j.b (nombre complexe) 
elt) Méthode classique (ordre < 2) s(t) = £- [S(p)] 


Transformée de 
Laplace du signal 
d'entrée 


Transformée inverse de 
Laplace du signal de sortie 


Calcul 
opérationnel 


of [e(t] = E(p) S(p) = [s(t)] 


Fig. 1-8 : Détermination de la sortie du système par la méthode 
classique et par le calcul opérationnel 


Pour un rappel sur l'utilisation de la transformée de Laplace, voir l'annexe A. 


1- 4.2 - Utilisation de la transformée de Laplace 


En appelant S(p) et E(p) les transformées de s(t) et de e(t), si on prend la Transformée de Laplace des 


deux membres de l'équation différentielle : 


an Hadui sa tags = b ERa sus. +b, À +bo € 


On aura : 


k 
d'où : S(p) = . E(p) 
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Si l'on connaît l'image E(p) de e(t), il est facile, grâce aux tables de transformées de Laplace, de 
revenir à l'original de S(p). 


D'une manière générale, cette notation n’est valable que si : 


le système est linéaire à coefficients constants, 
toutes les variables et leurs dérivées sont nulles pour t < 0 (le système part du repos absolu), 


le système est dissipatif, donc sa réponse tend, plus ou moins, vers un régime permanent 
indépendant des conditions initiales. 
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Chapitre 2 : NOTION DE FONCTION DE TRANSFERT 


2- 1 - Introduction 


Rappelons que : 


* SI nous considérons un système quelconque A, le plus général possible, possédant une entrée 


e(t) et une sortie s(t) (fig. 2—1) : 


Fig. 2-1 : Représentation d'un système quelconque à 1 entrée — 1 sortie 


* ALORS, Si on applique un signal à l'entrée, on recueillera, à la sortie, un signal qui sera liée au 
signal d'entrée par une équation différentielle de type : 
d's(t) ds(t) 


a, —— "+... +a +4s = b 
n dt" 1 dt 0 k dt* 


En appelant S(p) et E(p) les transformées Laplace de s(t) et de e(t), si on prend la Transformée de 
Laplace des deux membres de l'équation différentielle, on aura : 


a, p° S(p) +...... +a p S(p)+ao S(p) = bę p“ E(p)+...... +b; p E(p) + bọ E(p) 


d'où : S(p) = 


C'est aussi le rapport de la transformée de Laplace de la sortie à la transformée de Laplace 
de l'entrée quand toutes les conditions initiales sont nulles. Dans ce cas, on a : 


S(p) = F(p) . E(p) 


La Fonction de Transfert caractérise la dynamique du système. Elle ne dépend que de ses 
caractéristiques physiques. Ainsi, dorénavant, un système sera décrit par sa fonction de transfert et non par 
l'équation différentielle qui le régit. 


Notons enfin, que cette fonction de transfert est aussi appelée transmittance par analogie avec 
limpédance dans les systèmes électriques. 
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2- 2 - Fonction de transfert d'un ensemble d'éléments 
2- 2.1 - Eléments en série (ou cascade 


Soit n éléments de fonction de transfert G1 (p) -...... Gn (p) mis en série (la sortie du premier est reliée à 
l'entrée du second, etc...) (fig. 2—2). 


T E E En LS E s 
— G ı (p) Š LENEE G n (p) D E o S 


Fig. 2—2 : Connexion en série (ou cascade) de fonctions de transfert 


La fonction de transfert de l'ensemble est égale au produit des fonctions de transfert de chaque 
élément : 


2- 2.2 - Eléments en parallèle 


Soient n éléments de fonction de transfert G1 (p) ....... Gn (p) mis en parallèle (fig. 2—3). 


l E S 
------ - e re 


Fig. 2-3 : Connexion en parallèle de fonctions de transfert 


La fonction de transfert équivalente H(p) a pour expression : 


On peut considérer que S(p) est le résultat de la superposition des n sorties des n éléments, c'est-à- 
dire que : 


S(p) = Sı (p) + S2 (p) +... + Sn (p) (en vertu de la linéarité du système, les effets s'ajoutent) 


Chaque élément pris, indépendamment, donnera une sortie Si (p) quand on lui applique l'entrée E(p). 
Donc : 
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S(p) = D Si(p) = G1 (p) . E(p) + G2(p) . E(P) + …… + Gn (p) . E(p) 


S(p) = [ G1 (p) + G2 (p) + ....... + Gn (p) ] . E(p) 


d'où : H(p) = G1 (p) + G2(p) +... + Gn (p) 


2- 2.3 - Cas d'un système à n entrées indépendantes 


E1 


Système D 
En 


Fig. 2-4 : Système à n entrées indépendantes 


La fonction de transfert n’a de sens qu'entre la sortie et une entrée. Le système de la fig. 2-4 pourra 
donc se décomposer en n constituants ayant la sortie en commun et pour entrée chacune des n entrées. 


On calculera les fonctions de transfert Gi (p) de chaque élément en supposant nulles les entrées autres 
que Ei (p). Ceci n'est possible que si les différentes équations du système ne sont pas couplées entre elles. 


Dans ce cas, on peut écrire : 


I ny a pas de fonction de transfert globale pour le système. 


2- 3 - Fonction de Transfert en Boucle Fermée FTBF 


Soit un système asservi, le plus général, représenté par le schéma de la fig. 2-5. 


Fig. 2-5 : Schéma fonctionnel d'un système asservi (Boucle Fermée) 


Soit A(p) et B(p), respectivement, les fonctions de transfert des chaînes directe et de retour. 


Cherchons la fonction de transfert du système complet : H(p) = S(p) 


Nous avons les relations suivantes : 
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S(p)=A(p).E(p) S’ (p) = B(p). S(p) , E (p) = E(p) - S° (p) 
doù S(p}= —*(P)__ Elp) 


La fonction de transfert d'un système bouclé ou en Boucle Fermée (FTBF) est donc le rapport 
de la fonction de transfert de sa chaîne directe à 1 + A(p). B(p.: 


A(p) 
H(p) = 
| P= TEA) BP) | 
2- 4 - Fonction de transfert en boucle ouverte FTBO 


La Fonction de Transfert en Boucle Ouverte (également appelée F.T.B.O.) est la fonction de transfert 


qui lie les transformées de Laplace de la sortie de la chaîne de retour S’ (p) à l'erreur € (p). Elle correspond à 
l'ouverture de la boucle (Fig. 2-6 ): 


Fig. 2-6 : Schéma fonctionnel d'un système asservi en Boucle Ouverte 


Dans ce cas, € = E puisque le comparateur ne reçoit plus qu'une seule information. 


On a donc: S'’(p) = B(p). S(p) 


= B(p) . A(p) . E(p) 
= B(p) . A(p) . E(p) 


d'où : ——— = K(p) = A(p) . B(p) 


La Fonction de Transfert en Boucle Ouverte (ou FTBO) d'un asservissement est le produit 
des fonctions de transfert de la chaîne directe par la chaîne de retour. 


| 
| La fonction de transfert en boucle ouverte a une grande importance dans l'étude de la 
| stabilité des systèmes ; de plus, elle est directement accessible à la mesure. 
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2- 5 - Fonction de transfert d'un système à boucles multiples 


régulateurs ou servomécanismes dans une chaîne. La figure 2-8 en donne un exemple. 


Fig. 2-8 : Exemple de système asservi à boucles multiples 


son inverse 1 / Y(p). 


O Afp) a 1 
YP) = TTAB) Wa aE 


Dans notre cas : 
B(p) =B6(p) transmittance de la chaîne de retour 


A(p) transmittance de la chaîne directe A(p) = Aı (p). Yı (p) . Y2 (p) 


En appliquant la même procédure, on a : 


ep. f =B5 + l ; Lepp 
ETE ie ht it 
AT AYY, A À AzA; A, 


| 1 1 1 1 1 
Soit: —— =Be6+ B» + Bs + B4 + 
Y(p) al í z] É xl i | 


Soit un système asservi représenté par sa fonction de transfert de forme générale suivante : 


F(p) Bnp” Grise +Bo _ N(p) 
App’ Ses +A5  D(p) 
Si N(p) et D(p) ont des racines alors : 
N(p) = Bm [(p — Z1)(p — z2)....…..(p — Zm)] ou N(p) = Bm | [(p-z:) 
i=1 
D(p) = An [(p — p1)(p — p2)------- (p — pn)] ou D(p) = An ITe-P) 


2- 6 - Formes générales de la Fonction de Transfert d'un système linéaire 


Il existe des systèmes complexes où l’on rencontre, non seulement une chaîne de retour principale, 
mais un grand nombre de chaînes de retour secondaires. Dans ces asservissements, il y a plusieurs 


Le calcul de la fonction de transfert d'un tel système peut paraître compliqué. Pour mener à bien ce 
calcul, il faut utiliser l’artifice suivant : au lieu de considérer la fonction de transfert globale Y(p), on considère 
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F(p) s'écrit alors : F(p) = —© i- Avec: m<n pour un système réel 


* 


les racines du numérateur sont appelées " zéros de la fonction de transfert", 


* 


les racines du dénominateur sont appelées " pôles de la fonction de transfert ", 


2- 7 - Règles de transformation des schémas fonctionnels 


D'une manière générale, pour simplifier un bloc fonctionnel il est souvent plus judicieux de déplacer les 
points de connexion et les comparateurs (ou additionneurs), d’inter-changer ces derniers, puis de réduire les 


boucles internes. 


Schéma fonctionnel original Schéma fonctionnel équivalent 


A-B A-B+C 


A A.G: A.G1.G2 A A.G2 A.G1.G2 
3. G: G2 G2 |__| G 


4 
A A.G1.G 
GiG | 


A A.G1+A.G 
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10. 


11. 


A.(G1+G2) 
—> 


12. 


13. 


14. 
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Règle n° 14 : 


Règle n°13 : 


Règle n° 13 : 


2- 7.1.a - Exemple de réduction successive d'un schéma fonctionnel 


Soit à réduire le schéma fonctionnel suivant : 


GG;.G; 
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Chapitre 3 : METHODES D'E TUDES DES ASSERVISSEMENTS 


3- 1 - Introduction 


Nous avons vu, dans le chapitre précédent, qu'il était possible connaissant les équations différentielles, 
de déterminer la fonction de transfert d'un système. Mais il existe de nombreux cas où le système est un 
système industriel mal défini et dont, à fortiori on ne connaît pas les équations différentielles. Or, la 
connaissance de sa fonction de transfert est très importante pour déterminer ses performances et surtout sa 
stabilité. II est donc important de mettre au point des méthodes capables de résoudre le problème. 


En général, on applique cette procédure pour déterminer les fonctions de transfert des éléments qui 
entrent dans une chaîne. La connaissance expérimentale ou mathématique de toutes les fonctions de transfert 
des éléments permet alors de déterminer la fonction de transfert de l’ensemble. 


Ces méthodes sont basées sur l’utilisation d'entrées dites canoniques, faciles à mettre en œuvre dans 
toutes les techniques (électrique, mécanique, hydraulique). On en déduit alors les différentes constantes de la 
fonction de transfert. 


Certains appareils dit analyseurs de fonction de transfert facilitent les mesures. 


3- 2 - Entrées canoniques 
3- 2.1 - Echelon unité 
C'est une fonction nulle pour t < 0 et constante et égale à 1 pour 0 < t < œ (fig. 3—1). 


Cette fonction est appelée quelquefois u(t) (unité). Elle n’est pas définie pour t = O puisqu'il y a 
discontinuité à cet endroit. 


u(t) 


Fig. 3—1 : Fonction Echelon 


Sa transformée de Laplace est : 


rom 
p 


3- 2.2 - Echelon de vitesse ( rampe unité ) 


C'est une fonction nulle pour t < 0 et qui varie linéairement avec t pour t > 0 (fig. 3—2). 
On l’exprime parfois sous la forme r(t) = t . u(t). 


Cette fonction est appelée échelon de vitesse ou rampe, car sa vitesse de variation est constante et 
égale à 1. 
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Fig. 3—2 : Fonction Rampe 


On vérifie aisément que sa transformée de Laplace est égale à : 


r 
oË {r(t)}= — 
p 
En effet : 
L{r(t)}= fetat, on pose : u=t dv = e™ dt 
0 
-pt 
du = dt V= cE 
—P 
er | [eri + 
Donc {r(t)}=|f(t}=—| =0+ —|- = — 
=P Jo p P lo P 
3- 2.3 - Echelon d'accélération 
0 pour t<0 


Soit f(t) la fonction échelon d'accélération, définie par (voir fig. 3—3) : f(t) = < t2 
—u(t) pourt>0 
2 


f 


0 


Fig. 3-3 : Fonction Accélération 


2 2 2 
Æ { f(t) } = fert., on pose : ust dv = e™ dt 
3 2 2 
du = t.dt de 
—P 


Donc {f(t) } = [CES =0+ 1 Pir) = A 
-=P 


D'où : 
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3- 2.4 - Impulsion unitaire 


Une impulsion est une fonction du temps de durée très courte mais dont l'amplitude est suffisamment 


grande pour que l'effet en soit sensible. L'impulsion est dite unitaire si la surface est égale à 1. On la note ä(t) 
(fig. 3—4). 


0 pour t<0 et tr 1/t 
S(t) = 


lim 1 pour O<t <r 
t0 T 


Fig. 3—4 : Fonction Impulsion ä(t) 
Toutes les impulsions, dont la durée égale numériquement l'inverse de l'amplitude, sont unitaires si 
cette durée tend vers zéro. 


Pour t = 0, l'amplitude est théoriquement infinie. 


+00 
ô(t) est définie par : f Slt). dt= 1 (ce qui est équivalent à la surface unitaire) 


—00 


Elle est appelée aussi impulsion de DIRAC. 


Calculons sa transformée de Laplace. 


Pour cela, définissons la fonction f(t) (fig. 3-5) telle que : 


1 pour O<t <rt 0 z 
Fig. 3—5 : Fonction f(t) 


wh pour t<0 et t>t 


f(t) peut être considérée comme la différence entre deux échelons unitaires dont l'un est décalé de +: 


f(t = u(t) — u(t 7) > L {it }=L {u(t)} -L {ult- 1) }= Ia - e`?) 
p 
Or: ô(t) = lim 1 
t—0 T 
Donc ; LL) = lim -e™?) 
t0 tp 
= lim all-e) 
t0 (xp) 
D'où : 
| nE | 
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3- 2.5 - Entrée harmonique (D 
Elle est définie par (Voir fig. 3-6): 


+ 


0 pour t<0 
f(t) = > 
A sin(ot+@) pour t20 


Fig. 3—6 : Signal harmonique 


Sa transformée de Laplace est : 


_ p.sinp+o.cos o | 
E 2 


| Č { f(t) } 


p? +o 


2 et £ {cos ot} = —P (Voir annexe B) 
pta? pta? 


car č {sin ot} = 


3- 3 - Réponse d'un système asservi aux entrées canoniques 
3- 3.1 - Réponse du système à une impulsion unitaire : réponse impulsionnelle 


Soit un système de fonction de transfert H(p). 
Appliquons sur son entrée une fonction, e(t) = &(t), c'est-à-dire une impulsion unitaire. 


Sa sortie sera donnée par : S(p) = H(p) . E(p) 
Or: E(p)=1, puisque o {8(t)}=1. 


Donc la transformée de Laplace S(p) de la sortie correspond exactement à la fonction de 
transfert H(p). 


S(p) = H(p) 


C'est aussi une autre définition de la fonction de transfert. On voit donc qu'une méthode pour 
connaître H(p) est de mesurer la réponse à une impulsion unité. 


La Fig. 3-7 montre deux types de réponses impulsionnelles (selon la nature du système à exciter). 


s(t) 
s(t) 
t À ue 


Fig. 3-7 : Exemples de réponses impulsionnelles 


Du point de vue pratique, cette méthode présente quelques difficultés, car il est pratiquement 
impossible de réaliser physiquement une entrée ô(t). On se contente, en général, d'une impulsion de durée 
aussi courte que possible mais finie, d'où une certaine imprécision. Après avoir envoyé cette entrée &(t) 
approchée, on doit enregistrer, en fonction du temps, la réponse s(t). Ce qui donne une courbe qu'il faut 


Cours d’'asservissements linéaires continus (2013-2014) Licence et Master « Automatisme Industriel » (Prof. FELLAH M.K.) 25 


ensuite interpréter. Si on veut l'expression mathématique de la fonction de transfert, on approche cette courbe 
par des morceaux de courbes correspondant à des fonctions connues. Il faut alors prendre la transformée de 
ces fonctions du temps pour obtenir la fonction de transfert. 


Cette opération, facile à décrire, est incontestablement délicate à réaliser. Elle est bien entendu 
entachée d'erreurs, mais aucune autre méthode n'est parfaite. 


3- 3.2 - Réponse du système à un échelon unité : réponse indicielle 


Pour pallier aux inconvénients de la réponse impulsionnelle, il est plus facile, pratiquement, d'utiliser 

comme entrée, un échelon unité. 
L'entrée du système est donc : e(t) = u(t), d'où E(p)= 1 
p 


Sa sortie est alors : 


S(p) = Hp) C'est l'intégrale de la fonction de transfert. 


p 


Pratiquement, l'essai est très rapide, il suffit de faire passer l'entrée de O à une valeur constante 
pendant un temps (donné) T puis de cette valeur à 0, et d'enregistrer la sortie en fonction du temps. 


En principe, cette sortie doit être du même type que l'entrée si le système est linéaire, c'est-à-dire qu’au 
bout d’un certain temps correspondant à la durée du régime transitoire, la sortie doit rester constante. 


S'il en est autrement, le système n'est pas linéaire (la réciproque n'est pas forcément vraie). Donc ce 
test permet de savoir si on est en présence d'un système linéaire ou non. 


La Fig. 3-8 montre deux types de réponses indicielles (selon la nature du système à exciter). 


“4 . s(t) 
en 


Fig. 3-8 : Exemples de réponses indicielles 


L'exploitation des résultats pour déterminer la fonction de transfert est beaucoup plus laborieuse : Il 
faut décomposer l'échelon d'entrée en série de Fourier ainsi que la réponse de sortie en supposant que le 
phénomène a une période T. 


Pour chaque fréquence élémentaire, on connaîtra les amplitudes de l'entrée et de la sortie, ainsi que 
leur déphasage. 


On en déduit, à chaque fréquence, le ” gain "et la ” phase ” du système. A ce stade, on se retrouve 
dans le cas de la réponse fréquentielle que nous allons étudier. 


Cette méthode nécessite peu d'essais, mais beaucoup de calculs. La méthode suivante (réponse 
fréquentielle) nécessite beaucoup d'essais et peu de calculs. 
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3- 3.3 - Réponse fréquentielle 


Plutôt que de faire une décomposition mathématique des signaux d'entrée et de sortie en fonctions 
sinusoïdales élémentaires, il est apparu plus simple d'appliquer à l'entrée du système une entrée harmonique 
et de relever la sortie correspondante, une fois le régime permanent établi. 


On répète cette opération en partant des fréquences les plus basses jusqu’à des fréquences 
suffisamment élevées pour que la sortie du système soit négligeable. En effet, à partir d’une certaine fréquence 
le système ne 
" suit " plus les oscillations que l’on voudrait lui imposer, exactement comme le tympan qui ne vibre plus au- 
delà d'une certaine fréquence. 


De ces mesures, on déduit, en fonction de la fréquence : 
-So 
0 
(So et Eo amplitudes maximales des signaux sinusoïdaux d'entrée et de sortie à la fréquence o). 
e la phase ọ (a). 


e les amplitudes d'entrée et de sortie, donc le "gain" du système, qui est le rapport (G(w)| 


On représente ensuite ces fonctions de la fréquence par des courbes qui sont, soit : 
e en coordonnées logarithmiques (plan de Bode) 

e en coordonnées polaires (plan de Nyquist) 

e la courbe de phase (plan de Black) 


| De ces courbes, on peut déduire la fonction de transfert et bien d'autres paramètres. 


3- 4 - Représentation de la réponse fréquentielle 
3- 4.1 - Courbes de Bode et diagrammes asymptotiques 


La fonction de transfert pouvant être assimilée à un gain complexe, il suffit de tracer les courbes de 
variation du module (ou gain) et du déphasage en fonction de la fréquence (ou de la pulsation). 

En effet, tout nombre complexe G(o) = [Gloeo 

a pour logarithme : log {G(@)} = log {|G(@)|.e® } = log [Glo] +j elo) 


Le logarithme d’un nombre complexe se représentera donc par deux valeurs : le logarithme du module 
qui sera la partie réelle et l'argument la partie imaginaire. C’est sur ce principe qu'est basée la courbe de Bode. 


Par exemple, si la fonction de transfert est : G(p)=1 + Tp 


Alors, en régime sinusoïdal, p = jo, et on aura : 
G(o) = (G(o)Le =1+joT avec: [G(w)| = 41+(oT)? module ou gain 


et: ZG(o) = arctg (oT) phase ou déphasage 


Pour une représentation de Bode, on tracera les courbes suivantes : 


e Courbe de gain (Fig. 3-9) : Log |G(œ)| = log 41+(@T)? en fonction de (log œ) 
e Courbe de gain (Fig. 3—10): ZG(o) = arctg (œT) en fonction de (log œ) 
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[GI =1 
log IG] =0 Q= 
log © = œ% 
o>0 

log © >- œ% 
= > 
es N 
1 sA log œ 
ROS | o SION 
ee -> Affaiblissement :222°\ 
122 OCR OO GORE OPDOROORDOEERECN 

e | pe 
log |G| = — œ | 
Fig. 3—9 : exemple de courbe de gain pour une représentation de Bode 
À /G(a) 
> 
log œ 
o>0 
ZG > 0° o — 0 
ZG =>- 90° 


Fig. 3—10 : Exemple de courbe de phase pour une représentation de Bode 


3- 4.1.a - Diagramme asymptotique 


Dans la plupart des cas (nous préciserons par la suite), il est inutile de tracer complètement la courbe 
de Bode (point par point). On peut se contenter des asymptotes. 


| On appelle diagramme asymptotique d’une fonction de transfert l'ensemble des asymptotes à 
la courbe quand œ — 0 et quand œ — œ. Les points d'intersection de ces asymptotes entre elles sont 
appelés points de cassure. 


Etant donné l’utilisation des logarithmes comme échelle, la précision est suffisamment bonne pour que 
l'on puisse assimiler la courbe à ses asymptotes sauf au voisinage des points de cassure. 


Si la courbe de gain peut être avantageusement remplacée par ses asymptotes, il n’en est pas de 
même pour la courbe de phase qui s'en écarte beaucoup plus. 


3- 4.1.b - Rappels : 


e On appelle octave, l'intervalle qui sépare une fréquence fı d’une fréquence f2 = 2 . fı. Franchir une 
octave c’est doubler la fréquence initiale. 


e On appelle décade, l'intervalle qui sépare une fréquence fı d’une fréquence f2. = 10 . fı 


Cours d’asservissements linéaires continus (2013-2014) Licence et Master « Automatisme Industriel » (Prot. FELLAH M.K.) 28 


3- 4.2 - Courbe de Nyquist ou lieu de Nyquist 


Le problème est toujours le même, il s’agit de représenter un nombre complexe, variable avec la 
fréquence, c’est-à-dire trouver le lieu du point d’affixe [Goje quand o varie de 0 à œ. Cette courbe 
définit entièrement les propriétés du système. 


On la trace généralement en coordonnées polaires. Comme pour les courbes précédentes, il faut 
mesurer le gain et le déphasage du système. 


On porte alors, pour chaque valeur de œ, un rayon vecteur dont la longueur est égale au module du 


gain : G(w)| SE et qui fait un angle ọ (œ) avec laxe des réels (Fig. 3-11). Quand o varie, on obtient une 
S 


0 
courbe graduée en o. 


Fig. 3-11 : Courbe de Nyquist 


On complète généralement cette courbe par symétrie par rapport à laxe réel pour avoir sa 
représentation quand o varie de 0 à — o. 


Cette partie de la courbe n’a, physiquement, aucun sens, mais il est nécessaire de la tracer si on veut 
étudier mathématiquement la stabilité de l’asservissement. 


On peut ainsi compléter la courbe parce que les coefficients de l'équation différentielle qui régit le 
phénomène sont supposés constants et réels. 


Comme les courbes de Bode, le lieu de Nyquist peut se tracer soit expérimentalement, soit à partir de 
la fonction de transfert si elle est connue, c'est pourquoi on l'appelle aussi " lieu de transfert dans le plan de 
Nyquist ". 


3- 4.3 - Courbe amplitude - phase ou lieu de Black (ou Black - Nichols 


Le lieu de Black concilie les avantages des courbes de Bode et ceux du lieu de Nyquist : la fonction de 
transfert est représentée par une seule courbe et les échelles sont logarithmiques (Fig. 3—12). 


log|G| 


o= 10 


l l l 


—180°% -90° 0°) —90° 


œ= 104 


Fig. 3—12 : Courbe de Black-Nichols 
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Ce dernier avantage est très important du point de vue pratique car, généralement, le module du gain 
varie dans des proportions énormes ce qui oblige pour le lieu de Nyquist à prendre des échelles assez 
grandes, supprimant toute précision pour les gains faibles. Cet inconvénient disparaît avec des échelles 
logarithmiques. 


Le lieu de Black se trace de la façon suivante : 


On porte en ordonnée le module du gain { log G(w)|}, et en abscisse le déphasage { ọ (œ) }. On 
obtient une courbe paramétrée par œ. 


Cette représentation est la plus employée dans l'industrie. L'intérêt pratique de cette représentation est 
évident : connaissant les fonctions de transfert de plusieurs éléments en cascade, on aura la fonction de 
transfert de l'ensemble par une simple addition vectorielle : 


Soit: G=G1.G2.G3 


|Gļ.e}? = G|.e1°1 | Gale}? ; IG 3/93 


d'où: log |G| = log |G 1| + log |G 2| + log |G 3] 
et P =Pi+P2+9p3 
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3- 5 - Etude des systèmes du premier ordre 
3- 5.1 - Définition 
On appelle système du 1% ordre, un système régi par une équation linéaire différentielle du premier 


ordre telle que : 


TEU | {t) = Ke(t) 
dt 


ou encore, un système dont la fonction de transfert est du type : S(P) ss 


E(p) 1+Tp 


Ces systèmes sont encore appelés systèmes à une seule constante de temps, ou système à 
retard. 


lils sont très nombreux en physique. En dehors des circuits RC ou RL en électricité, on peut considérer 
qu'un amplificateur est un système du 1°’ ordre. En mécanique, tous les assemblages comportant un ressort et 
un amortisseur sont du 1° ordre. Nous les passerons en revue à la fin de ce paragraphe. 


3- 5.2 - Réponse indicielle 
La réponse indicielle nous renseignera sur le comportement du système en régime transitoire. 


K Elp) ici E1(p)=1/p sipe 


p) = —— 
1+Tp p(1+ Tp) 


En consultant une table de Transformée de Laplace, on voit que l'originale s:(t) de S:(p) est : 
sı (t) = K(1-e-t/T) Lim sı (t) = K 


t—> œ 


On constate donc que la sortie s1 (t) (Fig. 3—13) atteint pratiquement le régime permanent au bout d’un 
temps qui dépend de la constante T. 


Cette constante T, appelée constante de temps, caractérise donc la rapidité du système à atteindre 
son régime permanent. 


e(t) 


entrée 


Fig. 3—13 : Réponse indicielle d'un système du 1°% Ordre 


La pente de la tangente à l’origine est K/T, plus le système a une constante de temps faible, plus il 
"répond" vite. 


Au bout d’un temps t = T, la sortie est sı (T) = K(1 — 1/e), ce qui représente environ 63% de K. 
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Temps de réponse : Nous avons vu que le temps de réponse était le temps au bout duquel la sortie 
avait atteint son régime permanent à 5% près. Dans le cas du système du premier ordre, ce temps correspond 
à 3T environ. 


3- 5.3 - Réponse à une rampe ( échelon de vitesse) 


Dans ce cas, nous avons : e2(t) = t.u(t) où u(t) : échelon unitaire 
E2(p) = 1/p° 
Donc S2(p) = 3 a = S1(p) avec S:(p) : transformée de la réponse indicielle d’un 1° Ordre. 
p°(1+Tp) p 


t t 
d'où : sat = | s1(t) dt = | K(1—e-tUT}dt 
0 0 

sat) =K{t-T(1-e-tT)} 


La figure 3—14 donne la réponse à une rampe d’un système du 1°" Ordre. 


Fig. 3-14 : Réponse d'un système du 1° Ordre à une rampe (tracée pour K=1) 


En régime permanent : t > œ, alors st) =K(t-T) 


On met, ainsi, en évidence le retard T qui constitue une erreur permanente. 


Dans le cas où K=1, l'erreur T entre l'entrée et la sortie est constante en fonction du temps. C'est 
l'erreur de "traînage”. 


Donc, un système du 1°" ordre suit les variations linéaires de l'entrée avec un certain retard, d'où leur 
nom de système à retard. 


3- 5.4 - Réponse à une impulsion unité 


Dans ce cas, nous avons : est) = ô(t) où 6 (t) : impulsion unitaire 
1 
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On met ainsi en évidence la constante de temps sur le graphique (Fig. 3-15). 


Fig. 3—15 : Réponse impulsionnnelle d'un système du 1°" Ordre 


3- 5.5 - Réponse fréquentielle 


3- 5.5.a - Courbe de Bode 


Courbe de Gain : 


En régime permanent sinusoïdal, la fonction de transfert a pour expression : 


K 2 
. Glo)| = ==> e(t) = E, cos(ot) 
a a Gloeo avec eo) Lor et | | 
Eo 1+joT Sotool s(t) = So cos(ot + p) 


(Fig. 3—16) 
Log |[G(w 
|[G(o) | EC.) Log K DEN Asymptotes 
\ 
(a) 
0 
GI =1 Log œo =Log (1/T) 
Log 1=0 
x 


Traçons la courbe de gain |G(œ)| en coordonnées cartésiennes et en coordonnées logarithmiques 


pente ©) 


Log œ 


Fig. 3—16 : Courbe de gain en coordonnées cartésiennes et en coordonnées 
logarithmiques d'un système du 1°" Ordre 


Log (K/T) —/ 
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On voit que la courbe est fortement dilatée pour les valeurs faibles de œ, et condensée pour les valeurs 
élevées. On remarque également que la courbe diffère très peu de ses asymptotes qui sont ici : 


e Pour o>0, log | Glo) | = log K — log Ÿ1+ T°? > log K 
e Pour O>, log | G(o) | = log K — log œT = log (K/T) — log œ > -logo 


Puisque l’abscisse est donnée par log œ, quand œ> œ , alors log | G(œ) | > — log œ, c'est à dire vers 
une asymptote de pente (—1), qui : 


e coupe l'horizontale ( log K) pour œ = 1/T 
e et coupe l'axe des (log œ) pour œ= K/T. 


Quand o> 0 , alors log | G(o) | > log K = Constante, c’est à dire vers une droite de pente (0). 


Les 2 asymptotes se coupent en : Log K = log KT —log © > œ= 1/7 (point de cassure) 
L’axe des log œ (pour un gain unitaire) est coupé en: Log K/T -log œ = log 1 > o = K/T 
Courbe de Phase : 


Il faut représenter ọ = Arctg (-o«T) en fonction de log œ (Fig. 3—17) : 


° o>0, p=0 
e >, yp=-90° 
° œ@=1/T, ọ=-—45° 


0° : Log 1/T Log. 


_45° Point d'inflexion 
a PE correspondant au 
chanaement d’asvmotote 
—90° 


Fig. 3-17 : Courbe de phase d’un système du 1°" Ordre 


La courbe (Fig. 3-17) est loin d’être assimilable aux asymptotes surtout pour œ = 1/T. 
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3- 5.5.b - Courbe de Nyquist 


Il est facile de trouver les différents points de cette courbe (Fig. 3—18) : 


e o>0, [|G(0)I=K, 


. >, 


Imag 


Fig. 3-18 : Courbe de Nyquist d’un système du 1°" Ordre 


D'autre part, si on sépare les parties réelles et imaginaires de la fonction de transfert, on a 


So K K KoT 


-< - 
Eo 1+joT {49% 149% 
X Y 


Or: (X-K2)2+V2=(K/2 } 


équation d’un cercle de rayon K/2 et de centre (K/2 , 0). 


On a représenté un demi cercle en trait plein pour indiquer qu'il correspond au fonctionnement réel du 
système quand la pulsation œ croît de 0 à +, et en pointillés pour les pulsations "négatives" de —~ à 0. 


3- 5.5.c - Courbe de Black 


On trace cette fois log |G(œ) | en fonction de ọ(œ). La courbe a l'allure de la figure 3—19. On peut, 
également, la déduire des courbes d'amplitude et de phase du diagramme de Bode. 
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45° 


Arctg(—K) 
œ=K/T 


o> o 


Fig. 3—19 : Courbe de Black d’un système du 1°" Ordre 


3- 5.6 - Exemples de systèmes du 1°" ordre 
3- 5.6.a - Filtre "passe-bas" 


C'est un montage qui laisse passer les fréquences basses et atténue fortement les fréquences élevées. 


3- 5.6.b - Filtre "passe-haut" 


LEE T=RC 
Ve 1+RCp 
Le 1 T=L/R 
Ve 1+—p 
R 
1 
r-J to Q f 
dt P ado 
f 


f : coefficient de frottement visqueux 


T=J/f 


Ces dispositifs atténuent les fréquences basses et laissent passer intégralement les fréquences 
élevées. lls jouent donc le rôle inverse des filtres passe-bas. 
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vs. BCP T=RC 


Va 1+RCp 
R m 
p 
YR T=L/R 
Ve Vs f Ve 1 + R° 
L 
; f Tp 
Ces fonctions de transfert sont du type : G(p) = 
1+Tp 
oT 
|G(w)| = ———— 
| 272 
Dans ce cas : Lu T , et le diagramme asymptotique est celui de la figure 3—20. 
tg ọ (0) = —— 
oT 


Co) 


Fig. 3-20 : Courbe de Bode d'un système du 1°" Ordre de type 
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° Pour o>0: 


log | G(w) | = log œT -log V1+T?6? > logo +logT > log (asymptote de pente +1) 
tg p(@) > +œ, p(o) = 90° 


e Pour @ > oo: 
Log | G(o) | = log œT — log oT > 0 (asymptote de pente 0) 
tg p(o) > 0, p(o) = 0° 


On remarque que ces systèmes atténuent les fréquences basses et laissent passer intégralement les 
hautes fréquences. 


Le diagramme de Nyquist (Fig. 3-21) est le même que pour les filtres passe-bas, mais il est gradué en 
sens inverse. 


Imag 


© 


Tp 
1+Tp 


Fig. 3-21 : Courbe de Nyquist d'un système du 1°" Ordre de type 
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3- 6 - Etude des systèmes du second ordre 
3- 6.1 - Définition 


Les systèmes du second ordre sont régis par des équations linéaires différentielles à coefficients 
constants du 2è"e ordre, du type : 


Aoe(t) 
2 t 
6° _ ce CU LE pla fetat 
dt 0 
À de(t) 
dt 
A 
Leurs fonctions de transfert seront du type : S(p) _ 3 Ao 41/p 
E(p) B:p°+Bp+B 
p 


Le comportement du système sera extrêmement différent suivant que le degré qui figure au 
dénominateur aura des racines réelles ou imaginaires. 


On introduit les paramètres suivants : 


s(t) 


e Gain statique : K= — C'est le rapport —— en régime statique co ay ; — =0) 
Bo e(t) dt dt? 


e  Pulsation propre non amortie : Œn = FR rad/s 
2 


e Facteur d'amortissement (sans dimension) : Es & = 1 (valeur critique) 


e Constante de temps : T= —— 


La fonction de transfert s'écrit en fonction des paramètres ainsi définis : 


K 


14 2 p+ 7 pe 
On on 


Prenons le cas de H(p) = 


et étudions les différentes réponses. 


Le comportement dynamique d'un tel système dépend de la valeur des deux constantes œn et 
surtout de §&. 
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e  Sié>1: Le polynôme est décomposable, le dénominateur a 2 racines réelles (—p1 et -p2) : 
-pı =- on (E+ 467-1) 


-p2 =- on (6 - yE? -1) 


e Si 0O<E<1: Le polynôme n’est pas décomposable, le dénominateur a 2 racines complexes 


conjuguées (-po et -pġ ) (Fig. 3-22) : 


Spy lEt) 


—Po =- œn (—j41-¢° ) 


Imag 
—Po _ +2 
J@n 1 E 
FA a N a 
A --77 i On 2 
Pôles l @n y1-6 
N B 
EN Du Eon 0 Réel 
N e—a 
© 
— Po 


Fig. 3—22 : Relations entre les paramètres d’un système du 2" ordre dans le plan Complexe 


SOn 


On 


Remarque :  E=cosf = : coefficient d'amortissement 


e Sič=1: Les 2 racines sont égales ( —p12): 


— P1,2 = — @n > P1.2 = on = 1/t 


3- 6.2 - Réponse à un échelon unité 


K 1 
E(p)=1/p S(p) = ŒÆ 1 
P 1+ p+ pt 
On oh 
e  Si0<E<1: Les 2 racines imaginaires conduisent à une solution oscillatoire amortie. 
Le régime permanent est ici s(t) = K (lim s(t) = lim S(p) = K) 


t> 00 p>0 
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La solution est : 


D -xf _ +250, a 
p pP?+2Eonp+ On 


P p+ëo on 
p_(p+Ea)/+@f(1-E7) (p+Eu,)? +a(1-E?) 


Or z e) e “nt cos(œn 41-6? )t 
et zo| Onvi=E 
( 


- ea ét ei 7 {= 2\t 
E A 


En posant : 
+ sin a=4/1-E? (<1) et Cos a= E (<1) 


1- 2 1- 2 
6 > a = arctg ( 6 


> tga= 


-ont i= 2 
On aura : s(t) = K| 1- S vie 


.sin{o, 41- E°t + arctg 


1-2 
6 1444442444 448 
® 


e  Siëé=1: Les 2 racines (pôles) sont égales : système amorti critique. 


K A K A 
S(p) = 2- = in 
p pf+20P+onr PPan) 
On aura s(t) = Kh -e+ ot) (t>0) 
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e Sič >1: Les 2 racines (pôles) sont négatives et inégales : système apériodique. 


2 
K oʻ 


AR p (p+Eu, sonde? -1[p+ 0, one -1) 


1 SGIRE] Jont 1 L- 2 Jont 
st) =K 11—— e "TE (t>0) 
2f? if + fer) 2e -16-467 -1) 
-a4 -a2 z nl 6 yE? -1 
s(t) =K| 1+ Pn $ n -| (t20) avec Ae ) 
2JE2-1l a az az =on( č- 467-1) 


La Figure 3—23 donne les réponses indicielles (2D et 3D) d’un système du second ordre en fonction du 


coefficient d'amortissement &. 


Amplitude 
7 


Amplitude 


Temps 


Amortissement 


Fig. 3-23 : Réponses indicielles d'un système du 2"4 ordre en fonction de & 
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3- 6.3 - Réponse à une impulsion unité 


K 2 
E({p) = 1 S) = — 
p° + 2Ë@p + Of 
e Si0<E <1: s(t) = K ee sin ont 1-8) t>0 
1=E 
e  Sié =1: s(t) = Klotet] t20 


e SiE >1: s(t) = i „o t>0 


La Figure 3-24 donne les réponses impulsionnelles (2D et 3D) d'un système du second ordre en 
fonction du coefficient d'amortissement &. 


Amplitude 


Amplitude 


Amortissement 


Fig. 3-24 : Réponses impulsionnelles d'un système du 2" ordre en fonction de & 
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3- 6.4 - Réponse fréquentielle 


3- 6.4.a - Diagramme de Bode 
On se place en régime permanent sinusoïdal : 


e(t) = Eo.sin(t) 


s(t) = So.sin(ot+p) puisque le système est considéré linéaire. 


p=)o > 
Solo K K K 
|G()| = LOIS = - = 
o(® o. O . o l o 2 2 2 
1+2—j+|— j 1-| — | p+j426— 1-2 26 © 
On On (OM (OM o, o, 
25 2 
p(o) =/G(œ) = -arctg T 
1- 2) 
On 
Module : 
2 2 
log |\G(@)| = log K- log 1-11 | +26- 
On On 
% pour <<, log [G(w)| > logK (asymptote horizontale), 
2 
+ pourw>> a, log |[G(w)| > log K — log e = (log K + 2 log œn ) — 2 log © 
On 
(asymptote de pente —2). 
Les 2 asymptotes se coupent en : log K = log K — 2 log e > pour @= @n 
On 
Phase : 
26 © 
p(o) =- arctg h 7 
1-1 
On 
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# pour o—0, (0) =0° 
v À ©=GM, (@n) = — 90° 


% pour o—>co , ps) = — 180° 


La figure 3-25 donne le diagramme de Bode en fonction du coefficient d'amortissement &. 


Bode Diagram 


Magnitude 


Pente (—2) 


Phase (deg) 


= On 
Frequency (rad/sec) 


Fig. 3-25 : Diagramme de Bode d’un système du 2™ ordre en fonction de & 
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Phénomène de résonance : 


Pour œ = mn, le module peut tendre vers des valeurs très grandes quand le système n’est pas très 
amorti : 


Gae 
2% 


Sié > 0, alors  |G(@n)| > œ 


On voit donc que dans ce cas particulier, la courbe est très loin de l’asymptote et que |[G(o)| a une 
valeur bien supérieure à celle prévue par le diagramme asymptotique. Il faut donc tenir compte d’un 
phénomène nouveau appelé " résonance " qui se produit dans certaines conditions. 


Cette résonance a lieu à la pulsation œr pour laquelle |G(œr)| est maximum. On peut calculer œr en 
annulant la dérivée du dénominateur de [G(o)|. On trouve : 


OR = @n 1-2? ( pour 0 < é < 0.707 ) @r : pulsation de résonance. 


Dans ce cas, le gain maximum vaut : 


On définit le coefficient de résonance Q par le rapport du gain maximum au gain pour les fréquences 
très basses ( œ >0 ) : 


1 


Q = Gumu / K = —— 
2641- 6° 


3- 6.4.b - Courbe de Nyquist 


On trace cette fois |G(w)| el°{°) en coordonnées polaires. 


L’allure de la courbe (Fig. 3-26) peut être déduite des courbes de Bode. On voit notamment que : 


% Pour o=0, | G(o) | =K et (o) = 0° 
+% Pour, [G(æ)|>0 et ọ(@) = — 180° 


Il est difficile de mettre en équation la courbe, car elle n’a pas une forme classique connue. 


e Pour un système très peu ou pas amorti, la résonance a lieu pour œ = 1/T , donc pour tg ọ = — œ, 
soit ọ = — 90° 

e Plus 6 est grand, plus la courbe est " petite " et se rapproche de celle d'un système du 1°% ordre 
(demi cercle). 


e  Quandoœo >, p > -— 180°. Les courbes sont donc tangentes à l'axe réel. 


Cours d’asservissements linéaires continus (2013-2014) Licence et Master « Automatisme Industriel » (Prof. FELLAH M.K.) 46 


O —> oo @ — 0 


3 
> 
g 
F 
£ 
0 
Real Axis 
Fig. 3-26 : Lieu de Nyquist d'un système du 2" ordre en fonction de & 
Remarque : 


Nous avons mis en évidence différentes fréquences ou pulsations ; il convient de ne pas les confondre: 


e Pulsation propre non amortie @n = 1/T (pulsation de cassure) 
e Pulsation propre amortie Op = En 41- g? 
e Pulsation de résonance OR = On 41-26? 


Dans le cas d’un système très peu amorti ( & très petit), les 3 pulsations peuvent être confondues. 


3- 6.4.c - Courbe de Black 


Les courbes de Black permettent de mettre facilement en évidence la résonance car elle correspond à 
un maximum de la courbe ; Plus le système est amorti, plus le maximum est faible. 


L’allure du lieu de Black est donnée sur de la figure 3-27. 
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Gain 


-180 -135 


3- 6.4.d - Exemples de systèmes du 2°"° ordre 


Ve HE 
C 


Vs (p) 1 Vs (p) 


-90 
Phase (deg) 


RCp 


Va(p) 1+RCp =+LCp° 


ci: me et Le 
Te 2VL 


Va(p) 1+RCp+LCp? 


@ — 0 


Fig. 3-27 : Lieu de Black d'un système du 2" ordre en fonction de & 


ÉGp” 


V, (p 


— 


Vap) 1+RCp+LCp? 
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Chapitre 4 : STABILITE DES SYSTEMES ASSERVIS LINEAIRES 


4- 1 - Notion de stabilité d'un système 
4- 1.1 - Définition de la stabilité 


On dira qu'un système linéaire est stable si, après avoir soumis son entrée à une brusque 
variation (échelon unité, par exemple) : 


| ° le mouvement amorcé par sa sortie reste borné en amplitude (c'est à dire que la sortie 
| garde une valeur finie) 


° ce mouvement s'amortit plus ou moins vite et la sortie tend vers un état d'équilibre. 


Les réponses indicielles des figures 4-1 et 4-2 correspondent à celles de systèmes stables. Nous 
retrouvons les critères cités ci-dessus. 


s(b s(t) 


Fig 4-1: Système oscillatoire amorti (stable) Fig 4-2 : Système non oscillatoire amorti (stable) 


La figure 4-3 est un cas de système instable. Les oscillations sont de plus en plus importantes et le 
système ne retrouve pas son état d'équilibre. 


Physiquement, un système instable dont la réponse croit sans limite peut se causer des dommages ou 
en causer à autrui (danger pour l'être humain). En pratique, la majorité des systèmes sont conçus avec des 
dispositifs de limitation. 


Si on considère le cas où des oscillations persistent indéfiniment (cas du pompage de la figure 4-4), on 
peut considérer le système comme stable (système marginalement stable) puisque sa sortie garde une valeur 
finie, à condition que l'amplitude ne soit pas trop grande. 

s(t) 
s(t) 


Fig 4—3: Système oscillatoire divergent (instable) Fig 4—4: Système oscillatoire 
(marginalement stable) 
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La stabilité est une condition impérative. Pour que les systèmes soient utilisables en asservissement, il 
est absolument nécessaire que les fonctions de transfert en boucle fermée FTBF soient stables. Ceci 
n'implique toutefois pas que les FTBO soient stables. 


4- 1.2 - Aspect mathématique de la stabilité 


Considérons un système asservi quelconque dont la fonction de transfert est : 


Hp) = —2P 
1 + A(p)-B(p) 


Si on envoie sur l'entrée un échelon unité E(p)=1/p, alors : 
1 A 
S(p) (p) 


- p1+A(p)B(p) 


Nous avons vu au chapitre 3, que S(p) pouvait se mettre sous la forme d’un quotient de polynômes du 


type pa et que celui-ci pouvait se décomposer en une somme de fractions rationnelles : 
S(p) = Np) = i i + Gi 
D(P) p-p; P -Pn 


où les pi sont les racines réelles ou complexes de D(p). 


Prenons, par exemple, le cas où le dénominateur contient des racines nulles (pôles multiples), des 
racines réelles (pôles réels) et des racines complexes. C'est-à-dire qu'il est de la forme : 


D(p) = pP™%(P-P1)-----(P-Pn, {p -01)? +02]...{p-an,)? +02 ] 


La décomposition de S(p) en fractions rationnelles sera : 


i=no A. k=n; B, j=n Cpp+D, 


p k=1 P — Pk T p-a)? +o; 


Les racines complexes étant aj +j œj (œj partie réelle, œj partie imaginaire), cherchons l'original s(t) de 
S(p) qui est la réponse du système à un échelon unité. On trouve : 


i=no A, 
s(t) = A1+ D 
i=2 


tD 
ü- 


k=n; j=n i 
+ $B- + D Fe” .Sin(ot+6;) 
ki F1 


On constate donc que la sortie garde une valeur finie quand t > œ, si les conditions suivantes sont 
remplies : 


e Les px et les oj doivent être négatifs pour que les exponentielles correspondantes soient 
décroissantes. 


° Les Ai doivent être nuls sauf A1. 


Nous verrons dans la suite que pour certaines fonctions de transfert, la présence de pôles multiples 
nuls n’entraîne pas forcément une augmentation infinie de la sortie. 


En effet, les termes en 1/p ont une action d'intégration, leur influence peut être combattue par des 
actions de dérivation provenant de terme en p au numérateur. S'il n’en est pas ainsi, le système possédant des 
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4- 1.3 - Conditions de stabilité 


| transfert n'a de partie réelle positive. 


Cela exclut : 
e Les racines réelles positives. 


On peut formuler ceci autrement : 


gauche du plan complexe. 


oscillatoire à variation constante lorsque t > œ. 


Imaginaire 


(2 racines complexes à 
partie réelles < 0) 


pôles multiples à l’origine au dénominateur de sa fonction de transfert est dit " intrinsèquement instable ", c'est- 
à-dire qu'aucune modification des coefficients ne peut le rendre stable. 


Un système linéaire est stable si aucune des racines du dénominateur de sa fonction de 


e Les racines complexes à parties réelles positives. 


Un système asservi bouclé est stable si tous les pôles de la FTBF sont localisés dans le demi-plan 


Un système asservi bouclé est instable si sa FTBF comprend, au moins, un pôle localisés dans le 
demi-plan droit du plan complexe et/ou des pôles de multiplicité > 1 sur l'axe imaginaire. 


Si le système comprend une seule paire de pôle sur l'axe imaginaire ou un pôle unique à l'origine, 
le système est dit marginalement stable. Sa réponse sera oscillatoire non amortie ou non 


La figure 4-5 récapitule les cas possibles suivant le signe et la nature des racines. 


Régime oscillatoire divergent 
(2 racines complexes à partie 
4 réelles > 0) 


X Í 
Réel 


Régime non oscillatoire amorti 
(1 racine réelle < 0) 


Régime oscillatoire non amorti 
(2 racines imaginaires pures) 


Régime non oscillatoire divergent 
(1 racine réelle nulle) 


Fig. 4-5 : Récapitulatif des comportements des systèmes selon la position et le signe des pôles 
(selon les réponses indicielles). 


Régime non oscillatoire divergent 
(1 racine réelle > 0) 
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Mais les conditions de stabilité ainsi définies ne sont pas suffisantes pour caractériser un système 
asservi : un système très mal amorti sera inutilisable, il faudra donc toujours définir des marges dites de 
sécurité sur les coefficients d'amortissement. 


Remarque à propos des systèmes instables : 


Quand on a affaire à un système instable, sa sortie tend théoriquement vers l'infini si on soumet son 
entrée à une brusque variation. En réalité, sa sortie ne tend pas vers l'infini, mais vers une valeur qui 
correspond à la saturation. Cette valeur peut être très grande et conduire à la destruction du système. En tout 
état de cause, dans le cas où la fonction de transfert a des pôles à parties réelles positives, le système sort 
rapidement de son domaine de linéarité et ses équations ne sont plus valables. 


4- 1.4 - Etude de la stabilité d'un système bouclé 


Le système asservi bouclé de la figure 4-6 a pour fonction de Transfert : 


Fig. 4-6 : Schéma fonctionnel d'un système asservi bouclé 


Sa stabilité est conditionnée par le signe des parties réelles des racines du dénominateur. 


Il suffira, donc, d'étudier l'équation : 1 + A(p).B(p) = 0, et de chercher le signe de ses racines. 


Plusieurs moyens sont possibles pour y arriver : 


a) 1° moyen: Calculer les racines de 1 + A(p).B(p) = 0 


Cette méthode est bonne puisqu'elle nous donne également les valeurs des racines en 
plus de leurs signes. Mais elle est pratiquement inapplicable à cause de la grande difficulté 
qu'elle présente si le degré du polynôme est important. L'usage d’un ordinateur peut simplifier le 
travail, car il peut aussi tracer le lieu des racines quand on fait varier les paramètres. C'est une 
méthode très puissante. 


b) 2°" moyen: Discuter le signe des racines sans les calculer, à partir des coefficients du 
dénominateur (critère de Routh-Hurwitz ) 


Malheureusement, si le système trouvé est instable, on ne sait pas sur quel paramètre il 
faut agir pour le rendre stable. Il faut en plus connaître la fonction sous sa forme mathématique. 


c) 3" moyen: Utiliser le critère de Nyquist (méthode graphique). 


Cette méthode est intéressante car elle n’a pas les inconvénients du critère de Routh. A 
savoir, on peut utiliser directement les résultats expérimentaux sans connaître les équations du 
système et elle montre graphiquement sur quels paramètres on peut agir pour rendre le système 
stable. 
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4- 2 - Critère de Routh - Hurwitz 
4- 2.1 - Enoncé du critère 


Soit P(p) le dénominateur de la fonction de transfert en boucle fermée. P(p) peut être écrit 
sous la forme : 


P(p) = 1 + A(p).B(p) = ao p" + ar pr +... + &n-10 + än 
(Équation caractéristique de la fonction de transfert en boucle fermée) 


Pour que le système soit stable, il faut et il suffit que les racines de P(p) n'aient pas de 
parties réelles positives. 


4- 2.1.a - Critère d’Hurwitz 


| Ce critère (nécessaire mais pas suffisant) indique que le système est instable si les ai sont de 
| signes différents ou certains sont nuls. 


4- 2.1.b - Critère de Routh-Hurwitz 


| La condition nécessaire et suffisante de stabilité est alors que tous les termes de la 1° 
| colonne du tableau de Routh soient de même signe. 


On construit le tableau de Routh de la manière suivante : 


Pour P(p) = ao p? + a pt +... + An-1 p + an, le tableau de Routh est simplement une matrice 
carrée avec une ligne pour chaque puissance de p dans le polynôme de l'équation caractéristique. 


p’ :| e {1ère ligne: coefficients des termes en pr-2k (avec k = 0,1, 2,...) 
DEA : | e 2ème ligne : coefficients des termes en p™=2k+1) (avec k = 0,1, 2,...) 
p" :| e 3ème ligne : combinaison des 2 lignes précédentes 

p° i| e 4ème ligne : combinaison des 2 lignes précédentes 

N N 

N N 

p° ‘| e dernière ligne: combinaison des 2 lignes précédentes 
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Si, par exemple, P(p) = ao p” + a: pê + a2 p° + as pt + a4 p? + as p? + aep + a7, alors le tableau de Routh 
se construit comme suite : 


p: ao a2 a4 as | 
pê: ai a3 as a7 | 
pi b _ 442 -8083 b _ 2434 — dods b _ &êç — Aoa7 
| a, a, a: 
pi À _ bras — abs ca = 28s -abs à _ bsa7 -34.0 ” 
bı bi bı 
cb; —-b.c cb; —b.c 
p: di -SEa ORs das 1C5 
C; C] 
d:c; — c;d d4Cs — c:.0 
p: e - 23 =M ega MM5 TOTS Sis 
d; d; 
pi f _ €193 — des 
e] 
fe; — €4.0 
p°: AA 
fi 


| Routh a démontré que le nombre de " pôles instables " (c'est-à-dire le nombre de pôles à 
partie réelle positive) de la fonction de transfert en boucle fermée est égal au nombre de changement 
de signe que comporte la 1*® colonne, lue de haut en bas. 


Si ce nombre est différent de zéro, alors le système est instable. 
Remarques : 


Cette méthode a l'avantage d'être rapide est exacte, mais elle ne donne pas une mesure de 
la stabilité comme les autres critères ; car elle se borne à dire si le système est stable ou non. De 
plus, elle est inapplicable si on ne connaît pas l'expression mathématique de la fonction de transfert. 


Le critère de Routh est intéressant pour connaître le nombre de racines réelles positives, 
mais il est incapable de donner des renseignements sur l'amortissement du système quand celui-ci 
est stable. 


La méthode est cependant en défaut dans les 2 cas suivants : 
e Si tous les coefficients d'une ligne sont nuls. 


e Si un terme de la 1ère colonne de gauche est nul à l'exclusion des autres termes de la même 
ligne. 
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4- 2.2 - Exemple 1 


P(p) = p° + 3p? + 3p +1 +K (K gain variable) 
1ère condition (critère d'Hurwitz) : 1+K>0 > K>-1 
2ème condition (critère de Routh-Hurwitz) | tableau de Routh 
p8 : 1 3 
p? : 3 1+K 
p': 9-(1+K) 
3 
p° : 1+K 
9-(1+K) an 
Pour que le système soit stable, il faudrait que : 3 
1+K>0 
C'est-à-dire: -1<K<8 (condition nécessaire et suffisante de stabilité) 
Si cette condition n'est pas vérifiée, c'est-à-dire, si : 
e K«<-1, il y a 1 seul changement de signe dans la 1ère colonne; donc un seul 
pôle instable. 
e K>8, il y a 2 changement de signe dans la 1ère colonne; donc 2 pôles 
instables. 
e Si(K=-1 ou K= 8), (frontière entre la stabilité et l'instabilité) on dit que le système est 


oscillant (marginalement stable). 


4- 2.3 - Exemple 2 (ligne complète de zéros) 


Nous avons dit que si une ligne complète était composée de zéro, la méthode était en défaut. En fait, il 
est quand même possible d'en tirer des conclusions moyennant certains aménagements. 


Si P(p) = p° + 7p* + 6p? + 42p? + 8p + 56 


Alors, le tableau de Routh est : 


pë : 1 6 8 | 

p4 : 71 426 568 Division de la ligne par 7 
p° 0 0 

p? : z 2 

p! : z 

p° : 3 
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La 3ème ligne est nulle. On substitue alors à cette ligne les coefficients obtenus en différentiant une 
fonction fictive, appelée polynôme auxiliaire, construite sur la ligne précédant la ligne nulle. Le polynôme 


auxiliaire pour l'exemple en cours s'écrit : 


Q(p) = pf + 6p? + 8 


Si nous le dérivons, par rapport à p, nous obtenons alors : 


dQ(p) 


~~ = 4p? + 12p +0 


dp 


Les coefficients de ce polynôme remplacement ceux de la ligne nulle dans le tableau initial. Le tableau 


devient alors : 


1 6 

1 6 
4>1 12>3 

3 8 
1/3 

8 


Il ny a aucun changement de signe sur la 1? colonne du tableau, donc aucune racine à partie réelle 
positive. Le système est donc stable. 


4- 2.4 - Exemple 3 (un zéro sur la première colonne) 


Si le premier élément de la ligne est nul, la ligne suivante ne pourra pas être calculée car il y aurait une 
division par zéro. Pour éviter cela, on utilise un nombre de valeur très faible € (epsilon) pour remplacer le zéro 


de la première colonne. € peut tendre vers zéro par valeur positive ou négative, pour permettre par la suite le 
calcul du nombre de changement de signe de la première colonne. 


Considérons le système dont la FTBF G(p) = 


10 


p? +20" + 3p° + 6p°? + 5p +3 


P(p) = pē + 2p* + 3p° + 6p? + 5p +3 


Alors, le tableau de Routh est : 


p? : 


p! : 


p° : 


1 3 
2 6 
LEZ 7/2 
Ge — 7 3 
E 
42s — 49 - 6s? 
12514 
3 
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Considérons uniquement le changement de signe dans la première colonne et calculons le signe de 
chaque ligne dans les 2 cas (e > 0+ ets > 0): 


1ère colonne € > 0+ e >20- 
p: 1 + + 
p : 2 + + 
p°? 0e + - 
p? Ge —7 — + 
E€ 
p' : 42s -49 6e? + + 
12-14 
p’ : 3 + + 


Si e est choisi positif, il y a 2 changements de signe. S'il est choisi négatif, il y a également 2 
changements de signe. Le système a donc 2 pôles dans le demi-plan droit du plan complexe (2 pôles 
instables) et ce n’est pas important si nous choisissons d'approcher le zéro par valeur positive ou négative. 
Ceci est toujours le cas. 


4- 3 - Critère de Nyquist 


Le critère de Nyquist permet de déterminer la stabilité d'un système bouclé sur la base de sa réponse 
harmonique en boucle ouverte. 


4- 3.1 - Énoncé du critère de Nyquist 


| La condition nécessaire et suffisante de stabilité d'un système asservi linéaire est que son 

| lieu de transfert en boucle ouverte, parcouru de œ = — à œ = +œ, entoure le point critique (—1,0) dans 

| le sens trigonométrique un nombre de fois égal au nombre de pôles instables de la fonction de 
transfert en boucle ouverte. 


Étant donné un système asservi, défini par sa fonction de transfert en boucle ouverte 
FTBO(p) = A(p).B(p). 


La relation : Z=P-N donne le nombre Z de zéros instables de l'équation 
caractéristique1 + FTBO(p) = 0 et donc de pôles instables de la FTBF(p), avec : 


e P: Nombre de pôles instables de la FTBO(p), 


e N: Nombre de tours que fait le lieu complet de Nyquist ( œ variant de —~ à + ) autour du point 
critique (—1,0) dans le sens trigonométrique ( sens anti-horaire ). 


En particulier, le système asservi est stable, à condition que : Z=0 > P=N 
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En pratique, on retiendra les étapes suivantes pour appliquer le critère de Nyquist : 


e Étudier la stabilité de la FTBO > P: nombre de pôles instables de la FTBO. 


e Tracer le lieu de Nyquist complet de la FTBO ( œ variant de -œ à +œ). 


e Calculer le nombre de tours (comptés algébriquement dans le sens trigonométrique), soit N, que 
fait le lieu complet de Nyquist ( œ variant de -œ à + ), autour du point critique (-1,0). 


e En déduire Z = P — N = nombre de pôles instables de la FTBF. 


4- 3.2 - Exemple 


Soit un système asservi à retour unitaire dont la FTBO est : G(p)= —— (T>0) 


Discutons sa stabilité suivant les valeurs de K. 


e K>0 (fig. 4-7) 
v La FTBO(p) a un pôle instable p = +1/T. 


> P=1 
y Le nombre de tours autour du point (-1,0) est : 
> N=0 


Z=P-N=1#0 
1 pôle instable de la FTBF 


> Système instable en boucle fermée. 


v Ce système est instable en boucle ouverte et 


instable en boucle fermée. Fig. 4-7 


e K<-1 (fig. 4-8) 
{v P=1 
{y N=+1 
Ÿ Z=P-N=0 


Imag. 


O> —00 


v Pas de pôle instable de la FTBF 
> Système stable en boucle fermée. 


v Ce système est instable en boucle ouverte et 
stable en boucle fermée. 


> 
> +0 Réel 


Fig. 4-8 
e —1<K <0 (fig. 4—9) 
⁄ P=] /Mimag 
Ÿ N=0 Pet 
Ÿ Z=P-N=1 o=0 K 0 To> -0 
v ôle i i > 
1 pôle instable de la FTBF 1 o=0 > +0 Réel 
> Système instable en boucle fermée. 
v Ce système est instable en boucle ouverte et 
instable en boucle fermée. ; 
Fig. 4—9 
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Vérification la stabilité par le critère de Routh : 
K K 


FTBO(p) = —— > FTBF(p) = ———— (retour unitaire) 
1-Tp 1+K-Tp 
Équation caractéristique : —Tp+1+K=0 
p! : =T 
p? | 1+K 
Pour que le système soit stable, il faudrait que : 1+K< 0 > K<-i 


4- 3.3 - Critère de Nyquist simplifié (critère du Revers) 


Il est déduit du critère de Nyquist : 


Nyquist de la FTBO, décrit dans le sens des pulsations croissantes (œ variant de O à +), laisse le 


Un système stable en Boucle Ouverte, est stable en Boucle Fermée, si le tracé du lieu de 
| point critique (-1,0) à sa gauche. Il est instable dans le cas contraire. 


Le critère du Revers est d'un emploi plus commode que le critère de Nyquist, car il ne met en œuvre 
que le lieu physique de la FTBO (correspondant aux pulsations positives) 


Par contre, le critère du Revers est moins général, et il ne peut s'appliquer sans danger que lorsque la 
notion de "gauche" n’est pas ambiguë (Fig. 4—10). 


Imag. 


> 
-1 Réel 


Fig. 4—10 : Cas d'’ambiguïté dans 
l'application du critère du revers 


4- 4 - Marges de stabilité 


Pour un système stable en Boucle Ouverte, nous venons de voir que la stabilité en Boucle Fermée 
dépend de la position du lieu de Nyquist de la FTBO par rapport au point critique (-1,0). 


Le critère de Nyquist spécifie que le lieu de Nyquist doit laisser le point —1 à gauche lorsqu'on le 
parcourt dans le sens croissant des œ. Le cas où il existerait une pulsation à laquelle le lieu traverserait 
exactement ce point est un cas limite correspondant à un système en boucle fermée dont la stabilité serait 
marginale. 


Mais la tendance vers l'instabilité est graduelle : plus le lieu de Nyquist est proche du point critique, 
moins le degré de stabilité est bon, et plus on aura par exemple d'oscillations avant stabilisation en boucle 
fermée. 


De façon à quantifier le degré de stabilité d’un système asservi, il est donc utile de chiffrer la distance 
entre le lieu de Nyquist et le point critique (—1, 0). La mesure effective de la distance minimum n'étant pas 
chose aisée d’un point de vue mathématique, on préfère, de manière traditionnelle, évaluer indirectement cette 
distance par les mesures des marges de phases Av et de gain AG. Ces marges représentent des marges de 
sécurité par rapport à l’état instable. 
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Ces grandeurs sont définies de la manière suivante : 


e La marge de phase Ag d'un système est mathématiquement la différence entre la phase de 
FTBO(wcO) c'est-à-dire ọ (wc0), et —-180° : AG = 6 (wc0) + 180° 


La marge de phase permet de préserver la stabilité en dépit de la présence de retards parasites 
— par exemple dus à la transmission des signaux — dont on n’a pas tenu compte au moment de 
l'étude de la stabilité. 


1 


e La marge de gain AG a pour expression : ac x: 
FTBO(%,)| 


Elle permet de préserver la stabilité en dépit des fluctuations de gain, qui affectent, en particulier, 
les amplificateurs de la chaîne de puissance. 


| Pour qu'un système soit stable, il faudrait que : A®9>0 et AG>0 


Ces marges sont illustrées sur le lieu de Nyquist des figures 4-11 et 4-12. 


Fig. 4-11 : Illustration des marges de gain et Fig. 4-12 : Illustration des marges de 
de phase sur le lieu de Nyquist gain et de phase sur le lieu de Nyquist 
(cas d’un système stable) (cas d'un système instable) 


4- 4.1 - Valeurs usuelles de Ag et AG 

Les marges définies ci-dessus permettent d'évaluer la distance entre le point critique et le lieu de 
Nyquist en boucle ouverte. Imposer leurs valeurs revient à s'assurer que l’on ait jamais FTBO(o) = —1, c'est-à- 
FTBO(w)| =1 


dire jamais simultanément (pour la même pulsation) : 
ZFTBO(o) = -180° 


L'expérience montre que pour des systèmes classiques (notamment à phase minimale), un bon degré 
A = 45A 60° 


de stabilité en boucle fermée est obtenu si l’on est capable d'imposer : 
AG > 8A 15 dB 


Avec ces valeurs, on obtient dans la plupart des cas une paire de pôles dominants en boucle fermée 
caractérisés par un coefficient d'amortissement & = 0,5 ...0,707. 


Pour régler la stabilité d'un système, il est souvent délicat de raisonner en tenant compte des deux 
marges à la fois. Dans ce cas, on privilégie, en général, la marge de phase Av. 


Cours d’asservissements linéaires continus (2013-2014) Licence et Master « Automatisme Industriel » (Prof. FELLAH M.K.) 60 


4- 4.2 - Critère de Stabilité utilisant les courbes de Bode et de Black 


Les marges de gain et de phase définie précédemment peuvent également être représentées sur le 
diagramme de Bode (figures 4-13 et 4-14) et sur le lieu de Black-Nichols (fig. 4—15 et 4-16). 


AIG(o) | NG(o) | 


Fig. 4-13 : illustration des marges de gain et Fig. 4-14 : illustration des marges de gain et 
de phase sur le diagramme de Bode de phase sur le diagramme de Bode 
(cas d’un système stable : AG > 0 et Ap > 0) (cas d’un système instable : AG < 0 et Ap < 0) 

—180° 


-180° 


FTBO() 

Fig. 4-15 : illustration des marges de gain et Fig. 4-16 : illustration des marges de gain et de 
de phase sur le lieu de Black-Nichols phase sur le lieu de Black-Nichols 

(cas d'un système stable : AG > 0 et Ap > 0) (cas d’un système instable : AG < 0 et Aọ < 0) 


4- 4.3 - Marge de stabilité appliquée à la position des pôles de la FTBF 


La notion de marge de stabilité appliquée à la FTBF conduit à s’interdire un domaine pour la position 
des pôles dans le plan complexe. 


On s'impose, en général, B + 60 ° (ce qui correspond à un système du deuxième ordre 
avec & = cos p = 0.5) (fig. 4—17). 


Fig. 4-17 : illustration de la 
marge de stabilité imposée sur 
la position des pôles dans le 
plan complexe. 


Zone interdite 
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Chapitre 5 : PERFORMANCES DES ASSERVISSEMENTS 


5- 1 - Introduction 
Nous supposerons dans l'étude qui suit que les systèmes asservis étudiés sont stables. 


La précision d’un système est définie à partir de l'erreur € entre la grandeur de consigne E et la 
grandeur de sortie S (fig. 5-1). Nous distinguerons la précision statique qui caractérise la limite de l’erreur au 
bout d’un temps infini pour une entrée donnée, c’est-à-dire le régime permanent, et la précision dynamique qui 
tient compte des caractéristiques d'évolution du processus en régime transitoire. 


5- 2 - Performances Statiques des Systèmes bouclés 
5- 2.1 - Erreur statique 


La précision d'un asservissement, en régime permanent, est définie par l'écart permanent 
E(t) qui existe entre la sortie réelle et celle que l’on désire obtenir. 


Par définition, on dira qu'un système est d'autant plus précis que le signal d'erreur E(t) est 
plus faible. 


L'idéal serait que l'on ait : E() =0, Yt 


En pratique, il en est autrement, car : 
e La consigne peut varier : la recherche de la minimisation de e(t), en dépit de ces variations, 
constitue un problème de suivi (ou de poursuite). 


e Un signal de perturbation aléatoire (exemple : un bruit) peut venir de superposer au signal utile en 


un point de la chaîne : le maintien de e(t), en dépit de la perturbation, constitue un problème de 
régulation. 


Fig. 5-1 : Schéma général d'un asservissement 


Calculons l'erreur statique : 


> 2(p) = E(p) - e(p) - A(p) - B(p) > e(p}= — EP) (5-1) 


D’après le théorème de la valeur finale {lim f(t) = lim p.F(p)}, l'erreur statique €s (ou encore £+) est 
— 0 p> 


donnée par la relation : Es = lim e(t) = lim p.e(p) 
t= p—0 
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Es = lim p.—EP) = lim g= Erreur statique 
1+ FTBO(p) 


€s dépend donc à la fois : 


* du système considéré {présence de A(p).B(p)}, 


*% du signal d'entrée appliqué {présence de E(p)}. 


Adoptons la notation suivante : 


e n: le nombre dď’intégrateurs que comporte la FTBO(p) 


pr 1+ bp + bap? Pau 


e m:ordre du signal d’entrée canonique 


E(p) = E? feniera 
m 
p 
Sim=1 > échelon 
Sim=2 > rampe 
Sim=3 > Accélération 


5- 2.1.a - Erreur statique pour une entrée échelon 


C'est l'erreur qui subsiste en régime permanent sur la réponse indicielle (fig. 5-2). 


Entrée de 
UN référence 


e(t) = Eo.u(t) 


Erreur €s 
constante | 


+ 


Fig. 5-2 : Erreur statique pour une entrée échelon 


2 
FTBO(p) = A(p). B(p) = K 1+a,p+a,p° +... 0 


(5-2) 


n est appelé classe du système 
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Eo 


Si l'entrée vaut : E(p) = ra 


eai Ep) | _ Eo 
s= lim|p. = - 

p—0l 1+FTBO(p) 1+ lim FTBO(p) 

p—0 
Es = Eo Avec Ke = lim FTBO(p) = Constante d'erreur statique d'échelon 
1+K, p>0 
ou gain statique en Boucle ouverte 
P E 
e Pour les systèmes de classe 0 : Ke = lim K + ap +t sa ut =K > ss = — -cte 
p>0 1+b4p + b2p° +... 1+K 


2 
e Pour les systèmes de classe n > 0 : Ke = lim a ITAP EAP taris = 00 > es=0 


p>0 p? 1+b4p + bp? +... 


5- 2.1.b - Erreur statique ( ou erreur de traînage ) pour une entrée rampe (ou vitesse) 


C'est l'erreur qui subsiste en régime permanent sur la réponse à une rampe (fig. 5—3). 


; Erreur €s ļ 
/ Entrée de Les 


référence Ja 
e(t) = Eo.t.u(t) 


> 
0 
Fig. 5-3 : Erreur statique pour une entrée rampe 
Si l'entrée vaut : E(p) = 
p 

£s = lim — Eo o Eo 

p>0 P + TE E p + p.FTBO(p) lim p.FTBO(p) 

p> 

Es = Bi Avec Kv = lim p.FTBO(p) = Constante d'erreur statique de vitesse 

K, p> 

2 

e Pour les systèmes de classe 0: Kv = lim p.K. pepe = 0 D £s = 00 


p—0 1+ b4p + bap? ae 


K 1+8ap + Lu : 


e Pour les systèmes de classe 1:  Kv= lim p. —. z 
p>0 p 1+b;p+b:p° +... 


2 
e Pour les systèmes de classe n>1: K= lim p. E LEANE aN aa > 
p—0 p” 


1+ b4p + bap? + -nn 
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5- 2.1.c - Erreur statique pour une entrée parabolique (ou accélération) 


C'est l'erreur qui subsiste en régime permanent sur la réponse à une entrée accélération (fig. 5—4). 


Entrée de référence 


2 


e(t) = Ed. .U(t) 


+ 


Fig. 5—4 : Erreur statique pour une entrée accélération 


Si l'entrée vaut : E(p) = 


F 
£s = lim me Ey) lim Eo Eo 
S = = = 
p> 1+ HO p>0| p? + p? FTBO(p) lim p°FTBO(p) 
p= 
E | ; P 
Es = A = cte Avec Ka = lim p?.FTBO(p) = Constante d'erreur statique d'accélération 
a P= 
2 
e Pour les systèmes de classe 0: Ka = lim pê.K. APPAR iii =0 d Es = © 
p>0 1+ bp + Dap? + 
e Pour les systèmes de classe 1 : Ka = lim p°. aa À + ap + AAP? + onns =0 > Es =0 
p>0 p 1+ bp + bap? + 
2 
e Pour les systèmes de classe 2: Ka = lim Bb. = ITAP Ap tee =K > ss= Eo _ cte 
p? 1+b4p bop? a K 


2 
e Pour les systèmes de classe n> 2: Ka= lim p°. a : ia =0 gs=0 


p>0 p? 14+b4p+bop? +4 


5- 2.1.d - Récapitulatif des erreurs statiques 


Le tableau 5-1 récapitule les valeurs de l'erreur statique en fonction : 


e dela classe n du système 
e de l’ordre du signal d'entrée canonique 


e dugain K de la FTBO du système 
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CI d f | ; 
on de Constantes d'erreur statique Erreur statique £s 
Entrée Echelon Entrée Vitesse Entrée accélération 
n Ke K, K: a E 0 Eo Eg 
1+K, K, Ka 
0 K 0 0 Eo co co 
1+K 
1 Co K 0 0 Bo co 
K 
2 oo œ K 0 0 Eo 
K 
N 
co co 00 0 0 0 
N 


Tableau 5-1 : Récapitulatif des erreurs statiques 


Remarques importantes: 


Les constantes d'erreurs Ke, Kv, et Ka décrivent l'aptitude du système asservi à réduire ou éliminer 
l'erreur statique. Elles renseignent, par conséquent, sur les performances du système en régime permanent. 


Il est généralement préférable d'accroître les constantes d'erreurs, tout en maintenant la réponse 
transitoire dans des proportions acceptables ; En effet, l'erreur statique, lorsqu'elle est finie et non nulle, décroît 
lorsque le gain en boucle ouverte croît. Mais cette croissance du gain peut détériorer la stabilité du système. 
Cette propriété est connue sous le nom de " Dilemme Stabilité — Précision ", qui nécessite souvent un 
compromis. 


Il est à noter également que pour améliorer les performances en régime statique, nous pouvons 
augmenter la classe du système en ajoutant un ou des intégrateur(s) dans la chaîne directe du système. Ceci 
peut, cependant, engendrer des problèmes de stabilité supplémentaires. 


5- 2.2 - Gain statique en boucle fermée 


Pour un système stable, le gain statique en boucle fermée est défini par : 


Ke = lim SE) 2 im FTBF(p) 
( p—0 


t> e 


5- 2.3 - Exemple 


Soit le système asservi de la figure 5-5. Calculons ses différentes erreurs statiques pour différentes 
entrées canoniques : 
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Fig. 5-5 : Exemple d'un système asservi 


Il faut s'assurer, tout d’abord, de la stabilité du système. Utilisons pour cela le critère de Routh : 


2 
p 
1+5 
4 2 Pur 1 
FTBO(p) = = > FTBF(p) = D 
P+2)  pit+P) 1+ Pipoy 
2 
p? 1/4 1 
p! 1/2 0 Tous les termes de la 1è° colonne sont de même signe 
p° f 1 > système stable 
e Gain statique en boucle fermée : Ks = lim FTBF(p) =1 
p> 
e Entrée échelon: £s = lim pD =0 avec E(p) = 1/p (Eo = 1) 
p>0\ 1+ FTBO(p) 
ST $ Eo 
Ou alors : Ke = œ, puisqu'il s’agit d’un système de classe 1. Es = nK 1/00 = 0 
$ e 
e Entrée vitesse : €s = lim TERON = 1/2 avec E(p) = 1/p? (Eo= 1) 
p>0\ 1+FTBO(p) 
Ou alors : Kv = K = 2, puisqu'il s’agit d’un système de classe 1. €s = ro = 1/2 
v 
e Entrée accélération : €s = lim pee = 00 avec E(p) = 1/p° (Eo= 1) 
p>0\ 1+FTBO(p) 
Se TE na x Eo 
Ou alors : Ka = 0, puisqu'il s’agit d'un système de classe 1. Es = AE co 
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5- 3 - Performances Dynamiques des Systèmes bouclés 


Dans la majorité des cas pratiques, les caractéristiques de performances désirées pour les 
asservissements sont exprimées relativement au temps. Les systèmes qui emmagasinent de l'énergie ne 
peuvent répondre instantanément et présentent des réponses transitoires à chaque fois qu'ils sont soumis à 
des entrées ou perturbations. 


Généralement, le comportement dynamique d'un système peut être entièrement caractérisé par la 
réponse temporelle de ce système suite à une entrée échelon puisqu'elle est facile à générer (Fig 5-6). 


100%. K 


50%. K -- 


ts 


Fig 5-6 : caractéristiques de la réponse transitoire 


La réponse transitoire d’un système suite à une entrée échelon dépend des conditions initiales. Par 
commodité dans la comparaison des réponses transitoires de différents systèmes, il est plus pratique d'utiliser 
les conditions initiales standards (système au repos à l'instant initial et toutes les dérivées par rapport au temps 
sont nulles). Les caractéristiques des réponses peuvent alors être comparées. 


La réponse transitoire des systèmes asservis pratiques présente souvent des oscillations amorties 
avant d'atteindre le régime permanent. Les critères de performances, communément utilisés pour la 
caractérisation des systèmes asservis linéaires dans le domaine temporel, sont définis comme suit : 


e Temps de retard (time delay), ta : il est défini comme étant le temps nécessaire pour que la 
réponse atteigne la moitié de sa valeur finale. 


e Temps de montée (rise time), tr : temps nécessaire à la réponse pour évoluer de 10 à 90%, de 5 
à 95%, ou de 0 à 100% de sa valeur finale. Pour les systèmes du 2" ordre peu amorti, le temps 
de montée de 0 à 100% est plus généralement utilisé. Pour les systèmes très amortis, l’évolution 
de 10 à 90% est plus souvent choisie. 
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e Temps de pic (peak time), tp : temps nécessaire pour atteindre le 1° pic de dépassement. 


e Dépassement maximum, d : c’est la valeur du pic maximal de la réponse mesurée relativement à 
l'unité. Si la valeur finale du régime permanent diffère de l'unité, on utilise plus souvent le 
dépassement maximal exprimé en pourcentage. Il est défini par : 

s(t) — s(%) 
s(%) 
L'importance de ce dépassement maximum (en %) est qu'il renseigne directement sur la relative 

stabilité du système. 


d% = x 100 % (5-3) 


e Temps de réponse ou d'établissement (settling time), ts : c'est le temps requis pour que la 
courbe de sortie atteigne et reste à l'intérieur d'une bande, exprimée en pourcentage 
(généralement 5%), relativement à sa valeur finale. 


Ces 5 grandeurs donnent une mesure directe de la caractéristique transitoire du système asservi 
relativement à sa réponse indicielle. Cela veut dire que le système doit être modifié jusqu’à ce que la réponse 
transitoire soit satisfaisante. 


Il est à noter que ses grandeurs ne sont pas toutes, nécessairement, applicables à n'importe quel 
système. Pour un système très amorti (non oscillant), t et d ne sont pas définis. 


5- 3.1 - Remarques sur les caractéristiques de la réponse transitoire 


Excepté pour certaines applications pour lesquelles les oscillations ne sont pas tolérables, il est 
préférable que la réponse soit suffisamment rapide et suffisamment amortie. Ainsi, pour une réponse transitoire 
acceptable d'un système du 2™ ordre, le coefficient d'amortissement doit avoir une valeur 0.4<¢< 0.8: 


e De faibles valeurs de ¢ (6 < 0.4) produisent un dépassement excessif de la réponse transitoire. 


e Des valeurs importantes de & (£ > 0.8) donnent une réponse très lente. 


Nous verrons, plus loin, que le dépassement maximum et le temps de montée ne peuvent pas être 
faibles tous les deux, simultanément. Si l’un d'eux est diminué, le second croît nécessairement. 


5- 3.2 - Caractéristiques de la réponse transitoire des systèmes du 2" ordre 


Dans la suite, nous allons déterminer le temps de montée, le temps de pic, le dépassement et le temps 
de réponse ou d'établissement d’un système du 2" ordre en fonction de & et de œn. Le système est considéré 
comme étant peu amorti. 


5- 3.2.a - Temps de montée t, (rise time) 


Le système considéré présente une sortie s(t) telle que (voir § 3-6—2) : 


s(t) = K| 1- costa 85 a € e Et (5-4) 
1-6 


Pour t, la sortie vaut: s(t) = K 
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> cos(to,f-#)+-—— sinto,fi-E7)-0 puisque 6% 40 


1-6? 
1- £2 
> toft 1 a 5 
(a 
1 =E 
> tr = arctg (5-5) 
ny1- ë? | $ | 
Avec œp, B et o définis sur la figure 5-7, il vient : 
> tr -L arog 2) > rs ho Ta (5-6) 
Op -0 op  @y1-6? 


Il est clair que pour de faibles valeurs de tr , œp doit être importante. 


Imag 
Pôle , 
r @, | Jp 
Ne cos(f) = ë 
Op = On y1- E? 
B 
-0 0 Réel 
© = Ë On 


Fig 5-7 : définition de œp, B et o 


A partir de l'équation 5—6, la figure 5-8 montre les variations de (@n.tr) en fonction de &. 


On.tr 


Fig 5-8 : Variations de œn.tr en fonction de E. 
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5- 3.2.b - Temps de pic t, (peak time) 
Le temps de pic correspond au maximum de la sortie s(t). II est obtenu en écrivant que la dérivée de 
s(t) par rapport au temps est nulle en ce point. 


DR 1- estoni É stofe e 


1=¢7 


ds(t) _K Eye “*!| cos(op-t) + £ sin(o,.t) |+ eat a. sin(ot) -2 cos(a, t) 
dst) L K| — 2n etat sin(op.t) 
dt 1-6? 
DSU Los 
dt 
1. t = œ (cette solution correspond au maximum uniquement lorsque ¢ > 1) 
2. sin(œp.t) = 0 
sin(œp.t) = 0 > @p.t = N.T n=0, 1, 2, ...... 
tee n=0,1, 2... (5-7) 
6 
Le temps de pic correspond au premier dépassement (n = 1), donc : 
b= — = —— (5-8) 
Op On V1 - Ee 
Smax 
s41 FN A 
Dépassement maximum d 
K 
T, : pseudo période 
On Tp 
el 
On tp @nt 
0 T 2T 3r Ar 
1- E? 1- E? H-E? h-E? 


Fig. 5-9 : Périodicité des maximas et minimas de la réponse indicielle 
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En se référant à l’équ.(5-7) et à la fig. (5-9), les maximas de la réponse transitoire apparaissent aux 
valeurs impaires de n. Les minimas, à celles paires. 


Remarque : 


Il est à noter que bien que la réponse indicielle pour & + 0 ne soit pas périodique, les maximas et les 
minimas de la réponse apparaissent à des intervalles périodiques de période Tp (fig. 5—9). 


5- 3.2.c - Dépassement maximum d (maximum overshoot) 


Ce dépassement apparaît à t = p= © = 2 


d = = 
s(%) K 
= —| cos(tp-@n y1- E°) + T sin(t,.o, y1- E°) e “tp 
= —| cos(r) + É sin(x) le = 
1=E 
d eE > d% = e/V x100% (5-9) 


La figure 5-10 donne les variations du dépassement d (exprimé en pourcentage) relativement au 
coefficient d'amortissement &. 


d% 


100 


Fig. 5-10 : Dépassement d% en fonction du coefficient d'amortissement E. 
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5- 3.2.d - Temps de réponse ou d'établissement t; (settling time) 


La sortie s(t) du système considéré peut être écrite sous la forme suivante (voir § 3-6—2) : 


—Eot _ #2 
s(t) = K| 1- Fe nt +arctg -= i (5-10) 


1-6 


Ont ʻ . Pep : 
Les courbes 1 + f T sont les courbes enveloppes de la réponse indicielle (fig. 5—11). La 
-ő 


sortie demeure toujours à l'intérieur de cette paire d’enveloppes. La constante de temps de ces enveloppes est 
1 


T= ; 
Eon 


Amplitude 


Temps 


Fig. 5-11 : Paire d'enveloppes de la réponse indicielle d'un système du 2" ordre 


La vitesse de décroissance de la réponse transitoire dépend de la constante de temps T. Pour un @n 
donné, le temps d'établissement ts est une fonction du coefficient d'amortissement &. Les systèmes très 
légèrement amortis, présentent un ts plus important que celui pour les systèmes correctement amortis. Pour les 
systèmes très amortis (& > 1), le temps d'établissement devient important à cause du départ très lent de la 
réponse (voir réponses indicielles en fonction de & , § 3-6—2). 


La description analytique exacte du temps d'établissement est difficile à obtenir. Il est cependant 
démontré que, pour un critère de + 5 % et 0 < ¢ < 0.9, ce temps ts varie légèrement et reste approximativement 
égal à 3 fois la constante de temps T. Il atteint un minimum autour de & = 0.68 puis augmente, presque 
linéairement, pour les grandes valeurs de &. 
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Par convention dans la comparaison des réponses transitoires des systèmes, on adopte généralement 
le temps d'établissement suivant : 


b=3T= + 
bOn 


A titre d'information, si le critère adopté est celui de + 2 %, alors : ts = 4.T 


(critère de + 5 %) (5-11) 


Il faudrait noter que ts est inversement proportionnel au produit du coefficient d'amortissement & par la 
pulsation propre non amortie œn. Puisque la valeur de & est généralement déterminée compte-tenu des 
exigences sur le dépassement maximum admissible, ts est déterminé essentiellement par la pulsation propre 
non amortie œn. Cela veut dire que la durée du transitoire peut être variée en ajustant uniquement œn, sans 
modifier le dépassement maximum. 


A partir de cette analyse, il devient évident que pour une réponse rapide, œn doit être important. Pour 


limiter le dépassement maximum d et pour réduire ts, le coefficient d'amortissement é ne doit pas être trop 
faible (fig. 5—10). 


Si 0.4 < é < 0.8, alors 25 % >d > 2.5 % 


5- 3.2.e - Exemple 
Considérons le système asservi à retour unitaire de la figure 5-12 dont les paramètres sont : 
E=0.6 
@n = 5 rd/s 


Fig. 5-12 : Système de 2" ordre 


Déterminer tr, tp, d, et ts lorsque le système est sujet à une entrée échelon. 


op = En y1- E° = 4 rd/s 


o=6.@n=3 
e temps de montée tr : js TE 
Op 
or B = arctg 2 0.93 rd ou encore f = arcos ( & ) = 0.93 rd 
6 
dons pe 0 Lise 
4 
e  tempsdepictp: tp= Ta a = 0.785 s 
Op 4 
e dépassement d : d% = e/V x100% = 9.5 % 
e temps d'établissement ou de réponse ts à 5% près : ts = a nn is 
On o 


Cours dďd'’asservissements linéaires continus (2013-2014) Licence et Master « Automatisme Industriel » (Prot. FELLAH M.K.) 74 


5- 3.3 - Relations Boucle Ouverte - Boucle Fermée à retour unitaire 


5- 3.3.a - Système du 1°’ ordre 


Soit le système asservi à retour unitaire de la fig. 5-13, dont la FTBO(p) = Bo 
1+ TBOP 
Avec : Keo et tBo respectivement le gain et la Constante de temps en boucle ouverte. 
E 
RE A) 
+ 1+ TBoP 
Fig. 5-13: Système du 1°% ordre 
FTBF(p) = Ker 
1+ TBFP 
Avec: Ker et t8F respectivement le gain et la Constante de temps en boucle fermée. 
K ž 
Alors : Keser = — 5> <1 Gain en boucle fermée (5-12) 
1+ KB0 
ter =—BO < teo Constante de temps en boucle fermée (5-13) 
1+ KB0 


La figure 5-14 montre les réponses indicielles en Boucles ouverte et fermée (ex. pour Keo > 1). 


Step Response 


Amplitude 


Time (sec) 


Fig. 5-14 : réponses indicielles en BO et en BF d’un système du 1°% ordre à retour unitaire 
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Remarques : 
e — Le système en boucle fermée est donc plus rapide qu’en boucle ouverte. 
e La pente à l’origine est la même en BO qu'en BF. 


e Un système du 1° ordre en BO (teo, Keo) reste un système du 1°’ ordre en BF (ter, Ker). 


5- 3.3.b - Système du 2" ordre 


Keo 


Soit le système asservi à retour unitaire de la fig. 5-15, dont la FTBO(p) = E i 
1+ 20 p + —— p’ 


©nBO OnBoO 


Avec : Keo, éso, et œngo respectivement le gain, le coefficient d'amortissement, et pulsation propre non 
amortie en boucle ouverte. 


E 
+ 
Fig. 5—15 : Système du 2" ordre 
K 
FTBF(p) = E BF J 
1+2E p+ —— p? 
@nBF OnBF 


Avec: Ker, éer, et ner respectivement le gain, le coefficient d'amortissement, et pulsation propre non 
amortie en boucle fermée. 


K 
Alors : Kg = — 59 <1 Gain en boucle fermée (5-14) 
1+ KB0 
EBF = eo < BO Coefficient d'amortissement en boucle fermée (5-15) 
1 + Kpgo 
OnBF = @hgo 1 + Kgo > @nso Pulsation propre non amortie en boucle fermée (5-16) 


La figure 5—16 montre les réponses indicielles en boucles ouverte et fermée (ex. pour Keo > 1). 
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Step Response 


Amplitude 


Time (sec) 


Fig. 5-16 : réponses indicielles en BO et en BF d’un système du 2°" ordre à retour unitaire 


Remarques : 

e — Le système en boucle fermée est donc plus oscillant qu’en boucle ouverte (sF < &Bo), 

e le système en boucle fermée est plus rapide qu’en boucle ouverte (ta, tp, tr sont plus faibles), 
e le temps d'établissement à 5% près est identique, car le produit (œn.&) reste constant. 


e — Un système du 2ère ordre en BO (neo, Bo, Keo) reste un système du 2ème ordre en BF (œner, EBF, 
Ker). 


5- 3.3.c - Relations d'approximations 
En représentation asymptotique de Bode, le gain d'un système du 2è™e ordre a l'allure suivante : 


Gain (échelle log) 


Fig. 5-17 : Diagramme du 
module (représentation 
asymptotique de Bode) 

œ (échelle log) 
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logK -log1 


Pente = — 2 = 
log 80 — 109 Oco 


Ce qui donne : 


| Ge @nBo Keo (5-17) | 


Remarque : 
Cette relation est approximative puisqu'elle est obtenue à partir d'un diagramme asymptotique. 


Si ce système est à retour unitaire, alors d’après (5-16) : OnBF = Oo Ÿ1+KBo - 


Généralement, K est grand, de sorte que l’on peut écrire : 


OnBF = ®h80 V1+KBo = ORBO Ko Æ Mco 


d'où : 
| OBE + Go 6-189 | 
La marge de phase est : A = @ (@co) + T 
O 
2.Ë8o v 
P(@co) =| -arctg igg: 

Or: ra) | 

1- 

OnBO O=0co 


Oco ~ OnBF — ©ngo y1 T K80 


> E (@co) = -arctg = —arctg 
1-(1+Kao) -Keo 
Si Keo est grand, alors :  (@co) + -arctg 2680 | . arctg- 2.&gr) (Car tgp EO) 
VK8o V1+K80 


>o (co) + - [n -arctgl2.&pF )|= - x +arctg(2.£3r) 


D'où : AQ = -n+arctg(2.6nr)+1 = arctg(2.Epr) 


Conversion en degrés : Aw°. Z- arctg(2.£pr) ~ 2.Ë8r 
D'où : 
Ap° 
z 5-19 
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5- 4 - Effets de l'addition de pôles et de zéros aux fonctions de transfert 


Il a été montré dans les chapitres précédents que le comportement dynamique des systèmes asservis 
dépend énormément de la position des racines de l'équation caractéristique dans le plan complexe (pôles de la 
FTBF). 


Cependant, en pratique, un asservissement réussi ne peut dépendre uniquement du choix des valeurs 
des paramètres du système de façon que les racines de l'équation caractéristique soient correctement placées. 


Nous allons voir que : 


| Bien que les racines de l'équation caractéristique, qui sont les pôles de la FTBF, affectent la 
| réponse transitoire des systèmes asservis linéaires, en particulier la stabilité, les zéros de la fonction 
| de transfert, s'il y en a, sont également importants. 


Ainsi le rajout et /ou la suppression de pôles et de zéros indésirables de la fonction de transfert est 
souvent nécessaire pour obtenir des performances temporelles satisfaisantes. Dans ce chapitre, nous verrons 
que l'addition de pôles et zéros aux fonctions de transfert en boucle fermée peut avoir des effets différents sur 
la réponse transitoire des systèmes bouclés. 


5- 4.1 - Addition d'un pôle supplémentaire à la chaîne directe d'un système asservi à 
retour unitaire 


Pour étudier l'effet de l'addition d’un pôle à la chaîne directe d'un système asservi à retour unitaire, 
considérons le système du second ordre de la figure 5-18, auquel nous rajoutons un pôle (p = — 1 / Tp) 
supplémentaire sur la chaîne directe (fig. 5—19). 


5 p? + 260P + of 


Fig. 5-18 : système asservi de 2"d ordre à retour unitaire 


Tp? + (1+ 260, Tp) P? + 260p + of 


Fig. 5—19 : Rajout d'un pôle supplémentaire à la chaîne directe d’un système asservi de 2" ordre 
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La FTBO devient : 


2 
Gp) = ©n 5-20 
PIE pe 20, jte Top) a 


La FTBF s'écrit alors : 


1+G(P) Tp? +(1+ 2E0, 7, Jp? + 2E0,p + © 


F(p) = FTBF(p) = 


La figure 5-20 illustre les réponses indicielles d'un système en boucle fermée lorsque : 


an = 1 rd/s, 6=1,et Tp = 0, 1,2, et5s 


Step Response 


Amplitude 


Time (sec) 


Fig. 5-20 : Réponse indicielle de F(p) (addition d’un pôle à la chaîne directe). 


Ex : an = 1 rd/s, 6=1, et Tp=0, 1,2, et5s 


Ces réponses montrent que : 


L'addition d'un pôle à la chaîne directe d'un système asservi augmente, généralement, le 
dépassement maximum de la FTBF, ainsi que le temps de montée t, (rise time). 


80 
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5- 4.2 - Addition d'un pôle supplémentaire à la fonction de transfert en boucle fermée 


Considérons la fonction de transfert en boucle fermée d’un système du second ordre auquel on rajoute 
un pôle (p = — 1 / Tp ) supplémentaire (Fig. 5-21) : 


2 
On 


(p? +26o,p + of )(1+ T,p) 


Fig. 5-21 : Rajout d’un pôle supplémentaire à la fonction de transfert d'un système asservi de 2" ordre 


2 
F(p) = FTBF(p) = Pn 5-22 
PE TERRE DE ape oe}ie Top) de 


La figure 5-22 illustre les réponses indicielles du système en boucle fermée lorsque : 
@n = 1 rd/s, 6 = 0.5, et Tp = 0, 0.5, 1,2, et4s 


Step Response 


Amplitude 


Time (sec) 


Fig. 5-22 : Réponse indicielle de F(p) (addition d’un pôle à la FTBF). 
Ex : an = 1 rd/s, 6 = 0.5, et Tp = 0, 0.5, 1, 2, et4s 
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Ces réponses montrent que : 


L'addition d'un pôle à la fonction de transfert en boucle fermée décroît, généralement, le 
dépassement maximum de la FTBF, mais fait augmenter le temps de montée t, (rise time). 


| 
| Si lon ne considère que le dépassement, ajouter un pôle à une fonction de transfert en 
| boucle fermée a un effet opposé à celui obtenu lorsque le pôle est ajouté à la chaîne directe. 


5- 4.3 - Addition d'un zéro supplémentaire à la fonction de transfert en boucle fermée 


Considérons la fonction de transfert en boucle fermée d’un système du second ordre auquel on rajoute 
un zéro (p =— 1 / T2 ) supplémentaire (Fig. 5-23) : 


of(1+T,p) 


(p? + 2Ewp + @f) 


Fig. 5-23 : Rajout d’un zéro supplémentaire à la fonction de transfert d’un système asservi de 2" ordre 


2 
F(p) = FTBF(p) = — (+ TP) 5-23 
(P) = FTBF() = UP (5-23) 


La figure 5-24 montre les réponses indicielles du système en boucle fermée lorsque : 
Avec œn = 1 rd/s, 6 = 0.5, et Tz = 0, 1, 3, 6, et 10s. 


Ces réponses montrent que : 


L'addition d'un zéro à la fonction de transfert en boucle fermée décroît le temps de montée t, 
(rise time) et augmente le dépassement maximum de la réponse indicielle. 
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Step Response 


Amplitude 


Time (sec) 


Fig. 5-24 : Réponse indicielle de F(p) (addition d’un zéro à la FTBF) . 
Ex : on = 1 rd/s, € = 0.5, et T2=0,1,3,6, et 105 


5- 4.4 - Addition d'un zéro supplémentaire à la chaîne directe d'un système à retour 
unitaire 


Considérons qu'un zéro —-1/T; est ajouté à la chaîne directe d'un système asservi de 3è"e ordre de 
FTBO G(p) (Fig. 5-25). 


E 6(1+ T,p) 
5 p? +3p? +(2+6T,)p+6 


Fig. 5-25 : Rajout d'un zéro supplémentaire à la fonction de transfert d'un système asservi de 3ème ordre 


6(1+T,p) 


Gip = -2T R] 
P= pipa 


6(1+T 
FTBF(p) = F(p) = „Â +P) (5-24) 
p'+3p° +(2+6T,)p+6 
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La différence entre ce cas et celui d'ajouter un zéro à la fonction de transfert en boucle fermée est que 
dans le cas présent, le terme (1+T:p) apparaît non seulement au numérateur de F(p) mais le dénominateur de 
F(p) contient également Tz. 


Le terme (1+Tz p) au numérateur de F(p) augmente le dépassement maximum, mais Tz apparaît dans 
le coefficient du terme en p au dénominateur, ce qui a pour effet d'améliorer l'amortissement ou réduire le 
dépassement maximum. 


La figure 5-26, illustre les réponses indicielles lorsque Tz = 0, 0.2, 0.5, 2, 5, et 10 s. 


Step Response 


Amplitude 


Time (sec) 


Fig. 5-26 : Réponse indicielle de F(p) (addition d’un zéro à la chaîne directe d’une FT). 
Ex : T2 = 0, 0.2, 0.5,2,5,et 10 s 


Il est à noter que lorsque Tz = 0, le système est au bord de l'instabilité. Lorsque Tz = 0.2 et 0.5, les 
dépassements maximums sont réduits, principalement à cause de lamélioration de l'amortissement. Lorsque 
Tz croît au delà de 2, bien que l'amortissement continue à être amélioré, le terme (1+Tz p) au numérateur 
devient de plus en plus dominant, et le dépassement maximum devient de plus en plus important au fur et à 
mesure que Tz augmente. 


Une conclusion importante est à tirer de cette discussion : 


Bien que les racines de l'équation caractéristique soient généralement utilisées pour étudier 
le relatif amortissement et la relative stabilité des systèmes asservis linéaires, les zéros de la fonction 
de transfert ne doivent pas être négligés quant à leurs effets sur les performances transitoires du 
système. 


Cours d’asservissements linéaires continus (2013-2014) Licence et Master « Automatisme Industriel » (Prof. FELLAH M.K.) 84 


5- 5 - Pôles dominants des fonctions de transfert 


A partir des discussions précédentes, il devient clair que la position des pôles de la fonction de 
transfert dans le plan de Laplace affecte énormément la réponse transitoire du système. Pour les besoins 
d'analyse et de synthèse, il est important de trier les pôles ayant un effet dominant sur la réponse transitoire. 
On les appellera pôles dominants. 


Puisque la majorité des systèmes de contrôle rencontrés dans la pratique sont d’un ordre supérieur à 
deux, il devient utile d'établir des indications quant à l'approximation des systèmes d’un ordre important par des 
systèmes d'un ordre plus faible aussi longtemps que la réponse transitoire est concernée. 


En synthèse, nous pouvons utiliser les pôles dominants pour contrôler les performances dynamiques 
du système, tandis que les pôles négligeables ou insignifiants sont utilisés afin d'assurer que la fonction de 
transfert du régulateur peut être réalisée par des composants physiques. 


Pour tous les besoins pratiques, nous pouvons sectionner, qualitativement, le plan de Laplace en 
régions dans lesquelles les pôles dominants et les pôles insignifiants sont séparés comme sur la figure 5-27. 
Nous avons, délibérément, choisi de ne pas assigner des valeurs spécifiques aux coordonnées, puisqu'elles 
sont toutes relatives au système considéré. 


Les pôles qui sont proches de l'axe imaginaire du côté gauche du plan complexe 
donnent lieu à des réponses transitoires qui vont s'amortir relativement doucement. 


Les pôles qui se trouvent loin de l'axe (relatif au pôles dominants), correspondent à 
| des amortissements rapides des réponses. 


La distance D entre la région dominante et la région peu signifiante peut être sujet à discussion : il est 
établi en pratique et dans la littérature que si le module de la partie réelle d'un pôle vaut 5 à 10 fois celle d'un 
pôle dominant ou d’une paire de pôles complexes de pôles dominants, le pôle sera considéré comme étant 
négligeable relativement à la réponse transitoire. 


jo 


Plan - p 


Zone des 5 
D pôles H 
= dominants 


E A ET RE: 


Zone de stabilité 4 Zone d'instabilité 


Fig.5-26 : Zones des pôles dominants et insignifiants dans le plan - p 
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Chapitre 6 : LIEU D'EVANS 
(Lieu des racines ou Lieu des pôles) 


6- 1 - Introduction 


Au chapitre précédent, nous avons démontré l'importance des pôles et zéros de la fonction de transfert 
en boucle fermée des systèmes asservis linéaires sur les performances dynamiques du système. 


Les racines de l'équation caractéristique, qui sont les pôles de la FTBF, déterminent la stabilité relative 
et absolue des systèmes linéaires. Cependant, il faut se rappeler que les propriétés transitoires du système 
dépendent également des zéros de la FTBF. 


Pour les systèmes asservis linéaires, un point d'étude important est la recherche et l'investigation des 
trajectoires des racines (lieu des racines) de l'équation caractéristique lorsque certains paramètres du système 
varient. 


Si le système a un gain de boucle variable, la position des pôles de la FTBF dépend de la valeur du 
gain choisi. Il devient, alors, important que le concepteur sache comment se déplacent les pôles de la FTBF, 
dans le plan p, lorsque le gain varie. 


Du point de vue de la synthèse, pour certains systèmes, un simple ajustement du gain peut déplacer 
les pôles de la FTBF vers les positions désirées. Le problème de synthèse devient une simple sélection de 
valeurs appropriées du gain. Si l’ajustement du gain seul ne donne pas les résultats escomptés, il devient 
nécessaire d'ajouter un correcteur (ou compensateur) au système. 


Les pôles de la FTBF sont les racines de l'équation caractéristique. Trouver les racines d'une équation 
caractéristique d’un ordre supérieur à 3 est laborieux et nécessite l’utilisation d’un ordinateur. Cependant, 
trouver les racines de cette équation n’a pas de sens en soi, car lorsque le gain de la FTBO varie, l'équation 
caractéristique change et les calculs doivent être répétés. 


Une méthode simple pour trouver les racines de l'équation caractéristique a été développée par W.R. 
Evans (vers 1950), et utilisée largement dans le domaine de la commande. Cette méthode, appelée la 
méthode du lieu des racines, est l'une de celles qui permettent de tracer les racines de l'équation 
caractéristique pour toutes les valeurs d’un paramètre du système. Les racines correspondant à une valeur 
particulière de ce paramètre peuvent alors être localisées sur le graphe résultant. Il est à noter que ce 
paramètre est, généralement, le gain, mais tout autre paramètre de la FTBO peut être utilisé. Sauf indications, 
nous supposerons que le gain de la FTBO est le paramètre à varier de 0 à l'infini. 


En utilisant la méthode du lieu des racines, le concepteur peut prédire les conséquences, sur la 
position des pôles de la FTBF, de la variation du gain ou de l'addition de pôles et/ou de zéros de la FTBO. 
Cette méthode indique, également, la contribution de chaque pôle et zéro de la FTBO sur la position des pôles 
de la FTBF, et par conséquent, indique de quelle manière ces pôles et zéros peuvent être modifiés pour que 
les performances requises pour le système en BF soient atteintes. Cette méthode est vraiment efficace pour 
trouver rapidement des résultats approximatifs. 
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6- 2 - Propriétés de base du lieu des racines 
6- 2.1 - Condition sur l'angle et condition sur le module 


Considérons le système de la figure 6-1. 


Fig. 6-1 : Système asservi 


S(p) _ Gp) 
E(p) 7 1+ G(p) H) 


L’équation caractéristique est: 1 + G(p).H(p) = 0 > G(p).H(p) =- 1 (6-2) 


La FTBF est: 


(6-1) 


Nous supposerons que G(p).H(p) est un rapport de polynôme en p. Puisque G(p).H(p) est une quantité 


complexe, l'équation (6—2) peut être décomposée en : 


e Une condition d'angle : Z G(p).H(p) =+ 180° (2i + 1) (s0 hass] (6-3) 
e Une condition de module : | G(p).H(p) | = 1 (6—4) 


Les valeurs de p qui satisfont, à la fois, la condition d'angle et la condition de module sont les 
racines de l'équation caractéristique, ou les pôles de la boucle fermée. 


Dans beaucoup de cas, G(p).H(p) fait intervenir le paramètre "gain K", et l'équation caractéristique peut 
être écrite sous la forme : 


 K(P+2:)(p +22)... (P + Zm) 
(p + p1)(P + P2)....(P + Pn) 


=0 (6-5) 


| Le lieu des racines pour le système est alors le lieu des pôles lorsque le gain K est varié de 
| zéro à l'infini. 


Il à noter que pour entamer le tracé du lieu par la méthode du lieu des racines nous devons connaître 
la position des pôles et zéros de G(p).H(p). II faut se rappeler que les angles des quantités complexes résultant 
des pôles et zéros de la boucle ouverte relativement à un point test p sont mesurés dans le sens horaire. Par 
exemple, si G(p).H(p) est donnée par : 


K(p +z) 


G(p).H(p) = 
PHP = bep +) tpa) 


où — p2 et — ps sont des pôles complexes conjugués. 
Alors l'angle de G(p).H(p) est : 
Z G(p).H(p) = db; — 07 — 02 — 03 — 04 
où 67, 07, 02, O3, et 04 sont mesurés dans le sens horaire 


comme indiqué sur les figures 6-2(a) et 6-2(b). 
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Le module de G(p).H(p) pour ce système est : 


KB. 


G(p).H(p) |= T1 — 
GPH I= AAAA, 


où Bı, Aı, A2, A3, et A4 sont les modules des quantités complexes 


(p + z1), (p + p1), (p + p2), (p + p3) et (p + p4) respectivement, 
comme indiqué sur la figure 6-2(a). 


A noter que, du fait que les pôles complexes conjugués et les zéros complexes conjugués, s'il y en a, 
sont toujours localisés symétriquement par rapport à laxe réel, le lieu des racines est toujours symétrique par 
rapport à cet axe. || suffit, par conséquent, de construire la moitié supérieure du lieu, puis de compléter la 
moitié inférieure en projetant le tracé relativement à laxe réel. 


Point test Ajo O: zéro Point test 02 No 
X—— 
A X : pôle j ~ Gfk 
bat SS 2 i 
F Aa 62 -.::7@ 
a Re Ad. Og me A 0; (oj 
f Br: sA pak y ` PL : on 5 
M a de | - pà -Z -pP |o 
/ 101,43 N, 1 
S64 à < z i 
A aa Maod D DIE 
- p4 er a <b |! - P3 
“= 03 
(a) 
x p3 


Fig. 6-2 : Mesure des angles entre des pôles et zéros en BO et un point test p 


6- 2.2 - Exemple 


Avant de présenter la méthode de construction du lieu d'Evans, étudions d’abord l'exemple d'un 
système du second ordre (Fig. 6-3) : 


FTBO(p) = G(p).H(p) = TTE) avec G(p)= bat et H(p) = 1 
HrBrp - SP) -OM `K 
E(p)  1+G(p)H(p)  p?+p+K 


L'équation caractéristique est : p2+p+K=0 


Calculons le lieu des racines de cette équation lorsque le gain K varie de 0 à o. 
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Afin d'illustrer le lieu des racines de l'équation caractéristique pour cet exemple, nous allons, en 
premier lieu, calculer analytiquement les racines de cette équation en fonction de K et, ensuite, faire varier K de 
0 à w. 

Il est à noter que cette manière de faire n’est pas celle usuelle pour le tracé du lieu d'Evans. Celui-ci 
est obtenu en appliquant un certain nombre de règles, présentées ultérieurement. 


Il est évident que si la solution analytique est facile à trouver, l'application de la méthode du lieu 
d'Evans n'est pas nécessaire : 


p2+p+K=0 
: 1 
e Les racines sont réelles pour : K< > p12 = 3 T 7 1— 4K 
' We gl 
e Les racines sont complexes pour š K> > p12 = -a +j VE — 1 


Le lieu des racines, avec comme paramètre le gain K, a la forme de la figure (6—4) 


Fig. 6-4 : Lieu des racines de la FTBF avec K comme paramètre de réglage. 


Le tracé terminé et étant à l'échelle, nous pouvons immédiatement déterminer la valeur de K pour une 


racine (ou pôle) donnée (ex : pour K=17,, alors p1,2 = -2+2 ). 


De cette analyse, il ressort que : 


e Les pôles en BF correspondants à K = 0 sont les mêmes que ceux en BO. 


e Lorsque K augmente de 0 à 4, les pôles se déplacent vers le point (— 2, 0). Dans ce cas, tous les 
pôles (deux pôles) en BF sont réels. Ceci correspond à un système apériodique (la réponse 
temporelle n'est pas oscillatoire). 


e Pour K = 4, les 2 pôles en BF se rejoignent (pôles doubles). Ceci correspond à un système amorti 
critique. 


e Lorsque K dépasse 1⁄4, les pôles en BF deviennent complexes conjugués, quittent laxe réel, et 
bien que leurs parties réelles restent constantes, ces pôles en BF se déplacent le long de la ligne 
p=-"#2. Pour K > 4, le système devient oscillatoire amorti. 


e Si K continue à augmenter, les pôles complexes conjugués tendront vers (— 2 + jæ), le système 
risque, alors, de devenir instable. 
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Il reste à montrer que tous point M (fig. 6-5) sur le lieu d'Evans satisfait la condition sur l'angle. Cette 
condition est donné par : 


p-s =- Zp — Z(p+1)= 180° (2i + 1) (i=0,1,2,...) 
p(p +1) 
= 
or : Zp=86;: et Ip = |AM | 
> 
et: Z (p+1) = 682 et |p+1| = |BM | 
Alors : Z wa =— 0; — 62= — (0; + 02)= —-180° (puisque M se trouve sur la médiane) 


Fig. 6-5 : Angles et modules entre pôles en BO et un point test M. 


e Sile point M se trouve sur l'axe réel entre O et +, alors (8; + 62) = 0°. Cette partie de l’axe réel 
n'appartient pas au lieu. 


e De même, si M se trouve sur l’axe réel entre —œ et —1, alors (8, + 682) = 2 * 180° = 360°. Cette 
partie de laxe réel n'appartient pas au lieu. 


e Si, par contre, M se trouve sur laxe réel entre —1 et 0, alors (8; + 62) = 180°. Cette partie de 
l’axe réel appartient au lieu puisqu'elle satisfait la condition d'angle. 


Il est évident que si un point n'appartient pas au lieu, alors -(8; + 02) ++ 180° (2i + 1). Donc, tous les 
points qui n'appartiennent pas au lieu ne satisfont pas la condition d'angle et ne peuvent pas être les 
pôles en BF, quelque soit la valeur du gain K. 


Si les pôles en BF sont spécifiés sur le lieu, la valeur correspondante du gain est alors déterminée à 
partir de la condition sur le module. Par exemple, si les pôles en BF sélectionnés sont ( — {2 + 2j ), la valeur 
correspondante de K est alors : 


| FTBO(p) | = Ta =1 > K= [P+D] 1, =] (+ 2))(- 2+ 2j + 1) |= H 
2 


p(p +1) 
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appliquées. 


difficultés. 


A partir du lieu des racines, on peut maintenant étudier les effets d'un changement de la valeur de K 
sur le comportement de la réponse transitoire du système de 2° ordre : 


L'augmentation de K fait baisser le coefficient d'amortissement & (& = cos B), ce qui augmente le 
dépassement de réponse. 


Imag 


jop 


Fig. 6-6 : Pôle complexe d'un système du 2" ordre 


L'augmentation de K, fait augmenter la pulsation propre non amortie œn et la pulsation propre 
PR 
amortie œp = œn teg í 


Si K devient supérieur à la valeur critique (celle correspondant à un système amorti critique pour 
lequel & = 1), l'augmentation du gain K n'affecte plus la partie réelle des pôles en BF. 


Du lieu, il est clair que tous les pôles en BF sont situés dans le demi-plan gauche du plan 
complexe. Quelque soit l'augmentation de K, le système reste stable. Un système du 2" ordre est 
donc toujours stable. Cependant, il ne faut pas oublier que cette stabilité n’est que théorique, et 
qu'un système dont la marge de phase est faible est pratiquement instable, et donc inutilisable. 


Par ailleurs, il est à noter que si le gain devient très important, les effets des constantes de temps 
négligées se font sentir, et le système, supposé du 2"4 ordre mais en réalité d'un ordre plus élevé, 
peut devenir instable. 


6- 3 - Règles de tracé du lieu d'Evans 


6- 3.1 - Généralités 


Si l'équation caractéristique est de degré n, elle à n racines, qui peuvent être réelles ou 
complexes. 


Par conséquent, le lieu comprend n branches, chacune d'elles correspondant à une racine de 
l'équation caractéristique. 


Puisque les racines réelles sont représentées par des points de l'axe réel, les branches 
correspondant aux racines réelles seront des portions de l’axe réel. 


Les branches correspondant aux racines imaginaires peuvent être plus ou moins compliquées. 


Puisque les coefficients de 1 + G(p).H(p) sont réels, les racines imaginaires sont groupées par 
couples, et les branches du lieu correspondant aux racines sont symétriques par rapport à l’axe 
réel. Le lieu d'Evans tout entier est, par conséquent, symétrique par rapport à laxe réel. 


6- 3.2 - Règles de tracé du lieu d'Evans 


Pour un système complexe avec plusieurs pôles et zéros de la BO, construire le lieu des racines peut 
paraître compliqué, mais en fait cela n'est guère plus difficile si les règles de tracé sont correctement 


En se basant sur des points et des asymptotes particuliers, et en calculant les angles de départ à partir 
des pôles complexes et les angles d'arrivée au zéros complexes, la construction du lieu peut se faire sans 
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Si on considère le système asservi de la figure (6—1), l'équation caractéristique est : 
1 + G(p).H(p) = 0 


En premier lieu, il faut réarranger cette équation de façon que le paramètre à intérêt apparaisse comme 
facteur multiplicatif (ici le gain K) : 


1 + K(p+2z,)(p+22)....(P+2m) _ÿ (K > 0) (6-6) 
(P +P1)(P +p2)---- (P +Pn) 


1. Localiser les pôles et zéros de G(p).H(p) dans le plan p. Les branches du lieu d’Evans partent 
des pôles en BO et arrivent aux zéros (zéros finis ou zéros à l'infini). 


Les points de départ du lieu (K = 0) sont les pôles de la FTBO(p) = G(p).H(p) 


Les points d'arrivée du lieu (K > œ) sont les zéros de la FTBO. 


Soient : n = nombre de pôles finis de G(p).H(p) 


m = nombre de zéros finis de G(p).H(p) (avec m < n pour un système réel) 


On dit que si m + n, la FTBO comporte alors (n — m) zéros infinis implicites. Si nous comptabilisons 
l'ensemble des zéros finis et infinis : (nombre de pôles en BO = nombre de zéros en BO). 


2. Le nombre de branches infinies (ou directions asymptotiques) est égal à (n — m). 


3. Les branches asymptotiques présentent des déviations oi telles que : 


pe ei (i=0, 12.) (6-7) 


n-m 


i = 0 correspond à l’asymptote avec la plus faible déviation par rapport à laxe réel. Bien que i est 
supposé évoluer jusqu’à œ, langle de déviation se répète, et le nombre de branches est égal à n — m. 


4. Les branches asymptotiques se rejoignent, sur l’axe réel, au point d’abscisse xi tel que : 


ôles de G(p).H - áros de G(p).H 
, - 2 pôles e G(p) Dee e G(p).H(p) 


(6-8) 


5. Les branches de lieu sur l’axe réel sont déterminées par les liens entre les pôles et zéros de 
la BO. Pour cela, positionner un point test entre chaque pôle et zéro sur l'axe réel, puis 
calculer le nombre de pôles et zéros réels à droite de ce point test. Si ce nombre est impair, 
ce point test appartient au lieu. S'il est pair, le point test n'appartient pas au lieu. Par 
conséquent, la partie du lieu sur l’axe réel est formée de segments alternés. 


6. Points de séparation et de rencontre du lieu sur l’axe réel : 


e Les points de séparation sont les points qui correspondent aux valeurs critiques de K pour 
lesquelles deux racines réelles se confondent en une seule (double) avant de se séparer à angles 
droits en deux racines complexes conjuguées 


e Les points de rencontre sont les points qui correspondent aux valeurs critiques de K pour 
lesquelles deux racines complexes conjuguées se confondent en une seule (double) avant de se 
séparer en deux racines réelles. 
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Par conséquent : 


Si le lieu lie 2 pôles adjacents de la FTBO sur laxe réel, il existe alors, au moins, un point de 
séparation entre ces 2 pôles (fig. 6-7a).. 


De même, si le lieu lie 2 zéros adjacents de la FTBO (un des 2 zéros peut se trouver à —œ) sur 
l'axe réel, il existe alors, toujours, au moins un point de rencontre entre ces 2 zéros (fig. 6—-7b). 


Si le lieu lie un pôle de la FTBO et un zéro (fini ou infini) de la FTBO sur laxe réel, soit qu'il 
n'existe ni point de séparation ni point de rencontre, soit qu'ils existent tous les deux (fig. 6-7c). 


K 2 © K > œ% 


(a) (b) 


Fig. 6-7 : Quelques cas de figures de points de séparation et de rencontre sur l'axe réel 


Les points de séparation et de rencontre ont pour abscisses o qui peuvent être calculés : 


e soit par: y TIe =0 


(6-9) 
o -pi G—Zi 


(avec pi, zi respectivement les pôles et zéros de la BO) 


Dans ce cas, la relation (6-9) ne donne que les points de séparation ou de rencontre qui ne sont pas 
issus de pôles ou zéros multiples. 


e soit en écrivant que : 1 + G(p).H(p) = 1 + K.G1i(p).Hi(p) =0 (équation caractéristique) 
1 
> KP) =- ———— 
G (p)-H4 (p) 
puis en calculant : [e =0 (6-10) 
dp Le 
d 1 
c'est à dire : —| —— = 0 
some 


7. Angle de départ (respec. angle d'arrivée) du lieu des racines à partir de pôles complexes 
(respec. zéros complexes) : 


Afin de tracer le lieu avec la meilleure précision possible, nous devons déterminer les directions du 
lieu aux alentours des pôles ou zéros complexes. L’angle de départ (respec. angle d'arrivée) du lieu 
à partir des pôles complexes (respec. sur les zéros complexes) peut être déterminé en substituant 
de 180° la somme algébrique des angles de tous les vecteurs entre tous les autres pôles et zéros, 
et le pôle complexe (respec. zéro complexe) concerné (fig. 6-8). 
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Angle de départ d’un pôle complexe = 
180° — (somme des angles des vecteurs entre le pôle complexe concerné et les autres pôles) 


+ (somme des angles des vecteurs entre le pôle complexe concerné et les zéros) 


Angle d'arrivée d’un zéro complexe = 
180° — (somme des angles des vecteurs entre le zéro complexe concerné et les autres zéros) 


+ (somme des angles des vecteurs entre le zéro complexe concerné et les pôles) 


Angle de départ ^ mag 
À 
À ae. 
0“ : a 
| 02 
= 


Fig. 6-8 : Angle de départ du lieu à 
partir de pôles complexes 


8. Intersection du lieu avec l'axe imaginaire : 


Les points d'intersection du lieu avec l’axe jœ peuvent être déterminés facilement soit en utilisant le 
critère de stabilité de Routh, soit en posant p = jo dans l’équation caractéristique, en égalant les 
parties imaginaire et réelle de cette équation à zéro, et en calculant œ et K. Les couples de valeurs 
(œ, K) trouvés donnent les fréquences et les gains pour lesquelles le lieu coupe l'axe imaginaire. 


9. Détermination des pôles de la Boucle Fermée : 


Tous point du lieu d'Evans est un pôle en BF si la valeur du gain K en ce point vérifie la condition 
sur le module. Inversement la condition sur le module permet de déterminer la valeur du gain Ken 
tous point particulier du lieu. (le lieu peut être gradué en fonction de K, si nécessaire). 


La valeur de K correspondant à n'importe quel point pı du lieu peut être obtenu en utilisant la 
condition sur le module, c’est à dire : 


_ produit des modules des vecteurs entre p4 et les pôles 
produit des modules des vecteurs entre p4 etles zéros 


Cette valeur peut être évaluée soit graphiquement, soit analytiquement. 
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6- 3.3 - Exemple 1 


Soit à déterminer le lieu des racines du système dont la FTBO(p) = S 
p(p +2)(p +5) 

e Points de départ du lieu = pôles de la FTBO p=0, p=-2, p=-5 > n=3 
e Points d'arrivée du lieu = zéros de la FTBO: aucun > m=0 
e — Nombre de branches asymptotiques = n-m = 3 
e Directions asymptotiques 

ga a (i=0, 1, 2, ....) 

n-m 
gs ED (= 0,1, 2...) 


@= +60°, —60 °, 180° 


e Intersection des asymptotes sur l’axe réel 


Ÿ_ pôles de FTBO(p) - $. zéros de FTBO(p) 
Xi = 


n-m 


Root Locus 


+10 , K=70 


Imag Axis 


Real Axis 


Fig. 6-9 : Lieu d'Evans du système de l'exemple 1 
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e Points de séparation (ou de rencontre) o sur le laxe réel : 
En appliquant (6-9), on obtient : 
1 1 1 


' = > oc=—-0,88 et o=-3,79 
o o+2 o+5 


De même, En appliquant (6—10), on obtient : 
K(o) =- o0 (o +2)(o +5) 


Š Klo) =- (80? + 140 +10) =0 >  o=-0,88 et o=-3,79 
(0J 


La valeur o = — 3,79 n'appartient pas au lieu, elle est donc à exclure. 


e Valeur du gain au pointo (Ko): 
La condition sur le module donne la valeur de Ks. 
| FTBO(o) | = 1 > Ks = 4 


e Intersection du lieu avec laxe imaginaire 


Si on ne désire avoir que la valeur du gain lorsque le lieu rencontre laxe imaginaire, on peut appliquer 
tout simplement le critère de stabilité de Routh. 


En appliquant ce critère on trouve que le système est stable pour 0 < K < 70. 


La valeur limite K = 70, correspond donc au point recherché. 


Si on veut calculer également les pôles au point d’intersection, on part de l'équation caractéristique et 
on remplace p par jœ. Les deux équations obtenues (relatives à la partie imaginaire et réelle) donnent les 2 


paramètres K et œ. On trouve : K=70, © =+ 10 > p12=+jv10. 


6- 3.4 - Exemple 2 


Soit à déterminer le lieu des racines du système dont G(p) = ses et H(p)=1 
p +2p+3 
FTBO(p) = PAS 
p +2p+3 
e Points de départ du lieu = pôles de la FTBO: p=-1+j 2, p=-1-j 2 > n=2 
e Points d'arrivée du lieu = zéros de la FTBO: p=-2 > m=i 
e Nombre de branches asymptotiques = n-m= 1 
e Directions asymptotiques 
qe ere (i= 0,1, 2...) 
n-m 
gs SE (i=0, 1, 2, …..) 
E; = 180° 


e Une portion du lieu des racines se trouve sur l'axe réel : (entre — œ et — 2) 
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Imag 
X 
ns © & 
> -I O) Réd 
X 


Fig. 6-10 : Détermination du lieu sur l’axe réel 


e Points de séparation (ou de rencontre) o sur le lieu 
En appliquant (6-9), on obtient : 


to u d ho 
o+1-jV2 o+1+jV2 +2 


A . 2 
z (+2)o+1+j2+0+1-jV2)-(02 +20+3) Lo 


(o? +20 +3)(o +2) 


> 0?+40+1=0 > o=—3,73 et 


En appliquant (6-10), on obtient : 


2 
Ka +26 +3 
o+2 
n2 
> À kio) __ (20+2)(o+2) i +20 +3) -0 
do (o +2) 
> +40+1-0 > o=—3,/3 et 


La valeur o = — 0,27 n'appartient pas au lieu, elle est donc à exclure. 


e Valeur du gain Ksau point o = — 3,73 


La condition sur le module donne la valeur de Ks. | FTBO(o) | = 1 


e Angle de départ à partir des pôles complexes conjugués de la FTBO : 


pôle complexe migre vers l’axe réel ou s'étend vers l’asymptote. 


+ 180° (2i+1), avec i= 0, 1,2, ..... : 
01— O1 +902) = + 180° (2k +1) > 0: = 180° — 02° + þ1’= 180° — 62 +0: 


o =-— 0,27 


o =-— 0,27 


> Ko = 5,46 


La présence de la paire de pôles complexes conjugués de la FTBO requiert la détermination de l'angle 
de départ à partir de ces pôles. La connaissance de cet angle est importante puisque le lieu des racines, aux 
alentours d’un pôle complexe, donne l'information sur la manière avec laquelle le lieu prenant naissance du 


En se référant à la figure 6—11, si nous choisissons un point test et que nous le déplacions dans les 
alentours immédiats du pôle p = — pı complexe de la FTBO, nous trouvons que la somme des contributions 
angulaires du pôle p = — p2 et du zéro p = — Z; vis à vis du point test peut être considéré comme restant 


inchangée. Si le point test appartient au lieu des racines, alors la somme de 01°, -07, et -92° doit être égale à 


(P proche de pı) 
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L’angle de départ est alors : 
0; = 145° 


Point test P 
[a 


í 6: 
(is 


8; = 180° — 90° + arctg ( V2 ) = 180° — 90° + 55° 


également égale à 145°. 


e Tracé du lieu d'Evans (voir fig. 6- 12) 


Root Locus 


Imag Axis 


NS Imag 
ri a | S pi 
LL 1 o 
A =: > 
=Z] i i Réel 
$ 02 
a 02 


Fig. 6-11 : Angle de départ 0;du lieu à partir du pôle complexe 


Real Axis 


Fig. 6-12 : Lieu d'Evans du système de l'exemple 2 


Puisque le lieu est symétrique par rapport à laxe réel, l'angle de départ à partir du pôle p = — p2 est 
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6- 4 - Sensibilité du lieu d'Evans aux variations des Configurations pôles- 
zeros 


Il est à noter que quelques fois un léger changement dans la configuration pôle-zéro peut causer des 
changements significatifs sur l'allure du lieu des racines. La figure 6-13 montre par exemple qu'une 
modification de la valeur du zéro (de p = —3 à gauche à p = -2.5 à droite) donne des allures de lieu 
complètement différentes. 


5 5 
Fig. 6-13 : Sensibilité du lieu 
0 z 0 Ə 
-5 -5 
-8 -6 -4 -2 0 -8 -6 -4 -2 0 


6-5- 


Simplification d'un pôle de G ar un zéro de H 


Il est important de signaler que si le dénominateur de G(p) et le numérateur de H(p) ont des facteurs 
communs, les pôles et zéros correspondants de la FTBO s'élimineront entre eux, réduisant de ce fait le degré 
de léquation caractéristique de 1 ou plus. Considérons, par exemple, le système de la figure 6-14(a). En 
modifiant le schéma fonctionnel de la fig. 6-14(a) en celui de la fig. 6-14(b), puis en celui de la fig. 6-14(c) (les 
2 branches de retour sont en parallèle, leurs effets se rajoutent), il est clair que G(p) et H(p) ont un facteur 
commun (p + 1). 


Fig. 6-14 : Système asservi présentant des réductions pôles-zéros 
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La FTBF(p) = S(p) / E(p) est : 


K 

S(p) — G(p)  p(p+1)}(p+2) K 

EARS FES), MBTI (p+1)[p(p+2)+K] 
p(p+1)(p+2) 


L'équation caractéristique est donc:  [p(p+2) + K](p+1) = 0 


A cause de la réduction des termes (p+1) apparaissant dans G(p) et H(p), nous aurons cependant : 


K(p+1) p(p+2)+K 
1 + G(p}H Dj RD EPS 
ANR + Hp+1P+2 ppa 


L'équation caractéristique réduite est: p(p+2)+K=0 


Le lieu d'Evans de G(p)H(p) ne montrera pas toutes les racines de l'équation caractéristique, mais 
seulement celles de l'équation réduite. 


Pour obtenir l'ensemble des racines de la FTBF, nous devons rajouter le pôle réduit de G(p)H(p) aux 
pôles de la FTBF obtenus à partir de G(p)H(p). 


Ce qu'il faut retenir c’est que le pôle réduit de G(p)H(p) est un pôle en BF du système, comme cela est 
montré sur la fig. V-14(d). 


6- 6 - Quelques aspects importants de la construction du lieu des racines 


Un des aspects importants de la technique du lieu des racines est que pour la majorité des systèmes 
asservis, de complexité modérée, l'analyste ou le constructeur peut obtenir des informations vitales sur les 
performances du système en construisant uniquement un tracé rapide du lieu en utilisant quelques unes de ses 
propriétés. 


Il est, par conséquent, important de comprendre, du point de vue de la synthèse des systèmes, les 
effets sur le lieu de l'addition ou du déplacement de pôles et zéros de la FTBO dans le plan de Laplace. 


6- 6.1 - Effets de l'addition de pôles à G(p).H(p) 


Ajouter un pôle à G(p)H(p) a pour effet de pousser le lieu des racines vers le demi-plan droit du 
plan de Laplace. 


Les figures 6-15 donnent quelques exemples des effets de l'addition de pôles à la FTBO sur le lieu 
d'Evans. Il est clair que les pôles rajoutés peuvent avoir un effet déstabilisant sur le système. 
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Fig. 6-15 : Effets, sur le lieu d'Evans, de l'addition de pôles à la FTBO. 


6- 6.2 - Effets de l'addition de zéros à G(p).H(p) 


Ajouter des zéros à G(p)H(p) a, généralement, 
racines vers le demi-plan gauche du plan de Laplace. 


pour effet de courber et d'attirer le lieu des 


Les figures 6-16 donnent quelques exemples des effets de l'addition de zéros à la FTBO sur le lieu 
d'Evans. Il est à remarquer que la stabilité relative est accrue par l'addition de zéros à la FTBO. 
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Fig. 6-16 : Effets, sur le lieu d'Evans, de l'addition de zéros à la FTBO. 


6- 7 - Conclusion 


chapitre consacré à la correction des systèmes asservis, cette méthode sera largement utilisée. 


performances dynamiques du système. 


2 
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Dans ce chapitre, nous avons introduit la technique du lieu d'Evans pour les systèmes asservis 
linéaires. Cette technique constitue une méthode graphique d'investigation des racines de l'équation 
caractéristique des systèmes linéaires lorsque un ou plusieurs paramètres varient. Au cours du prochain 


Cependant, il faut se rappeler que les racines de l'équation caractéristique donnent une indication 
exacte de la stabilité absolue des systèmes linéaires mais ne donnent qu'une information qualitative sur la 
relative stabilité, puisque les zéros de la FTBF, s'il y en a, jouent également un rôle important relativement aux 
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Chapitre 7 : Correction des systèmes asservis linéaires 


7- 1 - Introduction 
7- 1.1 - Nécessité de correction dans les systèmes asservis 


Nous avons vu, dans les chapitres précédents relatifs à l'analyse des systèmes asservis, que, pour 
satisfaire aux spécifications de stabilité et de précision, on est amené à formuler des conditions sur la FTBO : 


1. Stabilité : le degré de stabilité est défini par : 
e la marge de gain: la stabilité est d'autant meilleure que : 
— le gain de la FTBO est plus faible, 
— donc, que la bande passante en BO est plus faible. 
e la marge de phase : la stabilité est d'autant meilleure que : 


— le déphasage de la FTBO est plus faible. 


2. Précision : son étude se décompose en deux parties : 


e Précision statique : l'annulation de l'erreur en régime permanent nécessite la présence, dans la 
FTBO, d'une ou plusieurs intégrations selon l'entrée canonique imposée. 


e Précision dynamique : elle est d’autant meilleure que le gain de la FTBO est plus élevé, c'est-à- 
dire que la bande passante est plus large. 


Pour simplifier à l'extrême, on peut retenir, en résumé que la précision et la stabilité sont quantifiées, 
dans le diagramme de Bode de la FTBO, de la manière suivante (fig. 7—1) : 


PEEBO œc : pulsation de 
coupure 
Précision 
1} l j 
Basses 
fréquences 
Z FTBO 
Stabilité (0) 


- 180° 
| Hautes 
fréquences 


Fig. 7-1 : Quantification de la stabilité et de la précision sur le diagramme de Bode 


Il devient clair que : 
e pour améliorer la précision, il faut pouvoir augmenter le gain de la FTBO 


e la stabilité diminue si ce même gain devient trop élevé. 
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Il semble donc difficile d'obtenir un système, à la fois, précis (grand gain) et stable (faible gain). Ce 
premier dilemme stabilité-précision impose donc l'emploi de systèmes compensateurs, correcteurs, ou 
encore régulateurs dont le rôle sera de relever le gain dans une certaine zone de fréquence et de la diminuer 
ailleurs. 


Le correcteur ou régulateur va permettre de satisfaire les contraintes suivantes : 
e trouver un compromis entre la stabilité et la précision, 


e si besoin, rendre stable en boucle fermée un système qui serait instable en boucle ouverte, 


e si besoin, et c'est en général le cas, introduire un intégrateur dans la boucle pour obtenir une erreur 
| statique nulle (e = 0). 


Mais, il existe d’autres incompatibilités qui nécessitent également l'emploi de correcteurs : 


Un bon asservissement doit être insensible aux perturbations et, en même temps, il doit répondre 
rapidement aux variations des diverses grandeurs d'entrée. Ces deux conditions sont incompatibles puisque : 


e |a réponse rapide nécessite une large bande passante, 


e tandis que l'insensibilité aux perturbations exige une bande étroite. 


Le rôle des correcteurs, qui peuvent être électriques, mécaniques, ou hydrauliques, est donc de 
déformer le diagramme asymptotique ou la courbe de Nyquist pour leur donner des marges de gain et de 
phase capables d'assurer la stabilité tout en conservant aux basses fréquences un gain suffisamment grand 
pour que la précision soit bonne. De tels "filtres" pourront également supprimer l'influence de certaines 
perturbations sans limiter la bande passante globale. 


7- 1.2 - Stratégie de correction (ou compensation) des systèmes asservis 


Les outils d'analyse étudiés dans les chapitres précédents conduisent tous au même objectif : la 
synthèse de systèmes asservis corrigés ou compensés, et par conséquent, la synthèse de correcteurs (ou 
compensateurs). 


ut)  _—+ Processus commandé ——> sí) 
entrée de Variable commandée 
commande 


Fig. 7—2 : Processus commandé 


En partant de spécifications sur le comportement final d'un système commandé, tel celui représenté 
par le schéma bloc de la figure 7-2, l'établissement d’un système de correction exige le suivi des 3 étapes 
suivantes: 


1. Déterminer ce que devrait réaliser le système et la manière d'y aboutir (spécifications du cahier des 
charges) 


2. Déterminer la configuration du correcteur en relation avec la manière avec laquelle il est connecté 
au système corrigé. 


3. Déterminer les valeurs des paramètres du correcteur de manière à atteindre les objectifs. 


7- 1.2.a - Cahier des charges 


On utilise toujours un cahier des charges pour décrire ce que doit réaliser le système. Ensuite, sur la 
base de ce cahier des charges, on décide de la manière d'y parvenir. 
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Ces spécifications sont uniques pour chaque application et incluent toujours des exigences sur la 
stabilité relative du système, la précision statique (erreur), la réponse transitoire, et sur les caractéristiques de 
la réponse fréquentielle. Pour certaines applications, il peut y avoir des spécifications supplémentaires sur la 
sensibilité aux variations des paramètres (robustesse, insensibilité aux bruits, etc. ...). 


La synthèse des systèmes de correction linéaire peut être réalisée soit dans le domaine temporel, soit 
dans celui fréquentiel. Par exemple, la précision statique est toujours spécifiée pour une entrée échelon, rampe 
ou accélération, et la démarche classique suivie pour répondre aux contraintes imposées se fait dans le 
domaine temporel. D’autres spécifications, telles que le dépassement, le temps de montée, et le temps 
d'établissement sont toutes définies pour une entrée échelon unitaire, et sont, par conséquent, utilisées 
spécifiquement dans le domaine temporel. 


La stabilité relative est également mesurée en terme de marge de phase, marge de gain et résonance. 
Celles-ci sont des spécifications du domaine typiquement fréquentiel et sont prises en charge par les outils tels 
que le diagramme de Bode, les lieux de Nyquist et / ou de Black-Nichols. 


Pour mener à bien la conception du correcteur dans le domaine temporel ou fréquentiel, il est pratique 
de garder à l'esprit que la synthèse dans le domaine temporel se fait généralement avec l’aide du plan-p de 
Laplace et du lieu d'Evans. La synthèse dans le domaine fréquentiel est basée sur la manipulation du gain et 
de la phase de la fonction de transfert de la boucle jusqu'à atteindre les spécifications voulues. 


Voici quelques principes et règles de base : 


e Les pôles complexes conjugués de la FTBF conduisent à une réponse indicielle oscillatoire. Si 
tous les pôles du système sont réels, la réponse indicielle est très amortie. Cependant, les zéros 
de la FTBF peuvent causer des dépassements importants même si le système est très amorti. 


e La réponse d'un système est dominée par les pôles les plus proches de l’origine du plan-p de 
Laplace. Les transitoires dues aux pôles éloignés de l’origine et à gauche du plan-p s'atténuent 
rapidement. 


e Plus éloignés à gauche dans le plan-p sont les pôles dominants du système, plus vite répondra le 
système et plus importante sera la bande passante. 


e Lorsqu'un pôle et un zéro de la fonction de transfert d'un système ont, plus ou moins, tendance à 
se simplifier entre eux, la zone de réponse du système associée à ce pôle aura un faible module. 


e Les spécifications des domaines temporel et fréquentiel sont étroitement liées entre elles. Le 
temps de montée et la bande passante sont inversement proportionnels. La marge de phase, la 
marge de gain, la résonance, et l'amortissement sont inversement proportionnels. 


Nous avons montré que pour les systèmes du second ordre, il existait des relations analytiques 
simples entre certaines de ces spécifications temporelles et fréquentielles. Cependant, pour les systèmes d’un 
ordre plus élevé, la corrélation entre les spécifications dans ces deux domaines est plus difficile à établir. 


7- 1.2.b - Configurations de correction 


En général, la dynamique d’un processus commandé peut être représentée par le schéma fonctionnel 
de la figure 7-2. 


L'objectif est que la variable commandée, représentée par la sortie s(t), ait un comportement désiré sur 
un intervalle de temps donné. Il s’agit alors de déterminer le signal de commande u(t) qui, dans cet intervalle, 
garantisse la sortie s(t) désirée. 


Plusieurs configurations de base sont possibles, et sont différentes les unes des autres selon la 
position relative du correcteur par rapport au système commandé. Une fois la configuration de correction 
choisie, il ne restera plus qu’à calculer les éléments du correcteur pour répondre au cahier des charges. 
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Nous allons, brièvement, passer en revue les configurations de correction les plus utilisées : 


e Correction cascade ou série : On peut réaliser la compensation en insérant, dans une chaîne, 
un correcteur directement en cascade avec les autres éléments (fig. 7—3). 


e(t) e(t) u(t) s(t) 
—> — 
F (me Correcteur Processus 


Fig. 7-3 : Correction en cascade ou série 


e Correction en réaction ou parallèle : On peut placer ces correcteurs en parallèle sur un élément 
d'une chaîne ; dans ce cas, c’est un correcteur en réaction qui constitue alors une boucle 
secondaire (fig. 7—4). 


| i s(t) 
e(t) e(t) u Processus 


Correcteur 


Fig. 7—4 : Correction en réaction ou parallèle 


D’autres configurations sont également très utilisées, telles que : 


e Correction série-réaction : Le correcteur série A est utilisée conjointement avec un correcteur 
en réaction B (fig. 7—5). 


e(t) e(t) u(t) 
- —p us >> Processus 


Correcteur 
B 


s(t) 


Fig. 7—5 : Correction série-réaction 


e Correction en anticipation : 


1. Dans la fig. 7-6-a, le correcteur anticipatif Gca est placé en série avec le système en boucle 
fermée qui dispose lui-même d'un correcteur Gc dans sa chaîne directe. 


2. Dans la fig. 7-6-b, Le correcteur anticipatif Gca est placé en parallèle avec la chaîne directe. 
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Pour la correction en anticipation, le correcteur Gca est à l'extérieur de la boucle système. Il n'affecte 
pas les racines de l'équation caractéristique du système original. Les pôles et les zéros de Gca peuvent être 
choisis de manière à ajouter ou éliminer les pôles ou zéros de la fonction de transfert en boucle fermée. 


s(t) 


t) 

e(t) Correcteur à C u(t) 
orrecteur 

— Goa —+ Ga —+ Processus 


(a) 

Correcteur 

Gca 
+ 
e(t) u(t) s(t) 
Correcteur 
7 — Gc 7 Processus 

(b) 


Fig. 7—6 (a) et (b) : Correction en anticipation 


Les configurations de correction des figures 7—3 et 7—4 sont toutes à un seul degré de liberté en ce 
sens qu'elles ne disposent que d’un seul correcteur dans chaque système, bien que ces correcteurs puissent 
avoir plusieurs paramètres à faire varier. 


L'inconvénient avec la correction à un seul degré de liberté est que les critères de performances 
réalisables sont limités. Par exemple, si les racines de l'équation caractéristique sont sélectionnées de manière 
à produire une certaine valeur d'amortissement, le dépassement obtenu pour une réponse indicielle peut 
encore être excessif à cause des zéros de la fonction de transfert en boucle fermée (voir chapitre 5). 


Les configurations des figures 7-5, 7-6-a, et 7-6-b sont toutes à deux degrés de liberté. 


7- 1.2.c - Correcteurs de base 


Après avoir opter pour une configuration de correction, le concepteur doit choisir le type de correcteur 
qui, une fois les valeurs de ses éléments déterminées, répondra aux spécifications du cahier des charges. Mais 
même là, il existe une multitude de correcteurs disponibles. En pratique, on choisira le plus simple. Plus le 
correcteur est complexe, plus son coût est élevé, moins fiable il est, et plus difficile il est à réaliser. 


Le choix d'un correcteur spécifique pour une application spécifique est toujours basé sur l'expérience 
du concepteur, et quelquefois sur l'intuition. 


Le correcteur choisi, la tâche suivante consiste à déterminer les valeurs de ses paramètres. Ce sont 
les coefficients d’une ou plusieurs fonctions de transfert composant ce correcteur. L'approche de base est 
l'utilisation des outils d'analyse discutés dans le chapitre 5 (Performances des asservissements) pour 
déterminer comment les valeurs de chaque paramètre individuel influent sur le comportement global du 
système, et par conséquent sur ses performances. A partir de ces informations, les paramètres du correcteur 
sont sélectionnés tels que toutes les spécifications soient atteintes. Bien que, dans beaucoup de cas, cette 
procédure donne directement les résultats escomptés, il faut très souvent la répéter plusieurs fois car certains 
paramètres interagissent entre eux, et influent sur le comportement global. Par exemple, une valeur particulière 
d'un paramètre peut être choisie de sorte que le dépassement soit satisfait, mais en essayant de varier la 
valeur d'un autre paramètre pour obtenir le temps de montée exigé, le dépassement n'est plus acceptable. Il 
devient clair que, plus il y a de spécifications, plus il y a de paramètres du correcteur, et plus compliquée 
devient la conception. 
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Les correcteurs industriels les plus utilisés peuvent être classés, selon leurs actions de correction, de la 
manière suivante : 


e Correcteur à action proportionnelle (P) 


e Correcteur à action intégrale (l) 


e Correcteur à actions proportionnelle et intégrale (PI) 


e Correcteur à action dérivée (D) 
e Correcteur à actions proportionnelle et dérivée (PD) 


e Correcteur à actions proportionnelle, intégrale et dérivée (PID) 


La majorité de ces systèmes utilisent l'électricité ou un fluide sous pression tel que l'huile ou l'air 
comme source d'énergie. Ils sont également classés en fonction du type d'énergie utilisée dans l'opération 
(pneumatique, hydraulique, ou électrique) et choisis selon la nature du système à commander (ou à asservir) et 
selon les conditions d'opération (considérations de sécurité, de coût, de fiabilité, de précision, de poids, et de 
dimension). 


7- 2 - Correction en cascade ou série 
7- 2.1 - Principes généraux 


Si, théoriquement, la place du correcteur dans le schéma global de principe importe peu puisqu'il 
modifie globalement la fonction de transfert en boucle ouverte, pratiquement, on ne peut pas le mettre 
n'importe où. 


En effet, à l'extrémité de la chaîne d'action, on rencontre les organes de forte puissance. Donc, pour 
des raisons d'économie, il serait inopportun de placer un correcteur à cet endroit. 


D'autre part, par la chaîne de retour transitent des signaux d’information provenant généralement 
d'organes de mesure. Donc, il sera possible de placer un correcteur dans cette chaîne. Il ne devra cependant 
pas trop modifier la sortie. 


Par conséquent, en général, les correcteurs en cascade se montent soit à la sortie du comparateur 
dans la chaîne directe (avant amplification), soit dans la chaîne de retour. 


Il ne faut pas perdre de vue que les correcteurs que nous allons introduire ne sont que des réalisations 
imparfaites. En effet, compte-tenu des performances à atteindre, on devra faire la synthèse d'un système 
connu mathématiquement. Or, dans bien des cas, sa réalisation pratique est impossible ; on doit se contenter 
de système ayant des caractéristiques approchant celles du modèle dans une zone de fréquence donnée. Pour 
les autres domaines, le comportement du système réel peut-être très différent de celui du modèle. On devra 


donc tenir compte de ces caractéristiques lorsqu'on insérera le correcteur dans la chaîne. 


Pour déterminer le type de correcteur à utiliser et la valeur des paramètres à adopter, on peut utiliser 
plusieurs méthodes : 


soit considérer les réponses temporelles et analyser les performances statiques et dynamiques du 
système avant et après compensation. 


e Soit, à partir de la courbe de Nyquist du système compensé et par comparaison avec celle que 
l’on doit obtenir, en déduire la structure et les paramètres du compensateur, 


e soit procéder de la même manière, mais avec le diagramme asymptotique de Bode, 


e soit utiliser le lieu d'Evans, le correcteur introduisant de nouveaux pôles et racines. 


Dans la majorité des exemples utilisés jusque là, le correcteur a été un simple amplificateur avec un 
gain constant K. Ce type d'action de commande est connu sous le nom de correction proportionnelle, puisque 
le signal de commande u(t) à la sortie du correcteur est simplement proportionnel au signal à son entrée e(t). 


Il est également possible d'utiliser la dérivée ou l'intégrale du signal d'entrée e(t), en addition avec 
l'action proportionnelle. Par conséquent, nous pouvons considérer, plus généralement, le correcteur comme 
étant un ensemble de composants tels que des comparateurs (additionneurs ou soustracteurs), des 
amplificateurs, des atténuateurs, des dérivateurs et des intégrateurs. La tâche du concepteur est alors de 
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déterminer lesquels de ces composants sera utilisé, dans quelles proportions, et la manière avec laquelle ils 
sont connectés. 


Par exemple, l'un des correcteurs les plus connus et les plus utilisés en pratique est le correcteur PID. 
Les actions intégrale et dérivée du PID ont des implications individuelles de performance et leur utilisation 
nécessite une bonne compréhension des effets de chaque élément de base. 


Pour maîtriser rapidement ce correcteur PID, considérons, séparément, chacune des actions P, |, PI, 
D, PD. 


7- 2.2 - Correcteur à action proportionnelle (P) 
7- 2.2.a - Principe 


La relation entre la sortie u(t) et le signal d'erreur e(t) est : 


u(t) = Kp. e(t) (7-1) 
c'est-à-dire, up) = Kp (7-2) 
e(p) 
avec Kp appelé "gain proportionnel" 


Quelques soient le mécanisme et la source d'énergie utilisés, le correcteur proportionnel est 
essentiellement un amplificateur à gain variable. Son schéma fonctionnel est celui de la fig. 7-7. 


e(t) e(t) K u(t) 
— p — 
+ 
Fig. 7-7 : Correction P. 


La figure 7-8 donne la réponse indicielle du correcteur P. 


e(t) u(t) 
échelon unitaire échelon de valeur Kp 
[i jÁ Kp 


Fig. 7-8 : Entrée et sortie du correcteur P. 


La figure 7—9 donne le diagramme de Bode du correcteur P. 
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1 Bode Diagram 


A 
O 
2 10° 
5 
œ 
z= 
-1 
10 
1 
0.5 
© 
o 
z 0 
(o) 
© 
à -0.5 
-1 
-1 0 1 
10 10 10 
Frequency (rad/sec) 
Fig. 7—9 : Diagramme de Bode du Correcteur P. 
(représenté pour Kh=1) 
7- 2.2.b - Effet 


L'action proportionnelle P crée un signal de commande u(t) proportionnel au signal d'erreur e(t). Elle 
agit donc principalement sur le gain du système asservi et permet d'améliorer notablement la précision. 


L'action proportionnelle : 


e entraîne une augmentation du gain, d'où une diminution de l'erreur statique (amélioration de la 
précision), mais 


e augmente la bande passante du système, ce qui 


e améliore la rapidité du système et, 


e augmente l'instabilité du système. 


Le correcteur proportionnel P n’est généralement pas utilisé seul. On verra que tout correcteur possède 
au moins l’action proportionnelle. 


7- 2.2.c - Réalisation pratique 


Une réalisation pratique de ce correcteur en utilisant des circuits passifs et actifs est montrée sur la 


figure 7—10. 
R2 
R4 
R: 
E R3 
U 


Fig. 7—10 : Réalisation du correcteur P avec 2 amplificateurs opérationnels. 
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R, R ; 
U(p) 2 4 2K, avec Kp = R2/R: si  R4= R3 


Le circuit de la figure 7—10 utilise 2 amplificateurs (le second servant d’inverseur, avec un gain de 
valeur 1 en prenant R4 = R3). 


7- 2.2.d - Exemple 
La figure 7—11 montre le schéma fonctionnel d'un exemple de correction Proportionnelle. 


Fig. 7—11 : Exemple de correction P. 


Pour différentes valeurs de Kp : 


e les réponses indicielles s(t) de la FTBF du système corrigé sont reportées sur la figure 7—12, 


les diagrammes de Bode et de Nyquist de la FTBO du système corrigé sont montrées, 
respectivement, sur les figures 7—13 et 7—14. 


Step Response 


Kp=19 (d%=46.4 ; Tr=0.295 s ; Ts=1.72 s) 
Kp=9 (d%=32.7 ; Tr=0.454 s ; Ts=1.77 s) 
Consigne 


Kp=7/3 (d%=10.6 ; Tr=1.04 s ; Ts=2.04 s) 


Amplitude 


Kp=1 (d%=2.84 ; Tr=1.56 s ; Ts=1.83 s) 


Time (sec) 


Fig. 7—12 : Correction P. 


Réponses indicielles pour différentes valeurs de Kp 
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Bode Diagram 


T 
Q 
e 

o 
v 
2 
E 

D 

© 
= 
F AD =180° 
TD 
= A® =28.1° 
8 

œ 
P= 
[a 

-1 0 1 
10 10 10 
Frequency (rad/sec) 
Fig. 7—13 : Correction P. 
Diagrammes de Bode de la FTBO, pour différentes valeurs de Kp 
Nyquist Diagram 

2 

x 

< 

> 

œ 

c 

D 

œ 

£ 


-2 0 2 4 6 8 10 12 14 16 18 
Real Axis 


Fig. 7—14 : Correction P. 
Lieux de Nyquist de la FTBO, pour différentes valeurs de Kp 


On constate que l'augmentation de Kp, entraîne : 


e une amélioration de l'erreur statique, 
e une décroissance du temps de montée, 
e une faible amélioration du temps d'établissement, 


e mais également une diminution de la marge de phase et une augmentation du dépassement 
(augmentation de l'instabilité du système). 
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La figure 7—15 représente le lieu d'Evans du système corrigé. Elle montre que les 2 pôles en boucle 
fermée sont réels jusqu'à Kp = 1/8, et qu’à partir de cette valeur, ils deviennent complexes conjugués, donnant 
lieu à des réponses oscillatoires plus ou moins amorties. 


Root Locus 


Imag Axis 


Real Axis 


Fig. 7—15 : Correction P. 


Lieu d'Evans, avec variation continue du gain. 


7- 2.3 - Correcteur à action intégrale (1) 
7- 2.3.a - Principe 


La relation entre la sortie u(t) et le signal d'erreur e(t) est : 


MO EK, e(t) 7-3) 
dt 
t 
ou encore u(t)=K; [etat (7-4) 
0 
c'est-à-dire, Mer = Ki 1 (7-5) 
ep) p Tp 
avec Ki appelé " gain intégral " 


Ti appelée " constante de temps d'intégration " 


Son schéma fonctionnel est celui de la fig. 7—16. 


Fig. 7—16 : Correction l. 
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La figure 7—17 donne la réponse indicielle du correcteur Intégrale. 


échelon unitaire 


entrée nulle 


Fig. 7-17 : Entrée et sortie du correcteur I. 


La figure 7—18 donne le diagramme de Bode du correcteur. 


Bode Diagram 


i 
10 
2 
O 
J ne 
z 
œ 
>= 
-1 
10 Š 
-89 
-89.5 
D 
v 
Z -90 
8 
œ 
à -90.5 
-91 
-1 0 1 
10 10 10 
Frequency (rad/sec) 
Fig. 7-18 : Diagramme de Bode du Correcteur |. 
(représenté pour Ki=1) 
7- 2.3.b - Effet 


L'intérêt principal de ce correcteur est d'ajouter dans la chaîne de commande une intégration. Nous 
savons que la présence d’une intégration dans la FTBO augmente la classe du système et réduit ou annule, 
selon le type d'entrée, l'erreur statique du système. 


L'action intégrale pure : 
e améliore la précision en réduisant ou annulant l'erreur statique, mais 


e introduit un déphasage de —90° qui risque de déstabiliser le système (diminution de la marge de 
| phase). 


Le correcteur à action exclusivement Intégrale n'est pratiquement jamais utilisé, en raison de sa lenteur 
et de son effet déstabilisant. Il est, en général, associé au correcteur Proportionnel. 
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7- 2.3.c - Réalisation pratique 


Une réalisation pratique de ce correcteur en utilisant des circuits passifs et actifs est montrée sur la 


figure 7-19. 
Fig. 7-19 : Réalisation du correcteur | avec 2 amplificateurs opérationnels. 
Gc(p) = U(p) _ 1 _K avec Ki= + 
e(p) RC p R;C2 


Le circuit de la figure 7—19 utilise 2 amplificateurs (le second servant d'inverseur avec un gain de valeur 
unitaire). 


7- 2.4 - Correcteur à actions proportionnelle et intégrale (PI) 
7- 2.4.a - Principe 


La relation entre la sortie u(t) et le signal d'erreur e(t) est : 


t 
u(t) = Kp.e(t) + Ki. f e(t).dt (7-6) 
0 
… K ; 
c'est-à-dire, MI Kp + Kiin p+ ei (7-7) 
e(p) p P Kp 
ou encore, U(p) = K, 1 pR ORES (7-8) 
e(p) KLP TP 
avec Kb "gain proportionnel", 


Ki "gain intégral", 


Ti= se "constante de temps d'intégration", 
i 


K 
Tn= — = Kp Ti "dosage de corrélation d'intégrale 


Ki 


Son schéma fonctionnel est celui de la fig. 7—20. 
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Fig. 7-20 : Correction PI. 


La figure 7-21 donne la réponse indicielle du correcteur PI. 


échelon unitaire 


Fig. 7-21 : Entrée et sortie du correcteur PI. 


La figure 7—22 donne le diagramme de Bode du correcteur PI. 


Bode Diagram 


Magnitude (abs) 


Phase (deg) 


= - 0 2 
10 10 10 10 10 


Frequency (rad/sec) 


Fig. 7-22 : Diagramme de Bode du Correcteur PI. 
(représenté pour Kh=2, Ki=1) 
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7- 2,4.b - Effet 


La réponse indicielle montre qu'un correcteur PI assure une transmission instantanée du signal 
d'erreur e, suivi d’une intégration de ce signal. 


e Ce correcteur sera utilisé chaque fois qu'une erreur permanente doit être annulée ou minimisée, 
c'est à dire une amélioration de la précision du système. En effet, il introduit une augmentation du 
gain global du système aux basses fréquences. 


e Par ailleurs, le correcteur PI a un effet déstabilisant en raison du pôle à l'origine (déphasage 
supplémentaire entre 0 et -90°). Mais, Le zéro supplémentaire introduit tend à minimiser cette 
instabilité. 

: ss IE ; 
e ll est recommandé de placer le zéro | —— | du correcteur aux basses fréquences de sorte que le 
p 
déphasage supplémentaire introduit par le correcteur n'affecte pas beaucoup le déphasage global 
du système corrigé. Cependant, s'il est très proche de l'origine, son effet sera compensé par le 
pôle à l'origine. 

e K, sera choisi de manière à modifier, éventuellement, la fréquence de coupure du système corrigé 

et donc sa marge de phase. 


; p -K; ai e | 
e Très souvent, le zéro < est choisi de manière à compenser la constante de temps 
K 
p 
dominante du système initial de sorte que la boucle fermée gagne en rapidité. 


Kp et Ki sont tous deux réglables. Ki ajuste l’action intégrale, tandis que Kp affecte à la fois les actions 
intégrale et proportionnelle. 


7- 2.4.c - Réalisations pratiques 


e Une première réalisation de ce correcteur en utilisant des circuits passifs et actifs est montrée sur 
la figure 7—23. 


R2 C2 


U(p) _ R2 1 Kp + Ki avec Kp =R2/R: et Ki | 


e(p) R; R;C:p p R;C: 


L'avantage du circuit de la figure 7-23 est qu'il n'utilise que 2 amplificateurs (le second servant 
d'inverseur avec un gain de valeur 1). Cependant, ce circuit ne permet pas une sélection indépendante de Kp et 
Ki, puisque ceux-ci dépendent tous les deux de R1. 


e Une deuxième réalisation, avec 3 amplificateurs, cette fois-ci, est celle de la fig. 7—24. 
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gipa VE a Pey l A avec Kp =Rə/Ri et Ki=—— 


e(p) R; R;Cip p R;C; 


Pour ce circuit, Kp et Ki peuvent être réglés séparément. Cependant, pour un tel circuit, Ki est 
inversement proportionnel à la valeur du condensateur. Malheureusement, en général, les corrections PI 
exigent de faibles valeurs de Ki. Cela conduit à de larges valeurs de Ci, ce qui n’est pas très pratique. 


7- 2.4.d - Exemple 


La figure 7-25 montre le schéma fonctionnel d'un exemple de correction PI. 


1 


(p +1)(0.5p +1) 


Fig. 7—25 : Exemple de correction PI. 


Nous fixons Kp = 5, et nous varions la valeur de Ki. 


Les figures 7-26, 7—27 et 7—28 donnent les réponses indicielles s(t) pour, respectivement, Ki = 0 à 2.5, 
K=25à5,etK=5 à 15. 


Les figures 7-29, 7-30 et 7-31 donnent, respectivement, le diagramme de Bode, le lieu de Black- 
Nichols, et le lieu de Nyquist de la FTBO corrigée. 
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Step Response (Kp = 5 = constant) 


Ki/Kp=0.5 (d%=19.5 ; Tr=0.701 s ; Ts=2.52 s) 
Ki/Kp=0.3 (d%=12.6 ; Tr=0.745 s ; Ts=2.81 s) 


Ki/Kp=0.2 (d%=9.04 ; Tr=0.776 s ; Ts=5.79 s) 


Ki/Kp=0.1 (d%=5.42 ; Tr=0.817 s ; Ts=13.1 s) 


Amplitude 


Ki/Kp=0 (correction P, d%=22.1 ; Tr=0.646 s ; Ts=1.51 s, £+ 0) 


Time (sec) 

Fig. 7-26 : Exemple de Correction PI. 

Réponses indicielles pour Ki = 0 à 2.5. 
Kp = 5 (constant) 


Step Response (Kp = 5 = constant) 
Ki/Kp=1 (d%=35.1 ; Tr=0.631 s ; Ts=2.52 s) 
Ki/Kp=—0.8 (d%=29.1 ; Tr=0.654 s ; Ts=2.48 s) 
Ki/Kp=0.6 (d%=22.8 ; Tr=0.683 s ; Ts=2.48 s) 


Ki/Kp=0.5 (d%=19.5 ; Tr=0.701 s ; Ts=2.52 s) 


Amplitude 


Time (sec) 


Fig. 7-27 : Exemple de Correction PI. 


Réponses indicielles pour Ki = 2.5 à 5. 
Kp=5 (constant) 
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Step Response (Kp = 5 = constant) 


Ki/Kp=3 (d%=80.3 ; Tr=0.508 s ; Ts=18.1 s) 
Ki/Kp=2 (d%=60.5 ; Tr=0.554 s ; Ts=5.33 s) 
Ki/Kp=1.6 (d%=51.1 ; Tr=0.579 s ; Ts=4.33 s) 


Amplitude 


Time (sec) 
Fig. 7-28 : Exemple de Correction PI. 
Réponses indicielles pour Ki = 5 à 15. 


Kp = 5 (constant) 


Bode Diagram Kp = 5 = constant 


Magnitude (abs) 


Phase (deg) 


10° 10° 10° 10 10 
Frequency (rad/sec) 
Fig. 7-29 : Exemple de Correction PI. 
Diagrammes de Bode de la FTBO, pour différentes valeurs de Ki 


Kp=5 (constant) 
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Nichols Chart Kp = 5 = constant 


Ki/Kp=0.7 (AD =41°) Ki/Kp=0 (AD=55.9°) 


Ki/Kp=1 (AD or 


T 
Fr 
z Ki/Kp=0.1 (AD =53.8°) 
ë Ki/Kp=0.3 (AD=49.5°) 
© Ki/Kp=0.5 (AD=45.2°) 
3 
= 
(0) 
Q 
O 
-180 -135 -90 -45 0 
Open-Loop Phase (deg) 
Fig. 7-30 : Correction PI. 
Lieux de Black de la FTBO, pour différentes valeurs de Ki 
Kp =5 (constant) 
Nyquist Diagram Kp = 5 = constant 

2 

x 

<{ 

> 

œ 

c 

D 

œ 

£ 


Real Axis 


Fig. 7-31 : Correction PI. 
Lieux de Nyquist de la FTBO, pour différentes valeurs de Ki 


Kp=5 (constant) 
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7- 2.5.a - Principe 


c'est-à-dire, 


avec 


On constate que l'augmentation de Ki, entraîne : 
e une amélioration de la précision par annulation de l'erreur statique, 
e une décroissance du temps de montée, 


e mais également une augmentation du temps d'établissement (correction lente), 


P 


7- 2.5 - Correcteur à action dérivée (D) 


La relation entre la sortie u(t) et le signal d'erreur e(t) est : 


u(t) = Ka dett) 

dt 
Up) _ Ka.p = Ta.p 
e(p) 


Ka appelé " gain dérivé", 


Ta appelée " constante de temps de dérivation" 


Son schéma fonctionnel est celui de la fig. 7—22. 


Fig. 7—22 : Correction à action dérivée pure (D) 


(7 —9) 


(7-10) 


e une augmentation du dépassement et une diminution de la marge de phase (augmentation de 
l'instabilité du système). 


Dans le cas de l'exemple étudié, le système initial est composé de 2 pôles (—1 et —2). Le correcteur PI 
a un zéro qu'il est préférable de positionner le plus proche possible du pôle à l’origine. En effet : 


e si le zéro est exactement à l'origine Zo , il compense le pôle qui s’y trouve. Cela revient à 


faire une correction de type proportionnel en faisant varier Kp. L'erreur statique est, alors, non nulle 
pour une entrée échelon. Le système n’a que 2 pôles en boucle fermée. 


e dune manière générale, plus le zéro se déplace vers la gauche du plan complexe, moins 
importante est la marge de phase du système corrigé (l'apport négatif de la phase du correcteur 
est de plus en plus important), plus importante sont les oscillations de la sortie, et plus de temps 
met le système pour s'amortir. 
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La figure 7-21 donne les réponses à un échelon et à une rampe du correcteur D. 


Impulsion infinie et 


échelon unitaire c 
de courte durée 


Fig. 7-21 : Entrée et sortie du correcteur à action dérivée (D) 


La figure 7-22 donne le diagramme de Bode du correcteur D. 


Bode Diagram 


2 
10 

T 
B 
O 
(o 
o 
35 
= 0 
c 10 
D 
œ 
z= 

91 
— 90.5 
D 
o 
TD 
pa 90 
2 
© 
í 89.5 

89 

si 1 
10 10 10 10 
Frequency (rad/sec) 
Fig. 7-22 : Diagramme de Bode du Correcteur D. 
(représenté pour Ka=1) 
7- 2.5.b - Effet 


du signal d'erreur ; il sera d'autant plus actif que la variation de e(t) est rapide. 


La réponse indicielle montre qu'un correcteur à action exclusivement dérivée ne permet pas la 
transmission d’un signal. L'action dérivée ne peut donc être utilisée seule. On fait appel à elle lorsque le signal 
de commande u doit être particulièrement efficace. En effet, ce correcteur permet de faire intervenir la dérivée 
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L'action dérivée pure : 
e améliore la stabilité du système par l'introduction d'un déphasage supplémentaire de + 90° 
(augmentation de la marge de phase), 


e mais fait diminuer la précision du système, 
e et amplifie les bruits de hautes fréquences. 


Le correcteur à action exclusivement dérivée n’est pratiquement jamais utilisé. Il est en général associé 
au correcteur Proportionnel. 


7- 2.5.c - Réalisation pratique 


Une réalisation pratique de ce correcteur en utilisant des circuits passifs et actifs est montrée sur la 
figure 7-28. 


Fig. 7-23 : Réalisation du correcteur D avec 2 amplificateurs opérationnels 


Gc(p) = P = Ra Cap = Kap avec  Ka= Ra Ca 
elp 


Le circuit de la figure 7-23 utilise 2 amplificateurs (le second servant d'inverseur avec un gain de valeur 


7- 2.6 - Correcteur à actions proportionnelle et dérivée (PD) 
7- 2.6.a - Principe 


La relation entre la sortie u(t) et le signal d'erreur e(t) 


u(t) = Kp.e(t) + g Se (7-11) 
dt 
c'est-à-dire, up Kp +Kap (7-12) 
e(p) 
ou encore, UD- K| 1 eij = Kp (1 + Tap) (7 — 13) 
elp) Kp 
avec Kp "gain proportionnel", 


Ka "gain dérivé " 


Ta= Ka » constante de temps de dérivation" 
p 
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Son schéma fonctionnel est celui de la fig. 7-24. 


e(t) e(t) u(t) 
+ 


Fig. 7—24 : Correction PD 


La figure 7-25 donne les réponses à un échelon et à une rampe. 


échelon unitaire 


rampe unitaire 


> 


Fig. 7-25 : Entrée et sortie du correcteur PD 


La figure 7-26 donne le diagramme de Bode du correcteur PD. 


3 Bode Diagram 


a A 
1 
10° 


Magnitude (abs) 
o 


Phase (deg) 


Frequency (rad/sec) 


Fig. 7-26 : Diagramme de Bode du Correcteur PD. 
(représenté pour Kp=2, Ka=1) 
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7- 2.6.b - Effet 


La réponse indicielle montre qu’un correcteur à action proportionnelle et dérivée (PD) assure une 
transmission instantanée du signal d'erreur e(t) augmenté de sa dérivée de(t)/dt. On verra qu'un tel correcteur 
est à utiliser à chaque fois que le système corrigé doit être plus rapide. 


La correction dérivée peut être vue autrement : Puisque de(t)/dt représente la pente de e(t), la 
correction PD est essentiellement une commande anticipative. Connaissant la pente, le correcteur peut 
anticiper la direction de l'erreur et l'utilise pour mieux contrôler le processus. Normalement, dans les systèmes 
linéaires, si la pente de e(t) ou s(t) est due à un échelon d'entrée important, il apparaît un dépassement 
conséquent. La commande dérivée mesure la pente instantanée de e(t), prédit un grand dépassement, et 
effectue la correction nécessaire avant que le dépassement n’apparaisse. 


Intuitivement, la commande dérivée n'affecte l'erreur statique du système que si cette erreur varie avec 
le temps. Si l'erreur statique est constante, sa dérivée, par rapport au temps, est nulle, l'élément de dérivation 
du correcteur ne produit aucune entrée pour le processus. Mais si l'erreur statique croît avec le temps, u(t) est 
de nouveau développé proportionnellement à de(t)/dt, ce qui réduit l'amplitude de l'erreur. L'équation 7-13 
montre également clairement que le correcteur PD n'altère pas la classe du système, l'erreur statique d’un 
système à retour unitaire dépendant directement de cette classe. 


L'intérêt principal de la correction dérivée est son effet stabilisant. En régime dynamique, elle s'oppose 
aux grandes variations de l'erreur (donc aux oscillations), et permet donc de stabiliser le système et d'améliorer 
le temps de réponse. 


Kp et Ta sont tous deux réglables. Ta ajuste l’action dérivée, tandis que Kp affecte à la fois les actions 
dérivée et proportionnelle. 


L’avance de phase produite par ce correcteur peut être utilisée pour améliorer la marge de phase du 
système asservi. Malheureusement, son gain pousse la fréquence de coupure vers les hautes fréquences. 


Ainsi, l’utilisation du correcteur PD consiste à placer la fréquence de cassure du correcteur, œ = Kp / KD, 
de sorte que l'augmentation effective de la marge de phase ait lieu à la nouvelle fréquence de coupure. 


Pour un système donné, il existe tout un domaine de valeurs optimales Kp / Kb pouvant améliorer 
l'amortissement du système. Une autre considération pratique entrant dans la sélection des valeurs de Kr et Ko 
est l'implantation physique du correcteur PD. 


Par ailleurs, compte tenu de ses caractéristiques fréquentielles de filtre passe haut, le correcteur PD 
accroît, dans la majorité des cas, la bande passante du système et réduit le temps de montée de la réponse 
indicielle. L'inconvénient pratique de cet effet filtre passe haut, est l'accentuation des bruits de hautes 
fréquences provenant de l'entrée. 


Récapitulatif des effets de l’action de correction PD : 


e Amélioration de l'amortissement et réduction du dépassement. 
e Réduction du temps de montée et du temps d'établissement. 
e Augmentation de la bande passante. 


e Amélioration de la marge de phase et de la marge de gain. 


e Possibilité d'accentuation des bruits aux hautes fréquences. 
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7- 2.6.c - Réalisations pratiques 


e Une première réalisation de ce correcteur en utilisant des circuits passifs et actifs est montrée sur 
la figure 7—27. 


Fig. 7-27 : Réalisation du correcteur PD avec 2 amplificateurs opérationnels 


Gap) = “P = f2 + Re O1 p = Ke + Kap avec Kp =R2/R et  Ka=P2C: 


L'avantage du circuit de la figure 7-27 est qu'il n'utilise que 2 amplificateurs (le second servant 
d'inverseur avec un gain de valeur 1). Cependant, ce circuit ne permet pas une sélection indépendante de Kb et 
Ka, puisque ceux-ci dépendent tous les deux de R2. D'autre part, une valeur importante de Ka exigerait une 
large valeur de C4, ce qui n’est pas très pratique. 


e Une deuxième réalisation, avec 3 amplificateurs cette fois-ci, est celle de la fig. 7-28. 


Gp) = LP = 2 + Ra Ca p = Kp + Kipavec Kp =R2/R1 et Ka= Ra Ca 


Pour ce circuit, Kp et Ka peuvent être réglés séparément. Les valeurs importantes de Ka sont obtenues 
en jouant sur Ra et en maintenant Ca dans des proportions raisonnables. 
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7- 2.6.d - Exemple 


La figure 7-29 montre le schéma fonctionnel d'un exemple de correction PD. 


1 


Fig. 7-29 : Exemple de correction PD 


Nous fixons Kp = 1. 


Pour différentes valeurs de Ki , la réponse indicielle s(t) est reportée sur la figure 7—30. 


Kp = 10 (constant) 


Consigne 


ur 


nr ne 7 


Sortie corrigée 


Fig. 7-30 : Exemple de Correction PD. 


Réponses indicielles pour différentes valeurs de Kd. 
Kp =10 (constant) 


On constate que l'augmentation de Ka, entraîne : 

e une diminution appréciable du dépassement (stabilisation du système), 
e une amélioration du temps d'établissement, 

e mais également une persistance de l'erreur statique, 


e et peu d'influence sur le temps de montée. 
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7- 2.7 - Correcteur à actions proportionnelle, intégrale et dérivée (PID) 
7- 2.7.a - Principe 


Le correcteur PID combine les actions des 3 correcteurs P, l et D. 


La relation entre la sortie u(t) et le signal d'erreur e(t) est : 


t 
u(i) = Kelt) + Ki e(tidt + Ka AU 
à dt 
c'est-à-dire, U(p) _ Kp + K + Kap 
e(p) p 
K 
ou encore, Ulp) _ “p CE i 
e{p) p (Kp Kp 


avec Kp "gain proportionnel ", 
Ki "gain intégral", 


Ka "gain dérivé " 


K à 
Ta= — "constante de temps de dérivation" 
p 


K 
Ti = — "constante de temps d'intégration" 
K 


Son schéma fonctionnel est celui de la fig. 7-31. 


Fig. 7-81 : Correction à actions proportionnelle, intégrale et dérivée (PID) 


(7-14) 
(7-15) 
(7-16) 
(7-17) 
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La figure 7-32 donne les réponses à un échelon et à une rampe. 


u(t) 


échelon unitaire 


Cai 
+7 P 
Cai 
seulement 
rampe unitaire 


Fig. 7—32 : Entrée et sortie du correcteur à actions proportionnelle, intégrale et dérivée (PID) 


7- 2.7.b - Effet 


La réponse indicielle montre qu'un correcteur PID assure une transmission instantanée du signal 
d'erreur e(t) augmenté de son intégrale et de sa dérivée. Ce correcteur, facile à réaliser, permet d'annuler le 
signal d'erreur statique £» et d’avoir une réponse relativement rapide et bien amortie. En effet, le correcteur 
PID, fait croître la classe du système d’une unité et introduit 2 zéros qui peuvent être utilisés pour améliorer la 
réponse transitoire (eq. 7-16). La méthode du lieu d'Evans peut être mise à profit pour localiser ces zéros dans 
le but de satisfaire un cahier des charges sur les régimes statique et dynamique. 


Nous avons vu qu'un correcteur P (Kp) apporte de la rapidité au système en réduisant le temps de 
montée. il réduit également l'erreur statique, mais ne l’élimine pas. L'action intégrale (Ki) aura pour effet 
d'éliminer l'erreur statique. Elle ramène donc de la précision, mais dégrade la réponse transitoire. L'action 
dérivée (Ka) améliore la stabilité du système, réduit les dépassements et améliore le régime transitoire. 


Les effets de chaque correcteur (Kp, Ki et Ka) sur la réponse en boucle fermée du système sont 
regroupés sur le tableau 7 —1 : 


Temps de montée | Dépassement! Temps d'établissement Erreur statique 
| Si Kp croît Diminue Augmente (Peu de changements) Diminue 
| Si Ki croît Diminue Augmente Augmente Eliminée 
| Si Ka croît | (Peu de changements)|  Diminue | Diminue (Peu de changements) 


Tableau 7-1 : Effets des correcteurs P, I, et D sur les régimes statique et dynamique 
du système en boucle fermée. 


Il est à noter que ces corrélations ne sont pas exactement précises, car Kp, Ki et Ka sont dépendants 
les uns des autres. En fait, le changement de l’une de ces variables peut modifier l'effet de l'autre. Le tableau 
précédent n’est à utiliser que comme référence lorsqu'il s’agit de déterminer les valeurs de Kp, Ki et Ka. 
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Il existe des méthodes analytiques permettant de calculer les composantes du correcteur PID, mais 
elles sont assez complexes et sont peu utilisées. Des méthodes empiriques existent et permettent de faciliter 
amplement la détermination du correcteurs PID (méthode de Ziegler-Nichols, méthode de Chien-Hrones- 
Reswick,....). Certaines d’entre-elles seront étudiées en TP et TD. 


Des méthodes pratiques de réglages permettent d'obtenir de bons résultats. Elles sont basées sur la 
connaissance des effets que procure chaque correcteur sur la réponse du système bouclé (tableau 7-1). Par 
ailleurs, ces méthodes font beaucoup intervenir l'expérience de l'opérateur dans ce domaine. Il ny a pas de 
réglage unique permettant d'atteindre le cahier des charges, mais il est nécessaire de suivre quelques règles 
d'ajustement de ces correcteurs : 


Obtenir la réponse en boucle ouverte et déterminer ce qui est nécessaire d'améliorer. 


2: Rajouter un correcteur P pour améliorer la rapidité du système : modifier K, pour obtenir le temps 
de montée voulu. 


3. Rajouter un correcteur | pour éliminer l'erreur statique : modifier K; pour améliorer les 
performances en régime statique. 


4. Rajouter un correcteur D pour réduire les dépassements et améliorer le temps d'établissement : 
modifier Ka pour améliorer les caractéristiques en régime transitoire. 


5: Ajuster Kp, Ki et Ka jusqu'à obtenir les performances voulues. 


Finalement, il faut se rappeler qu'il n’est pas obligatoire d'insérer les 3 correcteurs dans un même 
système si cela n'est pas nécessaire. Si un correcteur PI donne des performances satisfaisantes pour la sortie, 
il n'est alors pas nécessaire de rajouter un correcteur D au système. Construire le correcteur aussi simplement 
que possible. 


7- 2.7.c - Réalisations pratiques 


e Une première réalisation de ce correcteur, en utilisant des circuits passifs et actifs, est montrée sur 
la figure 7—33. 


R2 C2 


Fig. 7-83 : Réalisation du correcteur PID avec 2 amplificateurs opérationnels 


Ge(p) = U(p) _ R2 (1+R:C:pX1+R2C2p) 
e(p) R; R2C2p 
Gp) = AP) ) - [Bes sBGp=lés sp 
e(p) R; C, R;C2p p 
avec Kp= Ga “a ; Ket et Ka = R2 Ci 
R; C; R;C32 
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L'avantage du circuit de la figure 7-33 est qu'il n'utilise que 2 amplificateurs (le second servant 
d'inverseur avec un gain de valeur 1). Cependant, ce circuit ne permet pas une sélection indépendante de Kp, 
Ki et Ka. 


e Une deuxième réalisation, avec 3 amplificateurs cette fois-ci, est celle de la fig. 7—34. 


Ge(p) = UP) = P2 + L + Ra Cap = Kp + 1 ep 
e(p) R; RiC;p Kp 
avec Kp =R2/R1  , Ka = Ra Ca et Ki= RiC: 


Pour ce circuit, Kp, Ki et Ka peuvent être réglés séparément. 


7- 2.7.d - Exemple 


La figure 7-35 montre le schéma fonctionnel d’un exemple de correction PID. 


(p +1)(0.5p +1) 


Fig. 7—35 : Exemple de correction PID 


Avec G(p) = nn système à corriger et Gc(p) : Correcteur 


(p+1)(0.5p +1) 
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La figure 7-36 donne la réponse indicielle du système seul (FTBO sans correcteur). 


Le temps de montée est de l’ordre de 3.7 s 


FTBO sans correcteur 


1 — 


[®] 
EN 
N 
w% 
A 
o 
[e>] 
N 
œ 
© 
en 
[®] 


Fig. 7-36 : Réponse indicielle de la FTBO sans le correcteur 


La figure 7-37 donne la réponse indicielle du système seul (FTBF sans correcteur et à retour unitaire). 


FTBF Sans correcteur 


Erreur e (t) \ 


Consigne 


Erreur e (t) et Sortie s(t) 


Fig. 7-37 : Réponses indicielles de la FTBF et de l'erreur e (t) sans le correcteur 
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Le système asservi étudié est très lent (le temps de montée est de l’ordre de 1.83 s), et très peu précis 
(l'erreur statique est de 50%). 


A titre d'exemple, nous nous proposons d'améliorer les performances du système étudié, 
conformément au cahier de charge suivant : 


e Temps de montée : moins de 0.15 s 
e Dépassement : moins de 1% 
e Temps d'établissement : moins de 2 s 


e Erreur statique nulle. 


Pour améliorer, tout d’abord, la rapidité du système, nous insérons un correcteur P, et nous faisons 
varier le paramètre Kp. 


Galp) = AP) =K 


elp) 


La figure 7-38 donne la réponse du système corrigé pour différentes valeurs de Kp. 


Correction P (Kp variable) 


Consigne 


Sortie s(t) 


Fig. 7— 38 : Réponses indicielles du système corrigé par un correcteur P 


Le temps de montée s'améliore avec laugmentation de Kp, mais le dépassement croît fortement. 
Pour Kp = 20, le temps de montée est de l’ordre de 0.286 s, le dépassement vaut 47.4 %. 


L'erreur statique a été réduite mais n’a pas été annulée. 


Tout en maintenant cette valeur de Kp = 20, rajoutons, cette fois-ci, un correcteur intégral | pour 
éliminer cette erreur statique de la réponse indicielle. Nous varierons ensuite le paramètre Ki pour voir son effet 
sur les performances du système. 
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La figure 7-39 donne la réponse du système corrigé pour Kp = 20 et différentes valeurs de Ki. 


Correction PI (Kp =20, Ki variable) 


Sortie s(t) 


Fig. 7— 39 : Réponses indicielles du système corrigé par un correcteur PI 


L'erreur statique a été annulée, mais le temps d'établissement a augmenté. 


Le temps de montée s’est légèrement amélioré avec l'augmentation de Ki, mais le dépassement 
continue de croître. 


Pour Kp = 20 et Ki = 30, le temps de montée est de l’ordre de 0.269 s, le dépassement passe à environ 
69.5 %. 


Nous maintenons les valeurs de Kp = 20 et Ki = 30 et nous rajoutons l’action dérivée D pour améliorer 
le dépassement et donner plus de stabilité au système. Nous varierons ensuite le paramètre Ka pour voir son 
effet sur les performances du système. 

1 
G¢c(p) = U(p) = Kp + — + Kap 
e(p) Kip 


La figure 7—40 donne la réponse du système corrigé pour Kp = 20, Ki = 30 et différentes valeurs de Ka. 


Le temps de montée n’a pas beaucoup changé (entre 0.25 et 0.26 s pour Ka entre 1 et 5), mais le 
dépassement s’est nettement écrasé. Le temps d'établissement a également été réduit. 


Si nous nous contentons de Ka = 10, la courbe est très écrasée (dépassement très faible), mais le 
temps de montée est de l’ordre de 0.55 s. On pourrait alors augmenter Kp pour améliorer, de nouveau, ce 
temps de montée. 
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Correction PID (Kp = 20, Ki = 30, Kd variable) 


Sortie s(t) 


0 0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 
Temps (s) 


Fig. 7— 40 : Réponses indicielles du système corrigé par un correcteur PID 


La figure 7—41 donne la réponse du système corrigé pour, Ki = 30, Ka =10 et différentes valeurs de Kp 
à partir de la valeur 20. 


Correction PID (Kp > 20, Ki = 30, Kd = 10) 


Sortie s(t) 


0 0.5 1 1.5 2 2:9 3 


Temps (s) 


Fig. 7- 41 : Réponses indicielles du système corrigé par un correcteur PID 
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Le temps de montée et le temps d'établissement ont encore diminué, mais le dépassement a 
augmenté. 


Si nous nous contentons de Kp = 80, le temps de montée est passé à environ 0.105 s et le 
dépassement à 9.31%. On pourrait alors augmenter Ka pour améliorer ce dernier. 


La figure 7—42 donne la réponse du système corrigé pour, Kp = 80, Ki = 30 et différentes valeurs de Ka 
à partir de la valeur 10. 


Correction PID (Kp = 80, Ki = 30, Kd > 10) 


Sortie s(t) 


0 0.5 1 1.5 2 2.5 3 
Temps (s) 


Fig. 7—42 : Réponses indicielles du système corrigé par un correcteur PID 
Si nous nous contentons de Ka = 25, le temps de montée a peu changé (environ 0.12 s) mais il ny a 
presque plus de dépassement. 


Avec ces dernières valeurs des paramètres du correcteur PID utilisé, nous répondons aux 
spécifications imposées. 
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Annexe À : TRANSFORMEES DE LAPLACE 


L'analyse temporelle des circuits linéaires en régime transitoire nécessite la résolution d'équations 
différentielles. Pour cela, nous allons introduire un outil mathématique puissant, la transformation de Laplace. 


Cette transformation permet d'associer, à toute fonction f(t), une fonction F(p) d'une variable complexe 
p=o+jo. Elle permet de remplacer les opérations analytiques de dérivation et d'intégration par des opérations 
algébriques. Cette propriété facilite la résolution des équations différentielles. 


A.1 - Transformée Directe de Laplace 
A.1.1 - Définition de la Transformée de Laplace 


Par définition, f(t) étant une fonction réelle du temps (nulle pour t < 0), on appelle Transformée de 
Laplace de cette fonction, notée oÊ {f(t}, la fonction de la variable complexe F(p) telle que : 


F(p) = £L {f0} = f e Pt f(tjdt pour t> 0 
0 
avec : f(t) = 0 pour t < 0 
p: complexe indépendant du temps 


F(p) =£ {{(t)} : transformée de Laplace ou image de f(t) 


f(t) : originale ou fonction objet de F(p). 


A.1.1.a - Exemple 1 


La figure A-1 représente la fonction échelon unitaire u(t) ou Heaviside : u(t) = 


u(t) 
1 
Ei t 


Fig. A-1 : Fonction échelon unitaire 


0 pour t< 0 
1 pour t2 0 


L {u(t)} = U(p) = f e Ptu(t)idt = f e Pt dt = B3 
0 


0 0 E 


Cours d'asservissements linéaires continus (2013-2014) Licence et Master « Automatisme Industriel » (Prof. FELLAH M.K.) 138 


A.1.1.b - Exemple 2 
Soit à calculer À {f "(t)} connaissant oÊ {f(t)}. 


Ona: Ê {fÐ} = F(p) = f e Pt f(t)dt 


0 


En utilisant l'intégration par partie, on aura : 


u = f(t) dv = e™ dt 
e Pi 

du = f” (t)dt ÿ= =— Or : (uv =uv+uv => fuv=uv-fuv 
—P 


= Ê {f t)} =p F(p) — t(0) 


Si la condition initiale est nulle (f(0) = 0), alors : 


of {f° )} = p F(p) 


De même, si toutes les conditions initiales sont nulles (f(0) = f ° (0) = f " (0) =... = 0), alors : 


of {f "t)} = p” F(p) 


Dans ce cas là, l'équation différentielle (pour un système linéaire) liant l'entrée e(t) à la sortie s(t), 


s'écrit, en utilisant la transformée de Laplace, : 


anp” S(p) +..... + ap S(p) +a S(p) = bpp“ E(p)+..... +b,p E(p) +b, E(p) 


(p) avec S(p) = À {s(t)} 


A.1.2 - Propriétés Usuelles de la transformée de Laplace 


À.1.2.a - Linéarité 


Si A1 et À2 sont constants, on a : 


| À {au f1 (t) + ho f2 (t)} = m {fr (t)} + A2  {f2 0) 
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o0 Le] 


car : fe” [ha fa (t) + A2 f2 (D} dt = 4 fes (t) dt + A2 f e Pt fo (t) dt 
0 0 


0 


et 


en particulier : 


L {f(t) + j g(t)} = Lj LÉ {g(t)} 
L {k{(t}=k Lf 


À.1.2.b - Dérivation 


(Voir exemple 1, ci-dessus) 


of {f° {1} = p F(p) — f (0) 


pour t <0. 


D'une manière générale, on peut écrire : 


f (0) représente la valeur de f (t) lorsque t — 0O par valeur positive puisque f (t) n’est pas définie 


A.1.2.c - Intégration 


Ona: LIPA 


0 


=>  À{P(t}- 


en} 1P tn, 
> ieor [PoE] + tferra- D 


= LTD) Tan 


LES (t)} F( = 2n-r (4) avec ft —n-1) (0) dln- t) 
ar ne ta CH CL o 
Exemple : a (t)} = p? F(p)-pf(0)-f"(0) 
3 
MR = (t)} = pĉF(p) -p° f(0) - p f '(0) - f” (0) 


t 
Soit à calculer £i fi (t) dt} = £i P (t) }, P (t) désignant une primitive de f (t) pour t > 0. 
0 


}= f e Pt P(t)dt 


En utilisant l'intégration par partie, on aura : 


dv = e™ dt 


(uv) =uv+uv > fuv=uv-fuv 


En général : 


t 
E Pit) 
p 


0 


en supposant nulles toutes les primitives de f(t) quand t — 0 par valeurs positives. 
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A.1.2.d - Changement d'échelle 
Un changement de l'échelle des temps se traduit par le changement de la variable "t > kt ou t > k/t" 


dans la fonction f. Soit à calculer alors : £ if (kt}} connaissant £ {i (t}} 
Ona: oÊ {f (kb)} = f e™! f (kt)dt 
0 
du 
Posons : kt=u 5 du =k dt => dt = T 
Donc : À {f (kt)} = Ife” f(ujdu = 1 e2) 
| k4 k (k 
iA 1 p 
D'où : Lif kd} = — FE 
| où t = £ f2) | 
| De même que : L {f Eo =k Ffkp) | 


A.1.3 - Théorèmes relatifs à la Transformée de Laplace 


À.1.3.a - Théorème du retard - Translation 


Soit à calculer À {f (t — 1)}, c'est-à-dire la transformée de f (t) quand on fait un changement d’origine 


des temps (Fig. A-2). 
(t) i gt) = f(t- 7) | 
0 i 0 ” i 


Fig. A-2 : Changement d'origine des temps 


Ona : L {f ©} = F(p) = f e Pt f (t)dt pour t20 
0 


G(p) = f e™t g(t)dt = f e™! f (t-r)dt pour t>r 
0 0 
Car, par définition : f(t) = 0 pour t<0 
ou encore : f(t— T) = g{t) = 0 pour ter 


Effectuons le changement de variable x=t-7 


G(p) = Je #19 f (xax = e fe? f (xax = e? Éf (o) = e™ Fip) 
0 0 


L{f(t-1}= e™ Lia pour t > + 


Cours d'asservissements linéaires continus (2013-2014) Licence et Master « Automatisme Industriel » (Prof. FELLAH M.K.) 141 


A.1.3.b - Transformée d'une fonction périodique 


Utilisons la propriété précédente pour calculer la transformée de Laplace d’une fonction périodique à 
partir de t = 0. 


Considérons une fonction périodique de période T pour t > 0 et nulle pour t < 0 (Fig. A-3). 


(D) 


0 T 2T 3T 


Fig. A-3 : Fonction périodique de période T 


La fonction f(t) peut être vue comme une somme de fonctions définies chacune sur une période : 


f(t) = falt) + f(t) + falt) +... = DAG) 


La fonction fı(t) se confond avec la fonction f(t) sur la première période [0, T] et est nulle à l'extérieur 
(Fig. A—4) : 


f1(t) 


0 T 2T 3T 


Fig. A—4 : définition de la Fonction fi(t) 
La fonction f2(t) est définie sur la seconde période [T, 2T] (Fig. A-5) : 


f2(t) 


g T 2T 3T 


Fig. A—5 : définition de la Fonction f2(t) 


Elle se déduit de la fonction fı(t) par un décalage d’une période : 


f2(t) = fi(t —T) 
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De même pour fa(t) (Fig. A-6) : fa(t) 


0 T 2T 3T 


Fig. A-6 : définition de la Fonction f(t) 


ce qui peut être généraliser à : 


Nous pouvons donc encore écrire : 


Calculons la transformée de Laplace de cette expression : 


L{(b}= L [Snan] = > of {fi(t-KT)} 
k=0 k=0 


L{r(}= X e®™F(p) où:  F1(p) =L {fi(t)} 


Ms 


r 
IL 
© 


Nous pouvons donc écrire : 


LEO Fip Sent = FO 
| fo) D ES 


1 
car : ——=1+X+X?+.... 
1-x 


A.1.3.c - Théorème de la valeur initiale 
Soit à démonter que : f (0+) = lim {p F(p) } 
p> 


Ona: L {f (D} = F(p) = f e Pt f(t)at 


0 


Lii =É = fe t'at =p F(p) -t (0 


quand : p>, e 0, donc 0= p F(p) -f (0+ 
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| f(0*)= lim {p F(p)} 


p—co 


A.1.3.d - Théorème de la valeur finale 
Soit à démonter que : f (vo) = lim {p F(p) } 
p— 


Le] 


Ona: ef {f (D} = F(p) = f e Pt f (tdt 
0 


Le] 


d — pt g’ = Z + 
Aie (bat = p F(p) — f (0°) 


Si (p—0), Alors ( e” >1). d'où [rat = lim { p F(p) — f (01) } 
p— 
0 
co t 
Or: fi '(t)dt = Iiri f’ (t)dt = jo {f (t) -f (0+) } 
0 
Donc : D p F(p) — f (0*) } = m {f (t) —f(0*) } 
limf (t)= lim {p F(p) } 
t= p—0 


Ce résultat n’est valable que si { p F(p) } na aucun pôle (racine du dénominateur) dans le demi plan 
droit du plan complexe et aucun pôle sur l'axe imaginaire, à l'exception du pôle simple à l'origine. 


A.1.3.e - Théorème de Duhamel (ou de Borel ) - Intégrale de convolution 


t 
Soit à démontrer que : F;(p).F>(p) = £i filt- 1)-fa(ar J= £i fit) * f2(t) } 
0 


* : produit de convolution 


Posons : fz (t) = f fi(t-1).fa(1)dt 
0 


Par définition, on a: Ê {f (t}= f et ffi(t-1)fa(t)dt dt 
0 


o —— 8 


o0 


En inversant l’ordre des intégrations, on peut écrire : Lif (t)} = [ietode filt- neat 
0 


0 
En posant u=t-t, il vient : 
dits (te f f (dr f fi (uje Pt du 
0 0 


00 


Lin w= f fat) e Ptdr [fitue Pau 
0 


0 


D'où : oÊ {fa (D} = F(p)F(p) = F(p)F(p) 
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D'autre part, nous avons : fs (t) = fit-tali) 
0 


t œ 

Qui peut être écrite : fs (t) = f fi(t-7).fo(t)dt + f fit 0).f(r)dr 
0 t 

Et puisque : fi (t-r) =0 pour t<r 


Alors : fs (t) = f fi(t-7).fa(t)dr 


d'où : 


A.1.3.f - Fonction amortie 


Soit f(t) telle que oÊ {f (t)} = F(p), quelle est alors Ê { et f(t) }? 
Ona: Piet tt} f et f(tje Pidt = f f (hee +Ptdt =F (a+p) 
0 0 


| d{e “f(}=F(a+p) 


A.1.3.g - Multiplication par t 


Considérons la dérivée par rapport à p de f (t) : 


d 
TT =— | f(t)je”tdt 
ap J' 


Comme la variable p n'apparaît que dans l’'exponentielle, nous avons : 


dF(p)_f,. der f -pt 
k | f(t) T dt | filt) }e”'at 


Soit encore : 
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A.2 - Transformée Inverse de Laplace 


On peut exprimer la Transformée inverse, en utilisant les intégrales de Fourrier et de Melin-Fourrier. 


Si F(p) est la Transformée de Laplace d’une fonction f(t), on a : 


où c est une constante, appelée abscisse de convergence. 
Cette méthode est difficile à utiliser et on préfère généralement : 


* soit recourir aux tables de Transformées de Laplace. Dans ce cas, F(p) est immédiatement 
reconnaissable dans la table, 


* soit, lorsque la fonction F(p) n'apparaît pas dans la table, décomposer F(p) en fractions partielles 
et écrire F(p) en termes de fonctions simples de p pour lesquels la Transformée de Laplace est 
toujours connue. 


A noter que cette manière simple de trouver la Transformée inverse est basée sur le fait qu'il existe 
une correspondance unique entre la fonction temporelle et sa Transformée inverse de Laplace du fait de la 
continuité de la fonction temporelle. 


Soit F(p) = £ f} 


Si F(p) peut être décomposée en termes distincts : F(p) = Fi(p) + F2(p) + …… + Fn(p) 
et Si les transformées inverses sont disponibles, 
Alors L 1{ F(p)}} =% 1{ Fi(p)} + L 1{ Fa(p)} +... + L1{ Fn(p)} 
= fit) + fo(t) + ...... + fn(t) 
Remarque : 


Dans le domaine de la Théorie du contrôle, F(p) est fréquemment mise sous la forme : 


F(p) = Bip) avec A(p) et B(p) des polynômes en p, 


A(p) 
et degré B(p) < degré A(p) 


Cette méthode ne s'applique que si les racines du polynôme du dénominateur sont connues, 
autrement dit, que si le dénominateur est factorisable : 


_ B(p) K. (P +Z4)(P + Zo)... (P+Zm) 
A(p) (P + P1)(P + P2)-----.-(P + Pn) 


AL 
g 
l 
1 


OÙ P4,P»,.…...,Pn €t Z4, Z2,- Zm peuvent être des quantités réelles ou complexes. Mais pour chaque 
complexe p ou z, il apparaît un complexe conjugué de p ou z, respectivement. 
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A.2.1 - Si F(p) ne contient que des pôles distincts 


F(p) peut, alors, être décomposée en une somme de fractions partielles : 


B(p) a: 2 an 
F(p) = a + .. + 
A(P) P+P1 P+P2 P+Pn 
avec ai = EC po) ai : constante appelée " résidu au pôle p = pi" 
p p=-Pi 
A.2.1.a - Exemple 1 
Trouver la Transformée Inverse de F(p) = _— 2 pôles distincts: p=-1, p=—2 
(p+1)(p +2) 
p+3 
E (p ] =2 
T ' ma i 
(p) = —— => P 
p+1 p+2 -| P+3 p ] Snr 
>= = 
(p+1)(p +2) a 
Fo =- - — 5 f= LE FE) = PE) LUE 
p+1 p+2 p+i p+2 
=> f(t) = 2e ‘-e ? t>0 
> f(t) = e*(2-e*) t20 


Remarque : 


Dans le cas où le degré de B(p) > degré de A(p) dans F(p) = Fe , il faut alors diviser le 
| 


| numérateur par dénominateur, ensuite appliquer la méthode des fractions partielles. 


A.2.1.b - Exemple 2 


3 2 
Soit: G(p) = p? +5p° +9p+7 
(p +1)(p +2) 


En divisant le numérateur par le dénominateur, on obtient : 


G(p) = p+2+ —P+$ 


popa + 2+ 70) 


( voir F(p) dans l'exemple précédent ) 


=> G(p)=p.1 + 2.1 + F(p) 


g(t) = L { G(p) } = Sat + 2al) + L 1{ F(p)} t20 


g(t = L ‘{ G(p)} = Sat + 2.8( + e (2-68) t>0 


avec ô(t) : impulsion unitaire et £i ô(t)}=1 
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A.2.2 - Si F(p) contient des pôles complexes conjugués 


Soient p. et p, les 2 pôles complexes conjugués, alors : 


S E a a RS E 
A(p) (p+pı)(P+Pp2) P+P3 P+pPn 


Avec o4 et az les résidus aux pôles pietp2: (aP +a ) = ET +p4)(p + pe) 


(p) 


p=-p; ou p=-p2 


A.2.2.a - Exemple 1 

5 p+1 
Trouver la Transformée Inverse de F(p) = — 
p(p° +p+1) 


Ona: p?+p+1=0 pour p=-0,5 + j 0,866 


Donc : F(p) = ni P+@» | a 
(p + 0,5 + j 0,866)(p + 0,5 - j0,866) p 
+1 
Avec : (ap +a se E 
p p=-0,5 - j 0,866 
EE E E a 
-0,5 - j 0,866 


a4 (0,25 + j 0,866 — 0,75) + a (-0,5 — j 0,866) = 0,5 — j 0,866 


En égalant les parties réelles et imaginaires des 2 membres de léquation précédente, on obtient : 
Q4 + Co = —1 


y — Ep =—1 d'où : a =—1 et a) =0 


Déterminons ensuite la valeur de a : 


+1 
a= Fpp =] = 
p^ +p+1 = 


Donc F(p) s'écrit : 
"À 


4 
F un pe 
i P p?+p+1 


1 p 
p  (p+0,5+ j0,866)(p +0,5 - j 0,866) 


1 p 
p  (p=+0,5)? + (0,866)? 


1 p+0,5 0,5 
= + 
p  (p=+0,5)?+(0,866)? (p+0,5)? + (0,866)? 
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Ce qui donne pour f(t) (voir table des Transformées de Laplace) : 


f(t) = L 4 F(p)}=1- e %5t cos ( 0,866 t) + de © °° sin (0,866 t) (t20) 


A.2.3 - Si F(p) contient des pôles multiples 


Soit p, le pôle multiple de F(p), r étant l'indice de multiplicité de ce pôle. 


B 
Fo = D ALP) = (P+ P1)" (P + Peur )(P #Prz)e.(P + Pr) 
Alors F(p) s'écrit : 
F(p) = Ep)... 6 nt Ge ar tet bi pii yeaa 
A(p) (p+pı)  (P+p:) (P+p1) PPr  P+pPr2 P+Pn 
avec : ak = EC (k = r+1, 142, n.. , n) 
A(p) P=—PK 
, B(p) r 
et: b, =a REL p-p) | 
"LAR Les, 


(p+p:1)" 


A.2.3.a - Exemple 1 


ofa [BB n°] =“ e+p]  -(p+2) 1-0 
p=-1 


dp| A(p) dp = 
_ 1 | a? | Ep) | -Hie 2) sn > 
ee WP’ "2e (2p+ f 5 = 
. z 1 -f 2 -1 0 1 L _ (+2 -t 
Donc : f(t) =L ‘{ F(p) } =£ red red Pa (t? +1) e (t>0) 
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Annexe B : Table des Principales Transformées de Laplace 
et leurs propriétés 


Table | Propriétés des Transformées de Laplace | 
des Transformées de Laplace | | 
| ft) t>o | F(p)= Ëi fit) } | f (t>0) | F(p)= Ætt} | 
Impulsion unitaire Ter i 
1 F e ”' f(t)dt 
| de | | | J 
Echelon unitaire 1 | 
= | A Â(E) + A h(t) | AA (p) + XP (p) 
u(t) p | | TO | 
ED ne pF(p) — f(0) 
t p | dt | 
| d’ F(t) 
p al | — 5 p°F(p) — pf(0) — f(0) 
pl | dt 
ME (r-n-1) 
1 d” F(t) n a 2n-r d f0) 
p" | | a A DE Pi 
pra dt r=n+1 dt 
eat | 1 ltt t 
e 2 n 
(pta | | | | (€) .dt 
z |J- . 
P n 0 0 0 = 
te’ (p m ar! | ne conditions Pp 
a initiales nulles) 
ce | mes || 
p(p + a 1 
nl) w | vE 
e tt -e e*t / T Y m b) | | | 
PT a)\P t 
nlat a, | | 
mo | 7 | 
P 
i a 
cos(wt) | z? e | 
2pw 
t.sin(wt) 
(p? cu w? ja 
p —w 
t.cos(wt) 
(p? po Y 
at — U 
Baun (p + a) + w? | = f(t) fonction périodique de période T. 
e % cos{wt) P is a : = f(t) fonction définie sur la 1™ période de f(t). 
(p + a) +w F( ) = F(p) 
(n : entier positif) pP) = 1—e PT 
f(0*)= lim {pF (p)} f(œ) = lim {pF (p)} 
p— 00 p—0 


Si les limites existent 
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